AFWl-TR-87-19 


AFWL-TR* 

87-19 


AD- A 197  334 


/i 


A 


DEVELOPMENT  OF  A  RULE-BASED  STRUCTURAL 
DAMAGE  ASSESSMENT  CODE 

S.  J.  Savage,  et  al 


New  Mexico  Engineering  Research  Institute 
University  of  New  Mexico 
Albuquerque,  NM  87131 


April  1988 


Final  Report 

Approved  for  public  release;  distribution  unlimited. 


V 


V 


AIR  FORCE  WEAPONS  LABORATORY 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base,  NM  871 17-6008 


AFWL-TR-87-19 


This  final  report  was  prepared  by  the  New  Mexico  Engineering  Research 
Institute,  University  of  New  Mexico,  under  Contract  F29601-84-C-0080,  Job 
Order  37630340  with  the  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force 
Base,  New  Mexico.  Mr  Brent  Satterthwaite  (NTES)  was  the  Laboratory  Project 
Officer-in-Charge . 

When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  Government-related 
procurement,  the  United  States  Government  incurs  no  responsibility  or  any 
obligation  whatsoever.  The  fact  that  the  Government  may  have  formulated  or 
in  any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not 
to  be  regarded  by  implication,  or  otherwise  in  any  manner  construed  as 
licensing  the  holder  or  any  other  person  or  corporation;  or  as  conveying  any 
rights  or  permission  to  manufacture,  use  or  sell  any  patented  invention  that 
may  in  any  way  be  related  thereto. 

This  report  has  been  authored  by  a  contractor  of  the  United  States 
Government.  The  United  States  Government  retains  a  nonexclusive,  royalty-free 
license  to  publish  or  reproduce  the  material  contained  herein,  or  to  allow 
others  to  do  so,  for  United  States  Government  purposes. 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  and  is 
releasable  to  the  National  Technical  Information  Service  (NTI3) .  At  NTIS,  it 
will  be  available  to  the  general  public,  including  foreign  nations. 

If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing 
list,  or  if  your  organisation  no  longer  employs  the  addressee,  please  notify 
AFWL/NTES,  Kirtland  AFB,  NM  87117-6008  to  help  us  maintain  a  current  mailing 
list. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


BRENT  SATTERTHWAITE 
Project  Officer 


Lt  Col,  USAF 

Chief,  Applications  Branch 


FOR  THE  COMMANDER 

CARL  L.  DAVIDSON 
Colonel,  USAF 

Chief,  Civil  Engrg  Research  Dlv 


DO  NOT  RETURN  COFIES  OF  THIS  REPORT  UNLESS  CONTRACTUAL  OBLIGATIONS  OR  NOTICE 
ON  A  SPECIFIC  DOCUMENT  REQUIRES  THAT  IT  BE  RETURNED. 


UNCLASSIFIED 


pmajgasmaauEaasaiaE^ 


REPORT  DOCUMENTATION  PAGE 

la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

1b.  RESTRICTIVE  MARKINGS 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

2b.  DECLASSIFICATION  /DOWNGRADING  SCHEDULE 

Approved  for  public  release;  distribution 
unlimited. 

4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

NMERI  WALI-8  (11.04) 

AFWL-TR-87-19 

6a.  NAME  OF  PERFORMING  ORGANIZATION 

New  Mexico  Engineering 

Research  Institute 

6b.  OFFICE  SYMBOL 
(If  applkabl*) 

NMERI 

7a.  NAME  OF  MONITORING  ORGANIZATION 

Air  Force  Weapons  Laboratory 

6c.  ADDRESS  (G'ty,  Start,  and  ZIP  Coda) 

Box  25,  University  of  New  Mexico 

Albuquerque,  NJl  87 131 

7b.  ADORESS  (City.  State,  and  ZIP  Cod*) 

Kirtland  Air  Force  Base,  NM  87117-6008 

1  6a.  NAME  C  •  FUNDING /SPONSORING 
!  ORGANIZATION 

8b.  OFFICE  SYMBOL 

Of  applkabl*) 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F29601-84-C-0080 

1 8c  ADORESS  (City.  Start,  a,*  m  Cotie)  I 

10.  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM 


PROJECT 


ELEMENT  NO.  NO. 
647 1 IF  37 


It.  TITLE  (Intkds  Stci'rity  Claxuflcation) 

DEVELOPMENT  OF  A  RULE-BASED  STRUCTURAL  DAMAGE  ASSESSMENT  CODE 


ORK  Ul 

:cessio 

0 


D.  PERSONAL  AUTHOR(S)  ,  ,  , 

Savage,  Steven;  Ross,  Timothy;  Sorensen,  Harold;  Carson,  James;  and  Satterthwaipe,  Brent 


,5jggg%  Jan87  C0UNT 


1)«.  TYPE  OF  REPORT 
Final 


1C.  SUPPLEMENTARY  NOTATION 


IT.  COSATI  COOES  I  II.  SUBJECT  TERMS  (CondtHM  On  revane  if  nauuar y  and  kkntify  by  block  numb*r) 

QR0UP  |  sub-group - 1  Protective  structures  Damage  assessment  Survivability 

Dynamic  shear  test  Fuzzy  set  Expert  opinion 

Structural  damage  Expert  system  HEST  (over) 


1).  ABSTRACT  (Continue  on  reverie  If  nett  nary  anti  identify  By  bkxk  number) 

This  research  effort  uses  expert  system  techniques  in  the  day^ga  aaaeasmant  of  burled  pro¬ 
tective  structures.  An  expert  system'^pproacJ^ was  chosen'ifor  this  study^because  much  of  the 
Information; used  in  the  damage  assessment  process  is  obtained  from  engineering  judgment  and 
therefor^involves  linguistic  interpretation  of  vague  or  uncertain  terminology  and  ideas. 

The  development  of  this  damage  assessment  coda  is  based  on  a  series  of  -W  tests  conducted  in* 
*1981  and -49W-  on  1/4-scale  reinforced  concrete  boxes.  The  box  elements  in  these  testa  were 
subjected  to  simulated  high  explosive  blast  loads  and  thus  sustained  large  amounts  of  damage. 

Both  subjective  and  objective  Information  are  used  in  this  study The  subjective  information 
vae  derived  from  a  questionnaire  sent  to  experts  around  the  country  concerning  the  levels  of 
damage  and  modes  of  response  in  the  experimental  tests.  Because  the  knowledge  obtained  from 
the  experts  contained  ambiguous  and  uncertain  terminology,  the  theory  of  fuzzy  sets  (over) 


20.  DISTRIBUTION  I AVAILABIUTY  OF  ABSTRACT 
CB  UNCLASSIf UPrtJNI  IMITEO  C  SAME  AS  RPT 
22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 
Brent  Satterthvalte 


21.  ABSTRACT  SECURITY  CLASSIFICATION 
□  otic  USERS  UNCLASSIFIED _ _ ^ 

22b.  TELEPHONE  (Jnc/ytfe  Area  Cotie)  22c  OFFICE  SYMBOL 
(505)  846-6674  NTES 


DO  FORM  1473,14  MAR 


B3  APR  edition  may  be  utad  until  axhauued. 
AJJ  other  edition*  am  obaotete . 

i 


SECURITY  CLASSIC 
UNCLASSIFIED 


CONTENTS 


Section  Page 

I  INTRODUCTION  1 

PROBLEM  STATEMENT  1 

OBJECTIVES  AND  SCOPE  5 

SUMMARY  6 

II  THE  BURIED  BOX  STRUCTURE  7 

EXPERIMENTAL  TEST  SERIES  BACKGROUND  7 

TEST  DESCRIPTION  AND  CONFIGURATION  8 

AVAILABLE  EXPERIMENTAL  DATA  13 

III  ANALYSIS  OF  SUBJECTIVE  AND  OBJECTIVE  INFORMATION  15 

OVERVIEW  OF  DAMAGE  TO  CONCRETE  BOXES  15 

PARAMETERS  AFFECTING  DYNAMIC  RESPONSE  16 

SUBJECTIVE  INFORMATION  20 

Obtaining  Expert  Opinion  20 

Analysis  of  Expert  Opinion  22 

OBJECTIVE  INFORMATION  28 

Damage  Modes  31 

Damage  Level s  36 

IV  THE  STRUCTURE  OF  DAPS  42 

A  RULE-BASED  EXPERT  SYSTEM  APPROACH  42 

UNCERTAINTY  AND  THE  USE  OF  FUZZY  SETS  45 

Fuzzy  Set  Theory  as  an  Approximate  Reasoning  Tool  45 
Fuzzy  Sets  Defined  by  Example  46 

Multiple  Attribute  Decision  Analysis  48 

ORGANIZATIONAL  STRUCTURE  50 

Overview  50 

Structural  Integrity  Analysis  Using  Soft  Data  55 

Explanation  of  External  Programs  Used  in  Soft 
Data  Analysis  63 

Structural  Integrity  Analysis  Using  Hard  Data  67 

Early  Time  Analysis  68 

Late  Time  Analysis  70 

Combination  of  Early-  and  late-Time  Analyses  72 

Summary  72 

V  AN  EXAMPLE  SESSION  75 


iii 


CONTENTS  (Concluded) 

Section  Page 

VI  REVIEW,  CONCLUSIONS,  AND  RECOMMENDATIONS  94 

REFERENCES  99 

APPENDIXES 

A.  PHOTOGRAPHS  OF  OAMAGED  STRUCTURES  101 

8.  EXPLANATION  OF  EXSYS  FEATURES  115 

C.  RULES  119 

LIST  OF  SYMBOLS  153 


iv 


ILLUSTRATIONS 


Figure  Page 

1  Shallow  burled  protective  structure  2 

2  Linguistic  descriptions  of  various  levels  of  damage  4 

3  Dynamic  shear  test  configuration  0 

4  Test  element  construction  details,  Test  Groups  II  and  III  11 

5  Test  Instrumentation  layout.  Test  Groups  II  and  III  14 

fi  Dynamic  strain  rate  effects  on  strength  of  concrete, 

steel,  and  aluminum  17 

7  Typical  failure  curves  with  Group  III  type 

properties  and  fixed  ends  IQ 

8  Possible  sources  of  soft  and  hard  data  for  use  In 

damage  assessment  analysis  21 

8  Reducing  expert  opinion  to  lower  levels  of  knowledge  2Q 

10  Interface  pressure  plot  of  typical  flexural  failure  response  32 

11  Interface  pressure  plot  of  typical  shear  failure  response  33 

12  Deflection  profiles  of  Test  DS3  35 

13  Deflection  profiles  of  Test  0S2-5  35 

14  Typical  expert  system  configuration  43 

15  Classical  and  fuzzy  representations  of  a  moderate  strength 

concrete  4« 

Ifi  Representation  of  structural  damage  assessment  goal  52 

17  Organizational  structure  of  DAPS  53 

18  Logical  flow  of  structural  Integrity  module  54 

18  Matrix  relationships  for  subjective  Information  57 

20  Fxample  fuzzy  set  representing  "No  Relationship"  between 

a  shear  mode  and  structural  attribute  At  58 

21  Fuzzy  SPts  representing  linguistic  damaqe  levels  GO 

22  Calculation  of  difference  measure  81 

23  Organization  of  soft  data  external  functions  84 

24  Organization  of  hard  data  external  functions  8fi 

25  Questions  in  Program  DEFLECT  71 

28  Examples  rules  for  the  combination  of  early,  and  late-time 

analyses  73 

27  Title  block  screen  7fi 

28  Introduction  to  DAPS  78 


v 


ILLUSTRATIONS  (CONCLUDED) 


Figure  Page 

29  Determination  of  failure  77 

30  Introduction  to  FUZZSET  77 

31  Graphical  representation  of  linguistic  terms  79 

32  Fuzzy  set  representation  of  “Moderate  Damage"  80 

33  Introduction  to  subjective  questions  80 

34  Typical  format  of  structural  attribute  questions  81 

35  Preview  of  view/change  option  81 

36  Typical  fuzzy  set  relation  of  mode  versus  structural 

attribute  82 

37  Exiting  from  the  view/change  option  83 

38  Typical  format  of  an  EXSYS  question  83 

39  Results  of  Program  COMBIN  as  shown  by  COLORPLT  84 

40  Typical  rule  explaining  damage  related  to  a  given  mode  86 

41  Typical  rule  explaining  overall  damage  86 

42  Introduction  to  Program  TIM  87 

43  Introduction  to  Program  IFPLOT  87 

44  Interface  pressure  versus  time  at  three  locations  88 

45  Impulse  versus  time  plots  at  three  locations  88 

46  IFPLOT  question  concerning  interface  pressure  plot  similarity  89 

47  Introduction  to  Program  DEFLECT  89 

48  Typical  deflection  profile  at  3  and  15  ms  90 

49  Question  Number  2  in  Program  DEFLECT  90 

50  Typical  rule  interpreting  late-time  response  91 

51  Typical  rule  Interpreting  overall  damage  to  the  structure  as 

determined  by  hard  data  91 

52  Output  of  DAPS  for  Test  0S3  93 


t 


vi 


TABLES 


Table  Page 

1  DYNAMIC  SHEAR  TEST  MATERIAL  PROPERTIES  10 

2  DYNAMIC  SHEAR  TEST  GEOMETRICAL  PROPERTIES  10 

3  DYNAMIC  SHEAR  TEST  LOAD  PARAMETERS  12 

4  EXPERT  REASONING  GIVEN  FOR  SLIGHT  DAMAGE  24 

5  EXPERT  REASONING  GIVEN  FOR  MODERATE  DAMAGE  24 

6  EXPERT  REASONING  GIVEN  FOR  SEVERE  DAMAGE  25 

7  EXPERT  REASONING  GIVEN  FOR  SHEAR  DAMAGE  26 

8  EXPERT  REASONING  GIVEN  FOR  FLEXURAL  DAMAGE  27 

9  STRUCTURAL  ATTRIBUTES  DEVELOPED  FROM  EXPERT  RESPONSES  30 

10  CALCULATION  OF  VALUES  FOR  A  TYPICAL  DAMAGE  LEVEL 

INDICATOR  39 

11  EXAMPLE  RATINGS  AND  WEIGHTS  FOR  DECISION  ANALYSIS  49 


vll 


CONVERSION  TABLE 


To  convert  from 

To 

Multiply  by 

inch  (In) 

millimeter  (mm) 

25.4 

foot  (ft) 

meter  (m) 

0.3048 

pound  (lb) 

kilogram  (kg) 

0.453  592 

pound/square  inch 
(lb/in2,ps1) 

kilopascal  (kPA) 

6.894  76 

kip/sguare  inch 
(k/in%  ksi) 

megapascal  (MPa) 

6.894  76 

viii 


I.  INTRODUCTION 


PROBLEM  STATEMENT 

The  military  forces  of  the  United  States  use  a  great  nuaber  of  shallow- 
buried  reinforced  concrete  protective  structures  for  a  variety  of  purposes. 

As  such,  information  used  in  estimating  or  deteraining  the  survivability  and 
vulnerability  (S/V)  of  these  structures  is  needed  to  evaluate  structural 
safety  and  reliability.  The  quantity  of  data  available  on  which  to  base  a 
study  of  the  responses  of  shallow-buried  structures  to  blast  effects  is  very 
liaited.  Thus,  the  investigation  into  the  S/V  of  these  buried  structures  has 
depended  priaarily  on  conventional  explosive  airblast  simulation  techniques  m 
order  to  develop  a  data  base. 

Figure  1  is  a  scaled  representation,  of  the  shallow  buried  structure 
problem.  Again,  the  issue  in  this  problem  is  the  ability  to  estimate  the 
survivability  of  the  burled  box  to  the  conditions  imposed  by  an  extreme 
environment,  such  as  those  induced  by  an  impulsive  blast  load.  If  we  take  a 
closer  look  at  the  S/V  problem  definition,  we  see  that  the  idea  of  S/V 
analysis  involves  the  questioni  How  well  will  the  structure  continue  to 
function  (survive)  in  a  blast  environment?  By  studying  this  question,  we 
begin  to  envision  the  difficulty  of  the  analysis  process.  For  instance,  how 
does  one  go  about  applying  engineering  quantities  to  such  terms  as 
'functionality  remains  good,'  or  'poor  blast  resistance,*  which  obviously 
involve  degrees  or  levels  of  inforaation  that  are  subject  to  opinion. 

The  problem  is  complicated  by  the  fact  that  the  information  needed  to 
make  an  estimate  with  high  confidence  is  incomplete  and  involves  many 
uncertainties.  Evaluation  of  the  situation  becomes  even  more  difficult  when 
one  realizes  that  the  uncertainties  encountered  include  both  random  and 
nonrandom  kinds  of  data.  For  Instance,  very  little  is  known  about  the  effects 
of  modern  weapons  on  an  operational  facility  which  is  occupied  by  people,  and 
which  contains  sensitive,  technologically  complex  equipment.  Uncertainties 
such  as  aebiguitiet  in  explosive  effects,  variability  in  material  behavior, 
and  threat  scenario  could  be  classified  as  random  uncertainties.  On  the  other 
hand,  linguistic  data,  subjective  judgment,  and  imprecise  information  are 
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Figure  1.  Shallow  buried  protective  structure. 
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typical  examples  of  nonrandom  uncertainties  for  which  the  sample  space  is  not 
well  defined  and  the  mean  and  variance  are  not  meaningful. 

From  the  previous  discussion,  we  conclude  that  the  assessment  of  damage 
to  any  structure,  from  any  extreme  disturbance  (blast/shock,  wind,  earthquake, 
etc...)  is  a  difficult  process  in  which  human  judgment  plays  an  important 
role.  This  is  especially  true  of  protective  structures.  They  are  usually 
heavily  reinforced  and  yet  may  be  subjected  to  very  severe  if  not  total 
failure  levels.  Therefore,  because  of  its  inherent  dissimilarity  with  typical 
civil-engineering-type  problems,  we  cannot  rely  on  past  experience  in  the  area 
of  damage  assessment  with  any  confidence. 

Because  of  a  lack  of  complete  understanding  of  the  real  problem,  in  the 
past  the  typical  analysis  of  the  damaged  structure  would  simply  be  assigned  to 
one  of  two  groups  —  survival  or  failure.  If  we  take  a  closer  look  at  the 
problem,  however,  we  see  that  it  is  not  a  two-class  problem,  but  a  continuous 
one.  To  illustrate,  a  damaged  element  in  a  protective  structure  is  shown  in 
Figure  2.  The  evaluation  for  slight,  moderate,  and  severe  damage  differs 
among  experts)  but,  more  importantly,  there  is  overlap  between  different 
levels.  It  is  this  lack  of  crispness  (or  inherent  fuzziness)  in  the  problem 
that  causes  difficulty;  first,  in  determining  the  damage  level,  and  second,  in 
deciding  on  an  acceptable  level  of  damage. 

Even  when  we  assume  several  gray  levels  between  no  damage  and  total 
failure,  there  exists  possibilities  for  mi sinterpretati on.  The  reason  for 
this  is  that  the  process  of  damage  assessment  is  a  cause  and  effect  analysis. 
The  "cause"  usually  involves  engineering  quantities  (overpressures,  etc,..), 
while  the  “effect"  involves  less  crisp  subjective  information  in  the  areas  of 
functionality  and  repairabi 1 i ty  of  the  structure.  In  other  words,  no  matter 
what  the  cause  of  the  damage,  the  resulting  damage  level  chosen  is  susceotlble 
to  varying  expert  opinion  and  is  highly  dependent  on  the  particular  observer. 
Subjective  information  of  this  type  resides  in  the  minds  of  experts  who  have, 
over  time,  accumulated  much  experience  as  to  the  effects  a  given  loading  would 
produce.  And,  although  there  is  objective  information  available  to  the  expert 
(engineer  or  scientist)  in  the  foru  of  test  data  and  model  simulations,  the 
question  of  survivability  resides  in  the  minds  of  the  experts.  Thus  the 
quality  of  the  assessment  process  is  highly  dependent  on  an  expert’s  knowledge 
of  the  actual  situation  under  study. 
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VERY  EXTENSIVE  DAMAGE 

Figure  2.  Linguistic  descriptions  of  various  levels  of 
damage  (modified  from  Ref.  1). 


OBJECTIVES  AND  SCOPE 


The  discussion  in  the  previous  section  was  provided  in  order  to  introduce 
the  reader  to  the  overall  problem  of  S/V  assessment.  From  this  discussion,  it 
can  be  concluded  that  the  process  of  damage  assessment  can  be  divided  into  two 
subsets:  damage  descriptors  and  damage  levels.  Damage  descriptors  may  be 
either  numerical  engineering  quantities  or  subjective  linguistic  information; 
in  either  case,  they  are  used  to  describe  the  second  subset.  Damage  levels, 
on  the  other  hand,  involve  such  concepts  as  structural  integrity, 
functionality,  and  repairability  of  the  structure  and  the  degrees  to  which 
they  have  been  affected.  The  purpose  of  this  report  is  to  study  the 
feasibility  of  incorporating  the  information  outlined  above  into  a  damage 
assessment  code  in  the  form  of  a  rule-based  expert  system.  The  title  of  this 
expert  system  has  been  designated  DAPS,  which  stands  for  Damage  Assessment  of 
Protective  Structures. 

More  specifically,  the  objective  of  this  report  is  to  develop  a  framework 
by  which  the  concepts  and  information  developed  previously  may  be  converted  to 
a  Norkinq  computer  code.  This  framework  is  to  include  the  coding  of 
information  whether  it  be  in  the  for*  of  expert  opinion,  engineering  judgment, 
digital  wavefores,  etc...,  and  to  develop  a  method  for  the  efficient  retrieval 
of  this  information  for  use  In  an  expert  system.  The  expert  system  will 
combine  numeric  as  well  as  nonnumeric  inforaation  and  will  employ  the  recent 
theory  of  fuuy  set  logic  to  quantify,  combine,  and  interpret  linguistic 
damage  descriptors. 

In  a  broader  sense,  this  work  is  part  of  a  much  larger  project  in  which 
the  Air  Force  is  attempting  to  improve  the  methods  and  procedures  used  within 
the  field  of  information  management.  In  particular,  the  Air  Force  is 
interested  in  improving  five  specific  areas  related  to  struct jral  dynamic 
tests:  tost  design,  data  acquisition,  data  processing,  data  analysis,  and  data 
storage.  The  scope  of  the  work  in  this  report  falls  into  the  category  of  data 
analysis.  Uses  envisioned  for  the  expert  system  which  is  developed  include 
the  following:  design,  prediction,  diagnosis,  and  interpretation,  as  well  as  a 
training  aid  for  engineers  new  to  the  field. 
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SUMMARY 


In  order  to  make  this  oroject  manageable,  the  scope  of  work  was  narrowea 
to  include  only  one  specific  case  of  the  protective  structure  scenario.  In 
particular,  the  data  base  for  this  study  comes  from  a  series  of  eleven 
experimental  tests  which  were  performed  on  buried  reinforced  concrete  boxes 
subjected  to  explosive  pressures.  The  knowledge  base  for  this  work  comprises 
cisp  dava  in  the  form  of  instrumentation  waveforms  and  linguistic  data 
obtained  f^om  experts  through  questionnaires  based  on  the  ‘eleven  tests. 
Measured  data  come  from  active  or  passive  measurements  and  can  be  considered 
as  hard  (objective)  data.  Visual  images  and  linguistic  assessments  from 
expert  opinion  form  what  will  be  called  the  soft  (subjective)  data.  The  term 
"soft  data"  is  chosen  since  it  is  subject  to  individual  judgment  and  not 
readily  quantifiable  (even  though  it  may  contain  substantial  information 
pertaining  to  the  prohlem). 

In  Section  1!  the  eleven  tests  representing  the  knowledge  base  will  be 
reviewed.  The  objective  and  subjective  information  outlined  above  will  be 
thoroughly  examined  in  Socti-'  III  along  with  the  procedures  used  to  obtain 
and  evaluate  them.  Ir.  Section  IV,  introductory  material  on  expert  systems  and 
funy  logic  will  be  discussed.  Th  procedures  used  to  aggregate  and  combine 
the  information  desc-ibta  in  Section  II;  will  also  be  presented,  alonq  with 
the  structure  of  the  BAPS  code.  Section  V  contains  an  example  session  using 
DAPS.  Finally.  Section  VI  contains  a  summary  and  the  conclusions  of  the  work 
performed  }n  association  with  this  report,  ’long  with  recommendations  for 
future  work. 
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II.  THE  BURIED  BOX  STRUCTURE 


EXPERIMENTAL  TEST  SERIES  BACKGROUND 

Experimental  work  on  reinforced  concrete  structural  elements  subjected  to 
short-duration  impulsive  loading  is  limited  (Ref.  2).  The  lack  of 
experimental  data  and  the  need  for  improving  current  S/V  analysis  procedures 
of  shallow  buried  structures  (SBS)  led  to  the  creation  of  the  SBS  Test  Program 
at  the  US  Army  Engineers  Waterways  Experiment  Station  (WES)  in  Vicksburg, 
Mississippi  in  1978. 

The  SBS  program  was  initiated  at  WES  in  order  to  determine  the  magnitude 
of  the  overpressure  that  would  cause  collapse  of  a  structure.  Since  1978,  the 
UES  has  tested  89  one-way  slab  and  box  structures  subjected  to  dynamic  loads. 
All  test  programs  were  associated  with  the  investigation  of  the  behavior  of 
flat-roofed,  shallow  buried  structures  subjected  to  short  duration  (less  than 
3  ms),  high  overpressure  (greater  than  2000  lb/in2)  loads.  By  compiling  the 
information  obtained  from  these  tests  (design  parameters,  analytical  results, 
and  test  results),  a  common  structured  data  base  could  be  established  for  use 
in  characteriiing  concrete  slab  behavior  under  severe  loading  conditions. 

Early  tests  demonstrated  that  the  ability  of  these  structures  to  resist 
highly  impulsive  loads  was  normally  controlled  by  the  response  of  the  roof 
slab,  and  that,  typically,  the  slab  failed  in  a  flexural  mode.  More 
interestingly,  however,  it  was  found  that  direct  shear  failures  could  be 
induced  in  the  structure  by  a  sufficiently  high  overpressure  which  had  a  very 
short  rise  time  and  duration. 

Observations  of  structural  failures  resulting  from  the  Foam  REST  (High 
Explosive  Simulation  Technique)  Test  Series  (Ref.  3)  established  the  need  to 
develop  dynamic  shear  failure  criteria.  Thus,  a  series  of  eleven  tests, 
sponsored  by  the  Department  of  Defense,  was  conducted  in  1981  and  1982  on  one¬ 
way  reinforced  concrete  box  structures.  The  objective  of  these  experiments 
was  to  investigate  possible  shear  failures  resulting  from  highly  impulsive, 
uniformly  distributed  pressures.  Of  particular  interest  in  these  tests  was 
the  response  of  the  roof  element  ( length-to-depth  ratios  between  7  and  10)  of 
the  box  structure  near  the  walls  ( i . e . ,  regions  of  high  shear  stress).  Also 
of  interest  in  these  tests  was  the  ductility  of  the  structure  associated  with 
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dynamic  shear  failures  and  the  ability  to  document  the  response  of  the 
structure  with  high-speed  photography. 

TEST  DESCRIPTION  AND  CONFIGURATION 

The  eleven  dynamic  shear  tests  Mere  conducted  using  1/4-scale  reinforced 
concrete  box  structures  with  no  end  walls.  The  configuration  for  the  tests 
conducted  at  the  WE5  is  shown  in  Figure  3.  As  can  be  seen  in  Figure  3b,  each 
of  the  test  elements  was  placed  in  a  buried  reaction  structure  and  then 
covered  with  a  shallow  layer  of  soil.  A  foam  HEST  charge  cavity  was  placed 
between  the  soil  berm  and  the  soil  covering  the  concrete  structure.  When 
detonated,  this  HEST  charge  provided  the  high-intensity  blast  pressure  to  the 
test  specimen. 

The  eleven  tests  are  categorized  into  three  groups  as  shown  in  Tables  1 
and  2.  The  first  group  is  composed  of  the  five  tests  conducted  in  1981,  while 
the  second  and  third  groups  are  composed  of  the  six  tests  conducted  in  1982. 
Each  of  the  tests  was  identified  by  a  consistent  system  of  notation.  The  1981 
tests  were  given  the  designations  DS1  thru  DSS,  while  the  1982  tests  were 
designated  DS2-1  thru  DS2-6.  As  all  information  and  data  relating  to  these 
tests  is  labeleo  in  this  manner,  this  same  notation  is  adopted  in  this  report. 

All  elements  in  the  test  series  were  similarly  designed  and  constructed. 
Dimensions  and  reinf orcement  patterns  were  the  same  in  all  roofs,  walls,  and 
floors  and  are  shown  in  Figure  4.  The  slabs  had  equal  percentages  of  tension 
and  compression  reinforcement  and  contained  closely  spaced  stirrups 
throughout. 

The  major  variations  among  the  groups  of  tests  were  the  span-to-thickness 
ratio  and  the  reinforcing  ratio.  Properties  that  did  not  vary  within  test 
groups  included  overall  dimensions,  fabrication  scheme,  soil  cover  depth, 
design  concrete  strength,  and  design  steel  strength.  Haterial  and  geometrical 
properties  for  each  of  the  tests  are  listed  in  Tables  1  and  2,  respectively. 
Note  that  the  concrete  strengths  were  denoted  as  averages  of  cylinder 
strengths  at  or  near  the  date  of  testing. 

Although  the  design  load  (charge  density)  was  the  same  for  some  of  the 
tests,  variations  among  the  tests  included  peak  pressure  along  the  slab,  rise 
tine  to  peak  pressure,  and  decay  characteristics  of  the  pressure  pulse.  Load 
parameters  for  each  of  the  tests  are  shown  in  Table  3. 
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TABLE  3.  DYNAMIC  SHEAR  TEST  LOAD  PARAMETERS 


Element 

Group 

Charge 

density 

lb/in3 

DS1 

1 

1.37 

DS2 

1 

1.83 

DS3 

1 

0.91 

DS4 

?. 

1.37 

DS5 

1 

1.83 

DS2-1 

2 

2.29 

DS2-2 

2 

1.83 

DS2-3 

2 

1.14 

DS2-4 

3 

2.29 

DS2-5 

3 

1.60 

DS2-6 

3 

1.14 

Average  peak 

overpressure.  Rise  time, 
lb /in2  ms 


4000 

0.05 

4700 

0.05 

3300 

0.07 

3500 

0.05 

5000 

0.10 

6000 

0.05 

6000 

0.10 

3200 

0.05 

6000 

0.05 

5500 

0.10 

4000 

0.05 
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AVAILABLE  EXPERIMENTAL  DATA  . 

Information  deemed  important  for  structural  response  analysis  prior  to 
the  testing  of  the  structures  included  the  following  quantities:  airblast 
pressure,  interface  pressures,  active  and  passive  strains  in  the  longitudinal 
reinforcement,  soil  stress,  accelerations  (tests  DS3,  DS4,  and  DS5  only),  and 
slab  deflections.  In  order  to  acquire  these  data,  the  tests  were  instrumented 
as  shown  in  Figure  5. 

The  blast  pressure  gages  measured  the  actual  pressure  generated  by  the 
foam  HEST  charge,  whereas  the  interface  pressure  gages  measured  the  pressure 
transferred  from  the  soil  layer  above  the  roof  to  the  roof  slab  (i.e.,  the 
actual  load  the  slab  "feels").  Active  strain  gages  on  the  steel  were  used  in 
all  tests  to  measure  compressive  and  tensile  stresses  in  the  roof  slab, 
whereas  a  passive  measure  of  ductility  was  obtained  by  scoring  the  rebar  at 
specific  intervals  along  the  length  prior  to  load  application.  Finally, 
deflection  profiles  of  the  underside  of  the  roof  slab  element  were  obtained 
through  the  use  of  high-speed  photography. 

A  complete  summary  of  the  digitized  data  recovered  from  the  eleven  tests 
can  be  found  in  (Ref.  2).  In  general,  the  recovered  data  for  airblast 
pressure,  interface  pressures,  accelerations,  and  deflections  were  good. 
Unfortunately,  some  of  the  data  were  "clipped"  at  very  early  time  because  of 
extreme  pressures,  were  subject  to  recording  noise  due  to  very  high 
frequencies,  or  were  improperly  digitized  because  of  calibration  errors.  Due 
to  these  circumstances,  some  of  the  data  (active  strain  for  example)  were 
found  to  have  only  limited  use  or  increased  uncertainty. 

Despite  these  limitations,  however,  much  of  the  data  were  discovered  to 
contain  very  useful  information  pertaining  to  the  structural  response  of  the 
reinforced  concrete  box.  In  the  next  section,  the  exoerimental  information 
used  for  this  study  will  be  discussed.  The  subjective  counterpart  to 
experimental  information  in  the  form  of  expert  ooinion  will  also  be  described 
in  Section  III. 
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Figure  5.  Test  instrumentation  layout ,  Test  Groups 
II  and  III  (Ref.  2) . 


14 


III.  ANALYSIS  OF. SUBJECTIVE  AND  OBJECTIVE  INFORMATION 


OVERVIEW  OF  DAMAGE  TO  CONCRETE  BOXES 

In  order  to  have  a  complete  understanding  of  the  buried  protective 
structure  problem,  a  description  of  the  danage  sustained  by  the  test  elements 
is  presented.  Following  this,  a  brief  outline  of  paraaeters  affecting  the 
dynamic  behavior  of  the  structure’s  response  is  given.  By  first  previewing 
this  information,  it  is  possible  to  proceed  to  the  more  important  topics  of 
subjective  data  and  objective  data  used  to  qualify  and  quantify  nodes  and 
levels  of  daeage. 

Photographs  of  the  eleven  tests  are  shown  in  Appendix  A.  Since  the 
original  purpose  of  the  tests  was  to  study  direct  shear  failure  under  very 
high  overpressures,  it  was  only  natural  that  the  Majority  (7  of  11)  of  the 
tests  were  subjected  to  complete  collapse  of  their  respective  roof  element, 
whereas  only  four  concrete  test  boxes  were  loaded  to  noncollapse  levels. 
Although  four  of  the  eleven  test  elements  did  not  completely  collapse,  it  has 
been  reasoned  that  all  boxes  failed  Initially  in  a  shear  or  shear-f lexure  mode 
(Ref.  4).  However,  upon  Inspection  of  the  posttest  photographs,  the  variation 
in  post-failure  danage  is  found  to  be  fairly  large. 

As  was  predicted  prior  to  the  tests,  several  of  the  specimens  failed  in  a 
direct  shear  mode.  Inspection  of  element  DS3,  for  example,  reveals  that  the 
roof  slab  was  completely  severed  from  the  walls  of  the  structure  along  nearly 
vertical  failure  planes.  Approximately  60X  of  the  principal  reinforcement  was 
broken  along  these  two  faces,  together  with  slight  necking  down  of  the 
reinforcing  bars.  The  reeainder  of  the  reinforcing  bare  were  pulled  out  of 
the  roof  supporting  wail  during  failure.  Inward  translations  of  the  tops  of 
the  walls  were  3  to  4  in  froe  the  vertical.  Considerable  cracking  was  noted 
along  the  base  of  the  structure,  as  well  as  at  the  intersection  of  the  wall 
and  the  base  slab.  Most  of  the  concrete  in  the  roof  slab  which  fell  was  found 
to  be  crushed  and  easily  crumbled  except  in  the  middle  third  of  the  slab. 

In  a  slightly  different  scenario,  inspection  of  test  element  DS2-1 
Indicates  that  all  reinforcing  bars  along  the  west  wall  (see  Figure  A-7  in 
Appendix  A)  were  completely  severed  along  a  near  vertical  surface,  as  in  DSS 
above.  Unlike  DS5  though,  the  east  failure  plane  was  not  this  same  type. 
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Failure  obviously  occurred  along  the  west  wall  First,  Followed  by  rotation  and 
reinf orceaent  pullout  along  the  east  wall  until  collapse.  This  asymmetric 
Failure  could  have  been  caused  by  a  nonsymmetr ical  load,  and/or  poor 
workmanship  during  construction. 

Test  elements  DS3  and  DS2-5  represent  cases  in  which  the  structure 
sustained  severe  to  very  extensive  damage  but  had  enough  residual  strength  to 
resist  complete  collapse.  Permanent  deFlections  oF  10  and  12  in  For  test 
elements  DS3  and  DS2-S,  respectively,  were  recorded.  In  both  cases,  large 
vertical  oFFset  (abrupt  edge  displacement)  and  sharp  rotation  oF  reinForcing 
steel  at  the  slab/wall  interFace,  indicate  an  initial  shearing  action.  While 
the  more  heavily  loaded  DS2-5  is  curved  throughout  its  span,  the  middle  third 
oF  test  element  0S3  remained  relatively  Flat.  Note  also  the  large  inward 
translations  oF  the  tops  oF  the  walls  in  element  0S2-5.  The  concrete  cover 
was  broken  up  over  the  entire  span,  and  virtually  all  the  concrete  cover  on 
the  underside  oF  slab  032-5  spalled  oFF,  exposing  the  principal  reinForcing 
steel.  Similar  damage  can  be  seen  in  the  DS3  element  but  not  nearly  to  the 
same  degree.  Perusal  oF  the  load  parameters  in  Table  3  shows  that  the  roof 
element  oF  test  DS3  was  loaded  to  approximately  2000  ib/in2  less  than  the  0S2- 
3  element.  Thus,  by  Inspection,  it  can  be  seen  that  these  two  tests  responded 
similarly)  but  the  late-time  impulse  in  the  First  case  was  not  enough  to  Force 
it  to  the  same  level  oF  damage  as  In  the  second  case. 

Observations  oF  test  elements  DS2-6  and  DS2-3  show  instances  where  the 
level  o4  damage  attained  was  only  slight  to  moderate  when  comoared  to  the 
tests  discussed  previously.  Permanent  deFlection  is  3  to  4  in  For  each 
specimen.  Note  again,  the  response  oF  both  of  the  structures  appears  to  be 
predominantly  shear,  as  can  be  seen  by  the  vertical  offset  of  the  roof  slab  at 
the  slab/wall  interface.  The  slabs  remained  relatively  Flat  with  some 
cracking  on  the  underside  of  their  central  portions,  host  of  the  cracking  and 
crushing  of  the  concrete  appears  in  the  region  within  a  distance  equal  to  the 
slab  thickness  From  thm  wall. 

PARAMETERS  AFFECTING  DYNAMIC  RESPONSE 

The  response  of  slabs  u-der  dynamically  applied  loads  is  a  Function  of 
the  load  magnitude  and  distribution,  stab  character! sties,  and  boundary 
conditions.  Thus,  because  of  the  comolicated  interaction  of  variables 
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involved  in  slab  behavior  under  simulated  blast  overpressures,  there  is  no 
simple  way  to  predict  modes  of  structural  response.  Therefore,  before 
proceeding  with- an  investigation  of  soft  and  hard  data  that  are  available  for 
describing  aodes  of  response,  a  brief  overview  of  parameters  affecting 
structural  behavior  is  warranted. 

In  studying  the  response  of  these  structures,  it  is  important  to  have  an 
idea  of  whether  shear  or  flexure  typically  controls  response,  and  how  and  when 
shear  and  flexure  each  attain  their  failure  levels.  The  tern  failure  is 
defined  here  as  the  point  at  which  the  concrete  element  reaches  its  ultimate 
capacity  (either  in  shear  or  in  flexure).  It  is  also  important  to  note  at 
this  time  that  failure  in  a  given  mode  does  not  imply  collapse  of  the 
structure. 

It  is  known  that,  under  rapid  rates  of  loading  and  strain,  the  apparent 
material  strengths  of  both  steel  and  concrete  are  significantly  increased,  as 
shown  in  Figure  6.  Because  of  this,  the  determination  of  the  dynamic  ultimate 


Figure  6.  Dynamic  strain  rat©  effects  on  strength  of 
concrete,  steel,  and  aluminum  (Ref.  51. 
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capacity  of  a  member  is  a  difficult  task.  Therefore,  in  the  case  of  an 
impulsive  blast  load,  it  is  especially  important  to  understand  the  effects 
induced  by  changes  in  strain  rate  in  order  to  predict,  design,  and  protect 
against  failure. 

Ross  (Ref.  4)  used  a  viscoelastic  Timoshenko  beam  to  study  the  effects  of 
load  rate  on  incipient  elastic  failure  at  slab  supports  (assuming  that  one-way 
slab  response  ir  similar  to  that  of  a  bear).  More  specifically,  he  studied 
the  combination  of  conditions  which  were  necessary  to  produce  either  direct 
shear  or  flexure  failure  and  their  relative  likelihood  of  occurrence.  His 
findings  included  the  fact  that  shear  force  tends  to  be  amplified  more  than 
bending  moment  in  early  time  (less  than  1  ms)  because  of  strain  rate  effects. 
His  findings  indicate  that  time  to  failure  in  shear  is  significantly  decreased 
when  strain  rate  effects  are  included,  thus  predicting  shear  dominance  aver 
flexure  in  early  time. 

Ross  also  studied  the  early-time  response  effects  of  load  parameters, 
beam  end  restraint,  and  length-to-depth  ratios  under  the  same  types  of 
impulsive  load  conditions.  Excepting  load  duration,  all  parameters  were  f^und 
to  have  some  effect  on  the  response  nf  the  structure.  Direct  shear  failures 
precede  flexural  failures  at  early  time  for  certain  combinations  of 
parameters.  Failure  curves  were  produced  which  related  two  of  the  most 
significant  load  parameters,  peak  pressure  and  rise  time,  to  beam  parameters. 
Typical  curves  are  shown  in  Figure  7  for  slabs  with  parameters  in  group  3  (as 
given  in  Tables  1  and  2).  Note  that  points  falling  abjve  the  failure  curves 
correspond  to  direct  shear  failures  and  points  below  tho  curve  correspond  to 
flexure  failures.  By  comparison  with  the  actual  range  of  conditions 
experienced  in  the  test  series,  it  was  observed  that  peak  pressure  and  rise 
time  have  the  largest  effects  on  failure  mode. 

Another  dynamic  phenomenon  related  to  the  buried  boxes  is  pointed  out  at 
this  time.  Actual  testing  of  che  reinforced  concrete  buried  boxes  showed  that 
previous  analytical  methods  had  severely  underpredicted  the  actual  strength  of 
these  structures.  After  studying  the  pressure  gages  located  at  the  soil/roof 
interface,  it  was  discovered  that  significant  late  time  (greater  than  1  ms) 
transferral  of  pressure  from  the  flexible  central  portion  of  the  roof  to  the 
more  rigid  supports  was  taking  place.  This  is  due  to  a  phenomenon  known  as 
soil  arching.  Basically,  soil  arching  is  the  ability  of  a  soil  to  transfer 
loads  through  a  system  of  shear  stresses  from  one  location  to  another  in 
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(c)  Influence  of  L/d  on  failure  curves 
(1  KPa  -  H5  lb/in2). 

Figure  7.  Typical  failure  curves  with  Group  3-cype 
properties  and  fixed  ends  (Ref.  $). 
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response,  to  a  relative  displacement  between  the  locations.  Although  the 
effects  of  soil  arching  are  usually  neglected  when  the  depth  of  burial  is  less 
than  the  clear  span  of  the  roof,  its  effects  on  late-time  interface  pressure 
distribution  over  the  slab  length  can  be  rather  significant.  Although  the 
arching  phenomenon  does  not  influence  early  tine  incipient  failure  conditions, 
it  can  be  a  significant  feature  to  consider  in  late  tine,  if  there  is  a 
sufficient  interface  iapulse  to  cause  collapse. 

SUBJECTIVE  INFORMATION 

Obtaining  Exoert  Opinion 

As  mentioned  in  Section  I,  the  assessment  of  daaage  to  a  structure  aay  be 
a  combination  of  both  subjective  and  objective  inforaation,  as  seen  in  Figure 
B.  In  this  report,  the  tera  subjective  information  is  aeant  to  describe 
inforaation  that  resides  in  the  ainds  of  experts  in  a  particular  field  who 
have,  over  time,  accumulated  in-depth  knowledge  in  their  subject-field. 
Judgment  of  this  type  is  used  often  in  engineering  evaluations,  especially  to 
compensate  for  sparse  data  and  for  extrapolation  outside  a  given  data  range. 

For  this  study,  the  engineering  judgment  was  obtained  by  polling  experts 
in  the  area  of  structural  daaage.  k  questionnaire  concerning  the  daaage  to 
the  roof  eleaent  of  the  11  test  structures  described  previously  was  sailed  to 
60  structural  engineers  across  the  United  States  (33  responded).  The 
questionnaire  contained  posttest  photographs  (see  Appendix  A)  of  all  the  test 
elements. 

The  purposes  of  the  questionnaire  were  threefold.  First,  it  was  hoped 
that  inforaation  obtained  from  the  questionnaire  could  be  used  to  help 
determine  how  experts  arrive  at  their  evaluations  of  conceptually  difficult 
probleas.  The  second  purpose  was  to  determine  how  or  why  their  answers  were 
the  same  cr  different.  The  third  purpose  was  to  determine  if  it  is  possible 
to  aggregate  these  assessments  into  a  single  data  base. 

The  format  of  the  questionnaire  was  quite  simple.  For  each  picture,  the 
expert  was  asked  two  questions  dealing  with  the  daaage  shown-  in  the  given 
picture.  First,  "What  is  (are)  the  aode(s)  of  response  of  the  roof  element?'1; 
and  second,  "What  is  the  degree  of  daaage  to  the  roof  in  that  mode?"  Along 
with  their  responses  to  the  two  questions,  the  experts  were  asked  to  provide 
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their  reasonings  to  substantiate  their  answers.  The  first  question  was 
designed  to  be  qualitative  in  that  linguistic  responses  were  expected. 
Acceptable  answers  night  include  "direct  shear,"  or  "flexure  followed  by 
tension  membrane,"  etc.  The  second  question  was  intended  to  be  nore 
quantitative  in  content  (example:  slight  danage  =  1/10,  or  total  damage  = 
10/10),  although  linguistic  answers  with  appropriate  explanations  were  also 
acceptable. 

The  questions  asked  of  these  experts  were  relatively  vague  in  content  and 
were  designed  as  such  in  order  to  avoid  any  built-in  bias  on  the  part  of  the 
author.  Uncertain  warding  of  the  questions  allows  the  expert  to  interpret  the 
meaning  without  external  constraints  forcing  his  judgment  and  allows  the  most 
information  to  be  derived  from  a  single  question.  On  the  other  hand,  inherent 
vagueness  within  a  question  also  increases  the  possibility  of  confusion  and 
noncomprehension. 

Analysis  of  Expert  Opinion 

One  of  the  main  problems  encountered  with  the  analysis  of  the 
questionnaire  was  the  difference  of  opinions  expressed  in  their  reasoning  by 
the  experts.  Varying  backgrounds  among  the  experts  is  on*  possible 
explanation,  because  some  of  thee  were  familiar  with  the  experimental  test 
series  and  some  wsre  not.  Thus,  some  of  those  interviewed  had  unfair  insight 
into  the  questions  asked.  A  second  problem  encountered  when  analyzing  the 
responses  was  the  interpretations  of  ambiguous  technical  and  linguistic  terms, 
such  as  "severely  damaged,"  or  "deformation  near  wall  is  small."  Another 
problem  was  associated  with  a  lack  of  information.  Some  of  the  experts  were 
hesitant  or  confused  in  their  responses  simply  because  specific  information  on 
structural  parameters,  functionality  issues,  or  repairability  problems  were 
not  provided  to  them. 

Besides  the  large  amount  of  raw  data  obtained,  other  positive  aspects 
were  ascertained  from  the  questionnaire.  First  of  all,  it  was  discovered 
that,  in  general,  experts  were  comfortable  in  qiving  their  opinions  even 
though  the  pictures  submitted  to  them  were  poor  quality  and  the  information 
was  incomplete.  Secondly,  many  of  the  experts  were  able  to  predict  the  actual 
chronology  of  events  leading  to  deformation  and/or  failure  using  only  visual 
data. 


22 


Having  considered  the  initial  impressions  of  the  experts’  responses,  the 
most  important  questions  that  might  be  asked  next  are,  "What  do  we  do  with  all 
this  raw  information?",  and  "How  do  we  decompose  it  into  knowledge  tha-t  can  be 
implemented  into  a  computer  code?"  As  pointed  out  above,  the  greatest 
difficulty  arises  when  we  examine  the  data  and  discover  the  variety  of 
linguistic  terms  put  forward  by  the  experts  and  the  wide  variety  of  images  or 
situations  these  terms  describe.  Initially,  some  confusion  arises  in  trying 
to  decide  whether  different  terms  represent  different  explanations,  or  if  they 
are  using  different  terminology  to  describe  the  same  concept.  For  instance, 
some  of  the  terms  used  to  describe  a  shearing  action  include:  shear,  punching 
shear,  pure  shear,  direct  shear,  vertical  shear,  edge  shear,  and  sliding 
shear. 

After  careful  analysis  of  the  data,  however,  it  became  obvious  that  these 
problems  are  somewhat  superficial  and  that  there  was  a  tendency  (whether 
conscious  or  unconscious)  on  the  part  of  the  experts  to  break  their  reasoning 
up  into  smaller  components  of  the  problem.  For  instance,  on  the  topic  of 
damage  level,  the  experts  tended  to  place  their  reasonings  into  one  of  three 
categories;  structural  integrity,  functionality  or  use  of  the  structure,  and 
repairabil ity  of  the  structure.  Typical  responses  given  for  "slight," 
"moderate,"  and  "severe"  damage  levels  appear  in  Tables  4  through  6.  As  can 
be  seen  from  the  tables,  most  experts  grouped  their  reasoning  into  the 
category  of  structural  integrity.  This  is  understandable,  as  they  were  given 
no  information  other  than  the  photographs  and  thus  relied  heavily  on  the 
structural  portion  of  their  intuitive  knowledge. 

On  the  subject  of  modes  of  failure,  further  investigation  revealed  that, 
within  the  category  of  structural  integrity,  the  experts  tended  to  specify 
their  reasonings  depending  on  the  mode  of  response  they  believed  was  involved. 
Virtually  all  descriptions  were  closely  related  to  one  of  the  four  major  mode 
groups:  shear,  diagonal  tension,  flexure,  or  tension  membrane. 

From  this  point,  it  was  observed  that  the  data  could  be  further 
subdivided  into  specific  structural  attributes  or  parameters  dealing  with 
three  unique  locations  along  the  span  of  the  roof,  i.e.,  the  middle  portion  of 
the  slab,  the  slab  near  the  walls,  and  the  wall  supports.  Tables  7  and  0  show 
some  of  these  attributes  along  with  typical  descriptions  volunteered  by  the 
experts  for  shear  and  flexural  nodes.  Similar  tables  for  diagonal  tension  and 
tension  membrane  can  also  be  created. 
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TABLE  4.  EXPERT  REASONING  GIVEN  FOR  SLIGHT  DAMAGE 


minimal  permanent  deflection 

damage  localized  at  supports 

no  rebar  problems  and  minimal  dispacement 

small  support  rotations 

still  in  good  shape 

minimal  tensile  cracking  at  centerline  and  support 
considerable  capacity  remains 
no  plastic  hinge  formation 
small  crack  lines 

most  structural  resistance  remains 


-  still  usable 

-  system  remains  functional 

-  everything  inside  should  have  survived 

-  can  still  be  used  for  the  purpose  for  which  it  was 
designed 

-  reusable 


-  small  permanent  deflection  can  be  repaired _ 


TABLE  5.  EXPERT  REASONING  GIVEN  FOR  MODERATE  DAMAGE 


STRUCTURAL_INTEG.Rm 

-  end  roof  still  in  place  even  though  significant 
spalling  has  occurred 

-  plastic  hinge  not  fully  developed 

-  onset  of  membrane  action 

•  rebar  cage  still  intact 

-  lots  of  cracking  and  permanent  set 

•  some  tensile  pullout  and  shear  punch  on  right  side 

FUNCTIONALITY 

-  some  unobstructed  clear  space  is  provided 

-  slight  rotation  of  walls 

-  contents  probably  OK 

-  most  equipment  would  survive 

**  small  debris  would  fall  on  contents 

-  structure  may  still  be  useful  for  something 

REPAIRABILITY 

-  structure  serviceable 

-  minimum  repair  needed 

-  probably  repairable 
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TABLE  6.  EXPERT  REASONING  GIVEN  FOR  SEVERE  DAMAGE 


STRUCTURAL  INTEGRITY 

-  very  large  permanent  deformation 

-  severe  spalling  at  joints 

-  most  of  concrete  has  spalled  off  roof 

-  yielded  reinforcing,  spalled  concrete 

-  on  verge  of  collapse 

-  large  support  rotations 

-  no  longer  able  to  withstand  blast  pressures 

-  shear  hinge  at  support  is  fully  developed 

-  end  supports  nearly  broken 

EUKS.TIQNA.LITX 

-  only  short-term  use  is  possible 

-  reasonable  chance  for  survival 

-  spalling  would  have  harmed  contents 

-  contents  severely  shaken 

-  significant  debris  and  pressure  ingress 

-  associated  shock  and  vibration  would  have  severely 
damaged  contents 

REPAIRABI LITY 

-  needs  a  lot  of  work 

-  roof  needs  replacing 

-  could  possibly  be  repaired  for  temporary  use 
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TABLE  7.  ( EXPERT  REASONING  GIVEN  FOR  SHEAR  DAMAGE 


Location 

Attribute  or 
Parameter 

Expert  Description 

Main  slab 

Curvature 

-  relatively  flat 

(middle  1/3 
to  1/2  span) 

-  rebar  cage  appears  to  be 
flat 

-  almost  uniform  displacement 

Crack  Pattern 

-  not  many  cracks  visible 

-  no  yield  lines  @  centerline 

-  no  crushing  in  the  compres- 
!  sion  zone  @  top  centerline 

Slab  Near 

Displacement 

-  clean  vertical  break 

Wall 

-  local  vertical  deformation 

-  mostly  at  supports 

-  relative  displacement  at 
edges 

-  sharp  gradient 

Anchorage 

-  sharp  bending  of  bars  at 
support 

-  broken  bars 

-  apparently  yielded 
reinforcement 

-  rebar  severed  at  wall 

-  rebar  violently  ripped  out 

Crack  Pattern 

-  diagonal  cracking 

-  can  see  diagonal  struts 

spalling/ 

Crushing 

-  considerable  concrete 
crushing 

Walls 

Rotation 

-  minimal  inward  rotation 

Spalling/ 

Crushing 

-  wall  supports  remain  intact 
with  little  concrete 
crushing 
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TABLE  8 


EXPERT  REASONING  GIVEN  FOR  FLEXURAL  DAMAGE 


Location 

Attribute  or 
Parameter 

Expert  Description 

Main  slab 

Curvature 

-  smooth  top  bar  curvature 

(middle  1/3 

-  moderately  curved  roof 

to  1/2  span) 

-  interior  flexural  hinge 
formation 

Crack  Pattern 

-  cracks  indicate  3 -hinge 
mechanism 

-  most  cracks  are  on  the 
bottom  at  the  centerline 

-  crack  and  deformation 
pattern 

Spalling/ 

Crushing 

-  crushed  concrete  at  top 
centerline 

-  spalling  at  centerline 

Displacement 

-  large  displacement  at 
centerline 

Slab  Near 
wall 

Curvature 

-  large  rotation  near  wall 

Anchorage 

-  bars  pulled  out  without 
evidence  of  yielding 

Displacement 

-  lack  of  vertical  offset 
at  wall 

Crack  Pattern 

-  cracks  indicate  flexural 
hinge  formation 

-  tensile  cracks  at  top 
edge 

-  crack  lines  at  top  face 

Walls 

Rotation 

-  minimal  inward  rotation  of 
walls 

Spalling/ 

Crushing 

-  concrete  crushed  in 
compression  zone  of  wall 
face 
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Interestingly,  if  we  take  a  close  look  at  the  attributes  described 
previously,  we  notice  they  are  very  similar  to  those  used  to  describe  the 
levels  of  damage.  For  instance,  "large  deflection  at  wall"  could  be  used  to 
describe  a  shear  mode,  but  it  could  also  be  used  as  an  indication  of  severe 
damage.  Thus,  by  reducing  or  subdividing  the  expert  information  into  lower 
levels  of  knowledge  as  seen  in  Figure  9,  attributes  or  parameters  can  be 
defined  that  describe  both  modes  of  failure  and  levels  of  damage. 

Note  too  that,  by  carefully  choosing  the  wording  for  these  attributes, 
the  meaning  can  be  varied  by  changing  only  the  adjective  used  to  describe  it. 
In  the  example  given,  for  instance,  the  attribute  is  "deflection  at  the  wall." 
By  simply  changing  the  associated  adjective  from  "large"  to  "none  (no),"  we 
change  the  meaning  and  are  thus  able  to  use  the  same  attribute  over  a  wide 
range  of  conditions.  However,  an  obvious  problem  arises  because  the  meaning 
of  such  terms  as  "large,"  “slight,"  “none,"  etc.,  oust  now  be  defined. 

Once  we  have  solved  this  problem  of  terminology,  we  have  a  unique  set  of 
linguistic  attributes  that  we  can  use  to  describe  modes  of  response  and  levels 
of  damage.  In  all,  ten  of  these  structural  attributes  or  parameters  were 
identified  within  the  experts*  responses.  These  ten  attributes  are  listed  in 
Table  9.  The  use  of  these  attributes  in  the  knowledge  combination  process  is 
illustrated  in  Section  IV.  More  specifically,  expert  opinion  relating  modes 
of  damage  to  the  attributes  will  be  combined  with  (user  supplied)  observations 
of  attribute  levels  to  produce  output  identifying  modes  and  levels  of  damage. 

m 

OBJECTIVE  INFORMATION 

As  stated  in  the  introduction,  one  of  the  goals  of  this  study  was  to  be 
able  to  combine  objective  data  with  the  subjective  information  outlined  above. 
Before  this  can  be  done,  it  must  be  determined  which  data  are  the  most  useful, 
and  what  is  the  most  efficient  manner  in  which  to  implement  the  data  into  the 
knowledge  base.  Previously,  when  we  referred  to  objective  information,  it  was 
stated  that  we  were  talking  about  "hard"  data  such  as  digitized  waveforms, 
etc.  Although  these  data  are  crisp  in  the  sense  that  the  data  are 
quantifiable,  this  information  is  not  without  inherent  sub jecti veness  when  it 
is  applied  to  the  difficult  task  of  interpreting  modes  of  failure. 
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Hades/ Degree  of  Damage 
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Pigure  9.  Reducing  expert  opinion  to  lower  levels  of  knowledge. 


TABLE  9 

Attribute 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


STRUCTURAL  ATTRIBUTES  DEVELOPED  FROM  EXPERT 
RESPONSES 


Description 


Spalling  and  loss  of  concrete/ rebar  interaction 
on  the  bottom  side  of  the  main  slab  (central 
1/3  to  1/2  span) 

Crushing  of  concrete  on  the  top  side  of  the 
main  slab  (central  1/3  to  1/2  span) 

Vertical  offset  (step  function  shape  with  slab 
top  displacement  relative  to  edge  support)  in 
the  slab/vall  connection  region 

Rotation  of  the  roof  slab  in  the  slab/wall 
connection  region 

Curvature  of  the  main  slab  (central  1/3  to  1/2 
span)  with  corresponding  differential 
deflection  (centerline  minus  near  wall)  to  slab 
length  ratio 

Flexural  related  cracking  (including  yielding, 
hinge  formation,...  of  the  centerline  main  slab 
and  tensile  cracks  at  the  top  side  of  slab/wall 
connection 

Inclined  cracking  at  roughly  t/2  from  the  wall 
face  (at  approx.  45  degrees  to  the  horizontal) 

Spalling  and  loss  of  concrete/rebar  interaction 
in  the  slab/wall  connection  region 

Inward  rotation  of  the  tops  of  the  walls 

Crushing  of  concrete  and  loss  of  concrete  cover 
or  rebar  anchorage  in  the  wall  supports 
themselves,  or  in  the  slab/wall  joint 
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Damage  Modes 


The  amount  of  information  that  can  be  obtained  from  digitized  waveforms 
is  enormous.  From  interface  pressure  data  alone,  the  number  of  features  that 
can  be  derived  to  aid  the  mode  identification  process  is  almost  endless. 
Features  such  as  rise  time  and  peak  pressure,  frequency  in  terms  of  the  power 
spectral  density,  ratios  of  various  impulse  values,  and  decay  slope  of  the 
pressure  record  after  peak  pressure  are  but  a  few  of  the  possibilities  to  be 
investigated.  Whether  or  not  each  feature  contains  information  regarding 
modes  of  failure  is  a  separate  question. 

For  this  study,  experimental  test  equipment  provided  numerical  data  for 
analysis  of  roof  slab  rtsponse  modes  in  the  form  of  digitized  interface 
pressure  waveforms,  deflection  profiles  from  high-speed  photography,  passive 
and  active  steel  strains,  and  velocity/acceleration  versus  time  records.  The 
significance  of  each  of  these  objective  parameters  is  discussed  next. 

For  each  of  the  eleven  tests,  interface  pressures  were  available  at  three 
locations  along  the  roof  span  as  was  shown  in  Figure  5.  Ross  and  Krawinkler 
(Ref.  7)  pointed  out  that  the  pressure  decay  after  initial  peak  was  a  good 
indicator  of  early  time  response  (less  than  1  ms)  of  the  roof  slab.  They 
noted  that  pair-wise  comparison  of  the  three  pressure  records  wouid  indicate 
whether  shear  or  flexure  response  dominates  at  early  time. 

In  the  event  of  a  flexure  controlling  failure,  pressure  plots  are  similar 
for  locations  1  and  3  and  different  than  location  2,  as  shown  in  Figure  10. 
Rapid  decay  of  the  relatively  flexible  centerline  of  the  slab  indicates 
movement  downward,  while  readings  near  the  support  and  over  the  wall  remain 
higher  for  a  longer  period  of  time.  These  sustained  high  pressures  indicate  a 
rigid  boundary  condition  with  these  two  locations  (1  and  3)  "seeing"  little 
downward  movement. 

In  the  event  of  a  shear  controlling  failure,  pressure  plots  at  locations 
2  and  3  become  similar  in  shape  and  decay  much  more  rapidly  than  those  at 
location  1,  as  shown  in  Figure  11.  Here  again,  the  pressure  at  location  1 
stays  higher,  longer,  indicating  a  rigid  boundary  condition.  Rapid  decay  of 
pressure  at  location  2  and  location  3  indicates  downward  movemant  of  the 
entire  slab  as  a  rigid  body.  The  only  way  for  this  to. happen  would  be  for  an 
initial  slip  to  occur  at  the  support. 
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2000 


Figure  10.  interface  pressure  plot  of  typical  flexural  failure  response 


The 


As  mentioned  previously,  this  information  is  not  purely  “objective, 
graphs  in  Figures  10  and  11  are  for  two  tests  conducted  prior  to  the  11 
dynamic  shear  tests  (see  Ref.  3).  The  flexure  and  shear  conditions  are 
evident  in  these  two  cases.  But  in  the  event  that  the  similarity  of  the 
pressure  at  specified  locations  is  not  this  clear-cut,  how  does  one  go  about 
deciding  just  how  "similar"  the  graphs  are  at  each  location? 

Active  steel  strains  (shearing  and  normal)  are  useful  in  providing 
information  regarding  whether  a  shearing  action  or  flexural  action  has 
occurred.  As  pointed  out  in  Section  II,  however,  because  of  high  recording 
noise,  strain  measurements  at  early  time  were  indiscernible.  However,  late 
time  response  of  the  slab  element  in  the  form  of  flexure  or  membrane  action 
can  be  predicted.  In  the  case  of  a  fixed-fixed  beam,  flexure  is  denoted  by 
tension  in  the  top  reinforcement  and  compression  in  the  bottom  reinforcement 
at  the  support,  with  the  reverse  being  true  at  the  centerline.  A  complete 
membrane  response  is  indicated  by  tension  in  all  reinforcement.  If  initial 
shear  controls  because  of  dowel  action  at  the  support,  i^  is  possible  for  both 
top  and  bottom  reinf orcement  to  be  in  tension,  while  the  central  portion  of 
the  slab  responds  as  a  flexural  element  with  compression  in  the  top  and 
tension  in  the  bottom.  Passive  strains  obtained  from  scored  rebar  have  little 
use  except  for  measuring  overall  ductility  associated  with  total  collapse. 

In  general,  velocity/acceleration  graphs  can  provide  insight  to  trends  in 
the  response  when  combined  with  shear  strain  data.  For  instance,  modes  of 
damage  are  recognized  in  the  data  by  certain  features  in  the  midspan  velocity- 
time  histories.  A  dip  in  the  velocity  following  the  initial  peak  indicates 
flexure  is  dominant,  while  a  monotonic  increase  to  peak  velocity  indicates 
shear  is  dominant.  The  reason  for  these  trends  is  based  upon  rigid  body 
motion  concepts  similar  to  those  discussed  previously  for  the  interface 
pressures.  Because  of  the  fact  that  only  three  of  the  eleven  test  elements 
(DS3-DS3)  were  equipped  with  acceleration  gages,  more  velocity/acceleration 
data  are  needed  for  analysis  before  mode  predictions  can  be  made  with  any 
confidence. 

Typical  deflection  profiles  obtained  from  high-speed  photography  for 
tests  DS3  and  DS2-3  appear  in  Figures  12  and  13,  respectively.  Rotation  at 
the  support,  curvature  of  the  central  portion  of  the  span,  and  differential 
deflection  ratios  from  various  sections  along  the  roof  span  can  provide 
measures  of  mode  response,  large  vertical  offset  and  sharp  rebar  rotation  at 
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Figure  12.  Deflection  profiles  of  Test  DS3  (Ref.  2). 
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Figure  13.  Deflection  profiles  of  Test  DS2-5  (Ref.  2). 
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the  Mall  support,  combined  with  little  or  no  curvature  (i.e.,  -flat)  of  the 
eain  slab,  indicates  an  obvious  shear  response.  If  sufficient  residual 
strength  is  available  after  initial  shear,  the  slab  responds  in  a  flexural 
node  similar  to  that  of  a  sieple  bean.  And  finally,  continued  loading  nay 
force  the  entire  slab  into  a  membrane  response  where  all  fibers  are  in  tension 
and  a  large  curvature  prevails  at  the  centerline.  To  varying  degrees,  each  of 
these  responses  can  be  seen  in  Figures  12  and  13. 

Experimental  data  from  the  eleven  tests  provide  a  wealth  of  information 
for  determining  structural  modes  of  response.  Obviously,  not  all  data  are  of 
equal  value,  nor  is  it  within  the  scope  of  this  project  to  examine  and 
incorporate  all  this  information  into  the  computer  code.  Rather,  the 
objective  is  to  demonstrate  the  possibilities  that  are  available  for 
analyzing,  combining,  and  incorporating  different  types  of  data  into  a  rule- 
based  damage  assessment  system.  Therefore,  from  the  information  outlined 
above,  only  the  interface  pressure  data  and  deflection  profile  information 
will  be  used  in  the  code. 

The  reasons  far  choosing  these  two  features,  interface  pressures  in 
particular,  were  summarized  well  in  Reference  8.  First,  these  two  features 
were  the  most  complete  sets  of  data  from  the  eleven  tests.  Secondly,  they  are 
assumed  to  have  recorded  phenomena  significantly  above  the  noise  levels  of  the 
gages.  Thirdly,  these  two  features  were  used  in  previous  analytical  efforts 
in  attempts  to  characterize  failure  conditions.  Finally,  these  features  are 
an  actual  measure  of  the  overall  behavior  of  the  bean.  The  methodology  used 
for  implementation  of  these  data  into  the  code  is  detailed  in  Section  IV. 

Damage  levels 

A  significant  effort  has  bean  expended  over  the  last  several  years  to 
analyze  and  understand  the  structural  modes  of  response  controlling  the 
behavior  of  reinforced  concrete  boxes  to  impulsive  blast  loads 
(Refs.  3,5,8,9,10,11,12,13,14).  Although  nodes  of  response  are  important 
factors  in  this  study,  the  damage  level  imparted  to  the  structure  because  of 
these  various  modes  is  also  of  interest.  In  the  section  on  subjective 
information  it  was  found  that  careful  decomposition  of  expert  opinion  resulted 
in  linguistic  attributes  that  could  be  used  to  describe  both  modes  of  failure 
and  damago  levels.  Thus,  in  this  section,  we  seek  parameters  obtaineo  from 
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experimental  data  that  will  provide  a  measure  of  damage  imparted  to  the 
structure. 

Because  of  the  nature  of  the  loading  conditions,  the  structures  in  this 
study  can  be  expected  to  sustain  damage  within  a  broad  range  (possibly  even 
total  collapse)  as  seen  in  the  photographs  of  Appendix  A.  In  typical  civil 
engineering  structures,  it  is  generally  accepted  that  failure  occurs  in  a 
member  when  yielding  is  first  initiated.  Qn  the  other  hand,  a  protective 
structure  that  is  subjected  to  a  load  which  just  initiates  a  yield  condition 
in  the  roof  slab  might  be  considered  only  slightly  damaged  if  none  of  the 
contents  has  been  adversely  affected. 

From  the  above  discussion  we  can  conclude  that  the  measure  of  damage  in 
conventional  civil  engineering  structures  is  much  different  than  in  the  buriea 
protective  structure  scenario.  We  may  also  conclude  that  the  criterion  used 
in  the  measurement  of  damage  to  typical  structures  is  inadequate  for  our 
purposes.  Measures  of  damage  used  in  the  past  for  large  inelastic 
deformations  include  such  quantities  as  the  ductility  ratio,  rotation  ratio, 
or  curvature  ratio.  Damage  measures  such  as  these  simply  compare  the  maximum 
value  of  a  parameter  (deflection,  rotation,  or  curvature)  to  ti.e  first  yield 
value.  For  simple  models  and  static  loading  conditions,  these  quantities 
provide  adequate  damage  measures.  However,  experimental  comparisons  of  these 
damage  measures  in  structures  subjected  to  severe  load  conditions  show  they 
are  not  sufficient  to  produce  accurate  damage  assessment. 

Under  earthquake  loading  conditions,  for  example,  reinforces  concrete 
structures  are  damaged  by  a  combination  of  repeated  stress  reversals  (cyclic 
loading)  and  high  stress  excursions.  For  this  reason,  the  definition  of 
damage  in  terms  of  the  ductility  ratio  may  be  inadequate.  As  damage 
accumulates  in  a  structural  system,  its  strength  diminishes,  and  energy  is 
dissipated.  Energy  dissipation  then,  as  Banon  (Ref.  IS),  Park  (Ref.  14),  and 
Wang  (Ref.  16)  show,  can  be  an  effective  quantify  by  which  damage  is  measured. 
Calculation  of  the  dissipated  energy  produced  hy  earthquake  loads  is  a 
relatively  straightforward  procedure.  Single~degree*of-freedo«  (SQOF)  models 
can  be  used  to  calculate  the  incremental  energy  dissipated  from  hysteretic 
force-displacement  diagrams.  In  tin  problem  under  consideration,  however,  the 
loading  is  impulsive  rather  than  cyclic,  and  the  force-displacement 
relationship  is  much  more  complicated.  Thus,  application  and  implementation 
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of  similar  methods  to  the  buried  structure  problem  would  be  a  formidable  task 
in  itself. 

Rather  than  spend  unnecessary  amounts  of  time  developing  complicated 
damage  models,  the  authors  found  it  more  advantageous  to  investigate  the 
available  experimental  data  for  features  to  use  as  possible  damage  measures. 
Data  chosen  for  investigation  included  the  interface  pressure  data  and  the 
deflection  data  for  the  same  reasons  outlined  under  Damage  Modes.  Influenced 
by  the  discussion  above  on  energy,  initial  attempts  by  the  authors  to  find  a 
damage  level  indicator  were  concentrated  on  the  interface  pressure  data. 
Calculations  of  the  energy  per  unit  area  as  given  by  Equation  1  and  impulse 
per  unit  area  as  given  by  Equation  2  were  performed  for  each  test. 

P*d<$ 

Impulse  =  j"~  P  *dt 
fcG 

In  these  equations,  P  is  the  interface  pressure  at  a  specific  point  along  the 
roof  slab,  ^  is  the  deflection  of  the  roof  slab  at  a  specific  point  at  time 
tx  ,  and  62  it  the  deflection  of  the  roof  slab  at  the  same  point  at  time  t2 

Various  summations  and  ratios  of  these  values  were  tabulated  and  compared 
to  each  other  for  usefulness  as  damage  measures.  Little  correlation  was  found 
between  some  of  the  calculated  values  and  actual  damage  levels  of  the  tests. 
This  is  due  in  part  to  the  fact  that  pressure  data  are  available  for  only 
three  points  along  the  length  of  the  slab,  thus  making  it  difficult  to  obtain 
a  true  measure  of  the  entire  roof  response. 

Although  some  of  the  features  which  were  calculated  exhibited  fair 
within-group  trends  of  damage  (i.e.,  compared  well  with  linguistic  damage 
level  assessments),  none  of  the  features  studied  could  be  used  to  show  trends 
of  damage  for  all  groupings.  In  other  words,  none  of  the  features  was 
adequate  in  classifying  damage  for  all  It  tests. 

Deflection  profiles  were  examined  next  because  of  the  inherent 
information  contained  concerning  damage/def ormation.  The  calculation  for  one 
of  these  damage  features  is  shown  in  Table  10.  Note  that  the  aagnitudes  of 
the  normalized  values  compare  well  (within  the  group)  relative  to  the 
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(2) 


TABLE  10.  CALCULATION  OF  VALUES  FOR  A  TYPICAL  DAMAGE  LEVEL 
INDICATOR 


Damage  Indicator  (D.I.)  » 


(*10 

$Q»dt  + 

0 


rlO 

5wall*dt 

JO 


where  <S£  is  the  centerline  roof  slab  deflection  at  time  t, 
and  $wan  is  the  deflection  of  the  roof  slab  near  the  wall 
at  time  t. 


Test  No.  Group  Visual  D.I./  Normalized 

damage  in-ms  D.I. 

assessment 


FHla 

1 

Slight 

18.41 

0.22 

FH23 

1 

Total  (collapse) 

81.94 

1.00 

DS1 

1 

Total  (collapse) 

NAb 

DS2 

1 

Total  (collapse) 

74.69 

0.91 

DS3 

1 

Severe 

27.49 

0.34 

DS4 

1 

Total  (collapse) 

50.39 

0.61 

DS5 

1 

Total  (collapse) 

79.14 

0.97 

DS2-1 

2 

Total  (collapse) 

62.16 

0.76 

DS2-2 

2 

Total  (collapse) 

74.92 

0.91 

DS2-3 

2 

Moderate 

32.63 

0.40 

DS2-4 

3 

Total  (collapse) 

68.59 

0.84 

DS2-5 

3 

Very  Extensive 

48.01 

0.59 

DS2-6 

3 

Moderate 

NAb 

aData 

obtained  from 

Foam  HEST  Tests 

(see  Ref.  3) . 

^Deflection  data  not  recovered 

Note:  Values  were  normalized  by  dividing  by  the  largest 
D.I.  (Test  No.  FH2) . 
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linguistic  assessment  (i.e.,. total  damage  corresponds  to  values  approximately 
equal  to  1.0,  and  slight  or  moderate  damage  corresponds  to  values  much  less 
than  1.0).  As  with  the  energy  and  impulse  features,  few  of  the  deflection 
features  were  able  to  predict  damage  levels  for  all  11  tests. 

After  careful  study  of  all  the  features  calculated,  it  was  concluded  that 
the  phenomenon  of  within-group  agreement  is  caused  by  variation  in  the  major 
test  parameters  (f  ’,  P,  L/d).  Although  much  work  has  been  conducted  on  the 
interaction  of  these  parameters  with  respect  to  modes  of  response,  it  is  not 
known  with  confidence  how  these  parameters  affect  late-time  deformation 
response. 

Despite  all  the  problems  encountered  with  the  data,  a  rather  simple,  but 
effective  feature  based  on  deflection  profiles  was  chosen  for  inclusion  into 
the  computer  code  for  demonstration  purposes.  The  daaage  level  value  for  any 
given  test  is  computed  using  Equation  3. 


L 

f  5*dx 

Damage  Level  *  _0 _  (3) 

L  *t 

The  time  of  15  ms  was  chosen  because  most  deformation  has  taken  place  by  this 
time.  The  parameters  L  (slab  length)  and  t  (time)  are  inserted  into  the 
denominator  to  normalize  the  damage  indicator,  thus  simplifying  the  damage 
scale.  By  comparing  values  calculated  for  each  test  with  linguistic 
assessments,  a  scale  (shown  below)  was  created  for  classifying  damage  levels. 


None  Slight  Hoderate  Severe  Ex^en^iive  Collapse 

t . i . i . t . i . : . > 

0.0  0.2  0.4  0.7  1.0  1.3 


In  this  section,  all  the  information  needed  to  adequately  describe  the 
modes  of  response  and  levels  of  damage  present  in  the  buried  concrete  box 
structure  have  been  described.  The  only  remaining  task,  therefore,  is  to 
integrate  this  knowledge  into  an  efficient,  user-friendly  code.  In  the 
following  section,  the  methodology  used  to  accomplish  this  task  is  discussed 
in  detail.  This  discussion  includes  both  the  theory  behind  the  techniques  and 
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the  computer  algorithms  used. to  implement  them.  And,  in  order  to  help  the 
reader  understand  the  program  which  has  been  developed,  an  example  session 
using  the  DAPS  code  is  presented  in  Section  V. 
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IV,  THE  STRUCTURE  OF  DAPS 


Before  the  actual  organization  of  the  DAPS  code  is  discussed,  it  is 
necessary  to  present  the  underlying  theory  on  which  the  code  is  based.  The 
first  two  subsections  of  this  section  review  some  important  concepts 
concerning  expert  systems  and  fuzzy  set  theory.  This  review  should  help  the 
reader  understand  the  basis  on  which  the  code  is  developed. 

A  RULE-BASED  EXPERT  SYSTEM  APPROACH 

Expert  systems  are  an  outgrowth  of  a  subdiscipline  of  Artificial 
Intelligence  (AI)  research  begun  in  the  mid-1960s.  Early  researchers  in  the 
field  of  AI  were  interested  in  building  machines  capable  of  pattern 
recognition,  sensory  perception,  learning,  and  comprehending  the  semantics  of 
human  thought.  Although  the  major  use  of  computer  technology  over  the  years 
has  been  in  arithmetic  operations,  AI  focuses  on  the  use  of  the  computer  as  an 
idea  processor  rather  than  a  number  cruncher. 

Expert  systems  are  but  a  small  fragment  of  the  AI  field.  The  purpose  of 
an  expert  system  is  to  play  the  role  of  a  consultant  or  an  expert  operating  in 
a  very  restricted  domain  of  knowledge  who  gives  adyice  to  someone  with  a  task 
or  a  problem  within  that  domain.  In  other  words,  expert  systems  try  to 
capture  a  little  of  the  complex  process  of  human  reasoning.  In  its  simplest 
fora,  an  expert  system  is  a  collection  of  knowledge  in  the  fora  of  “rules-of- 
thumb,"  as  practiced  by  the  domain  expert,  together  with  special  techniques 
for  applying  the  proper  rule  at  the  proper  time. 

A  question  that  might  arise  at  this  point  is  “How  does  one  recognize  a 
good  candidate  problem  for  expert  system  implementation?H  A  response  to  this 
question  is  shown  below  as  a  list  of  identifiable  expert  system 
characteristics. 

•  A  relevant  body  of  knowledge  exists  and  is  available 

•  The  skill  involved  is  one  which  could  be  taught  to  a  new  employee 

•  The  knowledge  can  be  expressed  in  bite-sized  pieces  that  make  sense 
standing  alone 

•  Solving  the  problem  without  a  computer  takes  an  expert  no  less  than 
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a  few  Minutes,  and  no  nore  than  a  few  hours 

•  The  benefit  that  will  cote  froa  developing  the  systea  is  sufficient 
to  justify  the  cost  involved 

•  The  problea  contains  subjective  paraaeters  which  are  inherently 
difficult  to  quantify 

As  we  see  in  Figure  14,  the  two  aajor  parts  of  an  expert  systea  include  a 
knowledge  base  and  an  inference  aechanisa.  The  knowledge  base  is  the  storage 
facility  for  heuristic  rules;  it  atteapts  to  eaulate  the  coaplex  processes  of 
deductive  and  inductive  huaan  reasonirg.  The  heuristic  knowledge  in  the 
knowledge  base  has  been  collected  froa  one  or  several  huaan  experts  and  is 
coaposed  of  condition-action  stateaents  in  the  fora  of  IF-THEN  type  rules. 

The  inference  aechanisa  is  the  portion  of  the  code  separate  froa  the 
knowledge  base  which  uses  the  rules  together  with  data  froa  the  user  to 
“reason"  through  a  problea.  During  each  cycle,  the  conditions  of  each  rule 


Answer 

Figure  14.  Typical  expert  system  configuration. 
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are  Hatched  against  the  current  state  of  facts.  When  conditions  stated  in  the 
IF  portion  are  matched,  actions  are  executed  in  the  THEN  po'tion.  The  THEN 
statements  (consequents)  then  bee one  part  of  the  facts  base.  These  actions 
affect  the  current  state  of  facts,  Baking  it  possible  for  new  rule  conditions 
to  be  initiated.  The  inference  process  in  most  current  expert  systems  is 
keyed  to  a  search-and-aatch  strategy  of  this  type.  A  rule  Mil!  not  be 
triggered  if  there  is  the  slightest  discrepancy  in  Hatching.  Also,  one  rule 
aay  activate  several  other  rules,  and  the  propagation  of  rules  aay  proceed 
either  forward  or  backMard  or  a  combination  of  both. 

This  forMard  or  backMard  propagation  of  rules  describes  the  tMO  types  of 
inference  mechanisms.  The  first,  known  as  forMard  chaining,  is  simply  a 
bottoa-up  seat ch.  Starting  Mith  low-level  facts,  inforaation  is  deduced  until 
a  final  conclusion  is  reached. 


Example:  Rule  i 

IF: 

A 

THEN: 

B 

Buie  2 

IF: 

B 

THEN: 

C 

Given  A  is  true,  then  the  conclusion  C  is  reached.  The  second  method, 
known  as  backward  chaining,  is  a  goal-oriented  search.  Starting  with  an 
astuaed  final  conclusion,  the  systea  works  backward  through  the  rules  until  it 
finds  sufficient  data  to  establish  the  conclusion  or  goal. 


Example)  Find  out  about  C 

lFi  B  THEN:  C 

IF:  A  THEN:  B 

IF:  A  THEN:  C 


(goal  > 

(Buie  2) 

(Buie  D 
(Implicit  Buie) 


For  this  study,  various  inference  mechanises  were  investigated  for 
possible  use.  In  particular,  the  inference  mechanise  in  SPERll-1  (Ref.  17t 
was  thoroughly  examined.  The  inference  mechanise  used  in  SPEBIL-l  is  a 
backward  chaining  mechanism  using  the  Dempster-Shaf er  theory  of  evidence  to 
process  uncertainties  in  the  knowledge  base.  After  careful  analysis  it  was 
determined  that  an  expert  system  shell  would  be  used.  An  expert  system  shell 
is  essentially  a  complete  expert  systea  minus  the  knowledge  base.  The  shell 


chossn  for  this  project  is  called  EXSYS  and  is  available  for  use  with  personal 
computers.  A  brief  outline  of  the  features  and  capabilities  available  in 
EXSYS  are  presented  in  Appendix  B. 

UNCERTAINTY  AND  THE  USE  OF  FUZZY  SETS 

Fuzzy  Set  Theory  as  an  Approximate  Reasoning  Tool 

Uncertainty  is  associated  with  any  reasoning  system  and  it  comes  from  a 
variety  of  sources.  Reliability  of  the  information,  ambiguity  and  vagueness 
within  the  representative  language,  incompleteness  of  the  information,  and 
imprecision  in  aggregation  of  the  information  from  multiple  sources  are  but  a 
few  examples.  Zadeh  (Ref.  IS)  maintained  that  the  ability  of  humans  to 
extract  the  important  items  out  of  a  mass  of  information  is  derived  from  a 
human  tendency  to  think  approximately.  In  other  words,  human  thinking 
involves  a  gradual  transition  from  membership  in  a  class  to  non-membership  in 
a  class  rather  than  an  abrupt  change  between  classes. 

For  example,  most  peoole  would  agree  that  a  person  6’-6“  in  height  is 
"tall."  They  would  probably  also  agree  that  a  person  5’-6"  in  height  is  “not 
tall."  How  then  do  we  classify  a  person  who  is  5 ’ - 1 1 “  in  height?  Obviously, 
the  class  of  "tall"  people  is  not  a  "crisp"  set,  and,  thus,  there  is  a  gradual 
transition  from  being  a  member  of  the  class  and  not  being  a  member  of  the 
"tall"  class.  This  is  naturally  how  humans  think. 

Because  most  of  the  knowledge  in  an  expert  system  is  obtained  from  human 
experts  and  because  much  of  human  language  and  knowledge  is  vague,  it  is 
usually  true  that  facts  and  rules  are  neither  totally  certain  nor  totally 
consistent.  Because  of  this,  we  may  describe  the  reasoning  process  used  by 
experts  in  certain  situations  as  approximate.  In  this  report,  the  tneory  of 
fuzzy  sets  is  used  to  help  assess  uncertain  information  derived  from  this 
approximate  reasoni  '•  process. 

In  any  given  language,  the  values  of  a  linguistic  variable  are  words, 
phrases,  or  sentences.  For  example,  structural  damage  can  be  considered  as  a 
linguistic  variable  with  values  such  as  "severely  damaged,"  or  "moderately 
damaged."  These  are  meaningful  classifications  but  not  clearly  defined.  With 
the  use  of  fuzzy  sets,  however,  we  can  quantify  such  terminology  and  apply  it 
in  a  meaningful  way  to  help  solve  a  complex  problem.  An  evident  advantage  of 


the  fuzzy  set  approach  is  the  possibility  of  representing  numeric  and 
linguistic  variables  in  a  uniform  May  and  of  using  a  formalized  calculus  to 
manipulate  these  variables.  Simply  defined,  fuzzy  set  theory  is  a  theory  of 
the  mathematics  of  approximate  reasoning. 

Although  the  process  of  fuzzy  reasoning  is  expressed  using  mathematical 
equations,  it  is  not  a  statistical  method.  It  is  an  approximate  reasoning 
process  which  is  compatible  with  human  intuition..  The  advantage  of  the  fuzzy 
reasoning  is  that  it  can  yield  an  approximate  answer  even  when  there  is  not 
sufficient  information  to  support  a  statistical  method.  The  notion  of 
probability  stems  from,  and  depends  on,  the  idea  of  repeated  trials.  Under 
identical  and  repeatable  laboratory  conditions  conducted  on  simple  models, 
this  probabilistic  notion  readily  applies;  but,  in  real-world  (in  particular 
human)  systems,  experiments  are  rarely  identical  and  repestable.  Therefore, 
for  the  subjective  assessment  of  complex  systems,  probability  has  its 
limitations. 


Fuzzy  Sets  Defined  By  Example 


The  vagueness  or  uncertainty  as  to  whether  an  object  belongs  to  a  class 
or  set  is  a  question  of  membership,  In  classical  set  theory,  an  element  is 
either  within  the  domain  of  the  set,  or  it  is  not.  Mathematical ly,  this 
binary  notion  of  set  membership  is  handled  with  the  indicator  function.  In 
fuzzy  set  tnecry,  the  degree  of  meaberthip  of  an  element  x  in  a  set  A,  denoted 
can  be  any  value  in  the  interval  C0,l),  For  instance,  If  the 
membership  level  is  one,  then  the  item  or  object  is  definitely  a  member.  If 
the  membership  level  Is  zero,  then  the  item  is  definitely  not  a  member.  If 
the  membership  value  is  between  0  and  1,  then  the  value  stated  indicates  the 
belief  that  the  object  is  a  member  of  the  set. 

In  fuzzy  set  theory,  the  set  A  can  be  represented  in  terms  of  its 
membership  function  as  follows. 


Wft  I  *  Wa  ^  *2  1 

A  =  H  "  1  ♦ 
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where  each  x^  is  an  element  of  the  set  A.  When  x  is  a  continuous  variable, 
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the  set  A  is  denoted; 


a  .  /  jV!1 

x 

where  the  symbol  “  is  a  delimiter  which  denotes  the  association  of  the 

membership  value  y^(xj)  with  the  element  of  A,  and  the  symbols  /  and  "+" 
denote  the  union  of  all  elements  of  the  fuzzy  subset  in  the  continuous  and 
discrete  case,  respectively.  As  an  example,  suppose  the  set  A  represents  tne 
universe  of  discrete  concrete  strengths; 

A  =  [2.0,  3.0,  4.0,  5.0,  6.0,  7.03  (k/in2) 


Then,  for  this  example,  the  x.  are  the  discrete  values  of  concrete  strength. 
A  moderate  strength  concrete  may  be  expressed  in  fuzzy  terms  as: 


“Moderate  Strength" 


f  0.0  0.2  0.8 

L  2.0  +  3.0  +  4.0 


U0  CLi  0^0  I 
5.0  6.0  7.0  J 


In  other  words,  this  expression  means  that  5.0  k/in2  concrete  is  definitely  a 
member  of  the  set  “moderate  strength,"  2.0  and  7.0  k/in2  are  definitely  not 
members  of  the  set  “moderate  strength,"  and  3.0,  4.0,  and  6.0  k/in2  concrete 
are  somewhere  in  between  in  terms  of  membership.  A  classical  reoresentation 
using  crisp  set  theory  could  be  expressed  as: 


“Moderate  Strength*  » 


[ 


0 

2.0 


3.0  4.0 


1 

5.0 


1 

6.0 


JL. 

7.0 


] 


The  difference  between  a  “fuzzy"  representation  and  the  “classical" 
representation  of  “moderate  strength*  is  illustrated  in  Figure  15. 

The  construction  of  a  membership  function  is  a  strai ghtf orward  task  that 
can  be  accomplished  with  the  cooperation  and  assistance  of  a  panel  of 
experienced  engineers  or  experts  in  the  specific  field.  The  resulting 
membership  functions  can  then  be  manipulated  in  a  logical  manner  to  obtain  an 
answer  to  a  complex  problem. 

The  presentation  of  the  mathematics  of  the  fuzzy  set  theory  is  beyond  the 
scope  of  this  report.  For  more  detail*  the  reader  is  referred  to  any  of 
References  9*18,19,20,21,22,  and  23. 
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Multiple  Attribute  Decision  Analysis 


In  most  engineering  problems,  there  are  typically  many  solutions  or 
alternatives  available  to  the  engineer  lor  solving  a  problem.  In  order  to 
choose  the  most  optimum  solution  to  the  problem,  the  engineer  must  have  a 
method  by  which  to  assess  and  compare  each  of  his  alternatives.  One  such 
method  is  outlined  here.  The  usefulness  of  this  analysis  procedure  will  be 
explained  in  a  subsequent  section. 

First,  we  denote  the  alternatives  available  to  the  decision  maker  as  A  , 

A.,. ..A  .  For  each  of  these  alternatives  there  may  be  several  criteria 
2’  n 

(relevant  to  each  alternative)  available  lor  evaluation.  Each  ol  these 
criteria,  in  turn,  are  given  some  degree  ol  linguistic  rating,  denoted  r.^, 
such  as  “bad,"  “good,"  “very  good,"  etc.  Furthermore,  each  ol  these  criteria 
has  an  importance  factor  ( i . e . ,  not  all  criteria  are  equally  important).  We 
can  represent  the  importance  ol  a  criterion  with  a  linguistic  weighting  term, 
w..  The  relationship  among  these  linguistic  parameters  can  be  seen  in  Table 
ll. 


Strength  [ksi] 

Figure  15.  Classical  and  fuzzy  representations  of  a  moderate 
strength  concrete. 
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Assuming  that  the  linguistic  terms  r.^,  and  w.,  can  be  quantified,  the 
following  equation. can  be  used  to  determine  a  weighted  average  rating  for  each 
of  the  alternatives: 


Rating  for  Aj 


(4) 


where  n  equals  the  number  of  criteria.  Finally,  by  calculating  this  rating 
for  each  of  the  alternatives,  the  relative  merits  of  each  may  be  compared  to 
help  produce  an  optimum  decision. 

Dong  et  al.  (Ref.  20),  have  derived  an  efficient  computational  method, 
called  the  DSW  algorithm,  by  which  fuzzy  sets  are  integrated  into  this 
multiple-attribute  decision  process.  The  linguistic  terms  associated  with  the 
ratings  (r^.),  and  the  importance  weightings  (w^ )  are  assigned  fuzzy  set 
representations.  The  calculation  of  ratings  for  each  alternative  is  done  with 
Equation  4,  using  fuzzy  logic  equivalents  of  addition,  subtraction, 
multiplication,  and  division.  The  algorithm  is  well  suited  for  computer 
implementation. 

The  usefulness  of  this  algorithm  will  be  examined  later  where  it  will  be 
shown  how  the  DSW  algorithm  can  be  effectively  used  to  combine  the  subjective 
information  outlined  in  Section  III  into  a  preliminary  decision  relating  modes 
of  structural  response  to  degrees  of  damage. 


ORGANIZATIONAL  STRUCTURE 


Overview 

In  the  introduction  to  this  report,  it  was  stated  that  the  objective  of 
this  work  was  to  study  and  develop  a  computer  code  to  be  used  for  the  damage 
assessment  of  protective  structures.  In  Sections  1,  II,  and  III  the  problem 
was  introduced,  and  the  available  information  for  possible  inclusion  into  tne 
code  was  reviewed.  In  earlier  subsections  the  "tools'1  used  to  synthesize  all 
this  information  were  discussed  briefly.  The  purpose  of  later  subsections 
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Mill  be  to  show  how  all  this. material  is  combined  to  create  the  present' 
version  of  the  DAPS  code  (Damage  Assessment  of  Protective  Structures). 

As  discussed  previously,  the  major  factors  influencing  an  assessment  of 
the  buried  box  structure  included  structural  integrity,  functionality,  and 
repairabi 1 ity.  Given  sufficient  information,  it  would  be  possible  to 
determine  three  separate  damage  level  assessments  of  the  structure,  one  from 
each  factor.  If  each  of  these  assessments  is  then  given  a  weighting  or 
importance  value  by  a  decision  maker,  a  final  analysis  of  the  level  of  damage 
to  the  structure  can  be  computed  using  a  method  similar  to  the  one  outlined 
previously.  The  ultimate  goal,  therefore,  is  to  take  the  assessment  from  each 
of  these  factors  and  to  combine  them  into  one  meaningful  assessment.  This 
procedure  is  represented  in  Figure  16. 

The  organizational  structure  developed  to  accomplish  this  task  is  shown 
in  Figure  17.  The  level  of  damage  from  each  of  the  three  factors  can  be 
provided  using  an  expert  system  approach.  In  this  particular  case,  each  of 
the  factors  could  be  developed  as  a  separate  module  within  the  expert  system 
architecture.  Using  the  expert  system  shell  EXSYS  as  the  inference  mechanism 
and  knowledge  base,  assessments  for  structural  integrity,  functionality,  and 
repairability  can  be  calculated,  and  the  information  shared  and  stored.  With 
analyses  completed,  the  information  in  the  form  of  fuzzy  sets  is  passed  to  an 
external  program,  where  a  DSW  algorithm  is  activated  to  combine  the 
assessments.  Only  the  structural  integrity  module  of  Figure  17  was  developed 
for  this  report. 

The  logical  flow  of  the  structural  integrity  module  is  outlined  in  Figure 
IB,  There  are  three  major  parts  to  this  module  as  follows)  (!)  Rules 
concerning  total  failure,  (2)  Rules  concerning  soft  (subjective)  data,  and  (3) 
Rules  concerning  hard  (objective)  data.  The  first  part  of  the  code  deals  with 
the  trivial  case  of  total  failure  of  the  structure,  i.e.,  roof  collapse.  For 
this  particular  case,  several  rules  were  developed  from  the  questionnaire  sent 
to  the  experts  (rules  1-9  in  Appendix  C>.  If  total  failure  has  occurred, 
these  rules  are  activated,  and  the  mode  of  failure  is  determined.  The  level 
of  damage  is  set  to  total  failure,  and  the  remainder  of  the  structural 
integrity  module  is  bypassed.  An  example  rule  follows. 
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(a)  Damage  level  assessments: 

Structural  Integrity 
Lb]3x1  »  Functionality 

Repalrablllty 

(b)  Weighting  or  Importance  matrix: 


CA]ix3  ' 

(c)  Ultimate  goal: 

Final  answer: 

CC3W  •  Overall  damage 


Figure  16.  Representation  of  structural  damage  assessment  goal. 


CA]  x  [Q3  •  CCD 
Damage  level 


Damage  level 
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Expert  system  shell --EXSYS 


Figure  17.  Organizational  structure  of  DAPS. 
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Figure  18.  Logical  flow  of  structural  integrity  module. 
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RULE  NUMBER:  S 


IF: 

comolete  failure  of  the  structure  is  true  ti.e.,  the  roof  slab  has  been 
completely  separated  from  one  or  both  of  its  supports) 
and  the  failure  has  occurred  at  the  slab-wall  connection 
and  the  separated  slab  is  lying  flat  on  the  floor  (implying  a  symmetric 
failu  :) 

and  the  failure  surface(s)  at  both  supports  are  relatively  rough  with  many 
cracks  and  "concrete  teeth" 

and  inspection  of  the  failure  region  indicates  the  main  reinforcing  bars 
exhibited  rupture  after  significant  deformation  (note  “necking"  or 
stretching) 

THEN: 

damage  level  as  determined  by  structural  integrity  analysis  from  visual 
information  is  total  failure  -  Probabi 1 i ty=l 
and  the  rnode(s)  involved  in  the  deformation  process  as  determined  by  visual 
(subjective)  information  were  predominantly  shear  and  diagonal  tension, 
causing  rupture  of  reinforcement  after  significant  rotation  and 
deformation.  It  is  possible  that  the  failure  region  was  underreinforced 
-  Probabi 1 ity=l 

and  [FAILURE  MODES!  IS  GIVEN  THE  VALUE  "KNOWN" 


The  second  and  third  parts  of  the  module  deal  with  the  subjective  and 
objective  data  respectively!  as  outlined  in  Section  III.  Because  these  two 
parts  encoapass  the  eajor  effort  of  this  reporti  the  methodology  used  to 
develop  and  implement  these  two  parts  will  be  discussed  in  later  sections. 

Structural  Integrity  Analysis  Usioa.  So_f.t_D.aU 

Derivation  of  Analysis  Procedure.  When  the  data  from  the  expert 
questionnaires  were  analysed*  it  was  originally  intended  that  this  information 
would  be  codified  into  the  computer  in  typical  expert  system  fashion,  i.e..  in 
the  form  of  IF-THEN  production  rules.  Upon  further  investigation,  however,  it 
was  found  that  this  could  not  be  accomplished  easily.  The  reason  for  this  is 
due  to  a  combinatorial  explosion.  In  most  expert  systems,  the  whole  purpose 
of  a  rule  base  approach  is  to  use  intuitive  logic  in  the  fora  of  production 
rules  to  "prune"  the  number  of  paths  or  combinations  available  for  solving  a 
complex  problem.  In  the  case  under  study,  the  number  of  possible  structural 
attribute  combinations  (and  therefore  number  of  rules)  is  simply  much  too 
large  to  be  handled  in  the  normal  manner. 
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To  illustrate  this,  assume  that  each  of  the  structural  attributes  listed 
in  Table  9  can  attain  one  of  six  degrees  of  damage;  none,  slight,  moderate, 
severe,  very  extensive,  or  unknown.  Since  there  are  10  attributes,  the  number 
of  possible  combinations  for  the  full  space  is  6  =  60,466,176.  Obviously, 

even  though  all  these  combinations  are  mathematically  possible,  most  of  them 
would  have  no  real  physical  meaning  in  terms  of  describing  the  modes  of 
response  of  the  structure.  The  point  of  the  illustration  is  this;  the 
complexity  is  such  that,  to  describe  all  possible  responses,  the  number  of 
rules  required  renders  the  problem  intractable. 

Although  this  large  attribute-damage  space  could  be  pruned  using 
production  rules,  the  expertise  needed  to  develop  these  rules  is  not 
available.  HoNever,  it  is  possible  to  circumvent  this  problem  using  a  fuzzy 
logic  approach.  To  accomplish  this,  me  found  that  there  existed  two  very 
useful  conceptual  relationships  Mithin  the  expert  data  Mhich  could  be 
expressed  in  matrix  form.  The  first  matrix,  Mhich  Mill  be  denoted  the  user 
matrix  [111,  is  shown  in  Figure  19a.  The  ten  roes  represent  the  ten  structural 
attributes  identified  in  Table  9.  The  ten  slots  in  this  matrix  are  filled  by 
user  response  during  a  session  as  follous.  Given  some  type  of  visual  data 
(photographs,  visit  to  site,  etc.)  of  a  damaged  structure,  the  user  chooses 
one  of  six  degrees  (none,  slight, .. .unknown)  that  best  describes  the  level  of 
damage  of  each  structural  attribute.  For  example,  A^  might  be  slight,  Aj 
might  be  severe,  Aj  might  be  none,  and  so  on.  For  each  linguistic  term,  there 
is  a  corresponding  fuzzy  set  representation. 

The  second  matrix  relationship,  denoted  (El  in  Figure  19b,  is  slightly 
more  complex  and,  thus,  must  be  carefully  explainedm  This  matrix  is  termed 
the  expert  matrix  because  the  relational  contents  are  derived  from  expert 
opinion.  In  Figure  19b,  the  (El  matrix  is  shOMn  to  be  three-dimensional. 

Each  roe  represents  a  major  mode  of  response,  each  column  a  structural 
attribute,  and  the  third  dimension  normal  to  the  page  represents  the  six 
degrees  of  damage  attainable  by  each  structural  attribute.  Thus,  each  slot  in 
the  (A,  x  hode)  space  represents  the  relationship  between  a  given  mode  and  a 
given  attribute  Mhich  has  achieved  a  given  degree  of  damage.  Careful  study  of 
the  expert  reasoning  process  showed  that  most  experts’  judgments  of  the 
relationship  between  a  mode  and  an  attrioute  depend  on  the  degree  of  damage 
attained  by  that  attribute. 
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(a)  [U]  Matrix  produced  by  user  Input: 


Damage  level 


A,  *  jth  Structural 
attrl bute 


(b)  [E]  Matrix  produced  by  expert  opinion: 

Unknown 

Very  extensive - 

Severe 

iterate  — 

Shear 

Diagonal  Tension 
Flexure 
Tension  Membrane 


(c)  Product  of  C£J  x  CU] 


Shear 

Diagonal  tension 
FI  exure 
Tension  Membrane 


Figure  19.  Matrix  relationships  for  subjective  information. 


For  instance,  take  the  first  slot  in  the  upper  left-hand  corner  of  Figure 
19b.  Given  that  attribute  Aj  is  none  (with  respect  to  damage),  an  expert 
might  determine  that  there  is  no  relationship  between  shear  and  attribute  A  . 
Obviously,  the  relationship  can  vary  between  no  relationship  (i.e.,  attribute 
Aj  has  nothing  to  do  with  a  shear  mode  response)  and  complete  relationship 
(i.e.,  shear  mode  response  is  significantly  dependent  on  the  degree  of 
attribute  Aj).  Since  the  idea  of  relationship  is  a  continuous  function 
between  no  relationship  and  complete  relationship,  the  use  of  fuzzy  sets  to 
.represent  these  vague  terms  is  ideal.  Figure  20  shows  a  typical  fuzzy  set 
representation  of  no  relationship  for  the  example  above. 

Given  that  the  CUl  matrix  is  filled  by  user  response,  and  the  CE1  matrix 
is  filled  by  expert  opinion,  the  two  may  now  be  combined  to  produce  an 
analysis  of  the  damaged  box.  The  process  used  to  combine  the  two  matrices  is 
the  OSH  algorithm  discussed  previously.  In  this  case,  however,  the  criteria 
are  the  structural  attributes,  the  ratings  are  the  values  given  in  (Ul,  and 
the  weightings  are  the  values  given  in  (El.  Using  the  DSU  algorithm,  the 


Degree  of  asiociatlon/reUtton  between 
•  given  nodi  and  a  glean  attribute 


Figure  20.  Example  fuzzy  set  representing  "No  Relationship" 
between  a  shear  mode  and  structural  attribute  . 


damage  level  of  t.he  structure  associated  with  each  of  the  four  modes  can  be 
calculated  using  Equation  5. 


(Damage  Level ) . 


(5) 


It  is  important  to  remark  here  that  although  the  [El  matrix  is  three- 
dimensional  conceptually,  when  used  in  the  calculation  process,  it  has  been 
reduced  to  a  two-dimensional  matrix.  This  reduction  process  results  when  the 
user  chooses  one  of  the  six  levels  of  damage  for  each  of  the  structural 
attributes.  By  doing  so,  the  other  five  columns  (in  the  third  dimension  of 
Figure  19b)  for  each  attribute  are  eliminated.  In  other  words,  only  one 
column  in  the  third  dimension  is  used  for  each  attribute. 

Once  CE1  has  been  reduced  to  a  two-dimensional  matrix,  it  may  be  combined 
with  1U].  The  result  of  this  operation  is  a  4x1  matrix  denoted  [PI,  as  shown 
in  Figure  19c.  Each  slot  of  CPI  is  a  fuzzy  set  describing  the  level  of  damage 
associated  with  the  given  mode  of  damage. 

Before  continuing  on  with  an  explanation  of  how  the  results  of  this 
process  are  interpreted,  an  explanation  of  how  fuzzy  sets  are  chosen  for 
linguistic  terms  is  necessary.  Because  the  linguistic  terms  themselves  are 
vague,  the  choice  of  shape(s)  used  to  describe  these  terms  is  also  rather 
arbitrary.  In  the  literature,  the  shapes  most  often  chosen  include  the 
triangle  and  the  trapezoid.  Shapes  with  other  than  straight  lines  may  be 
chosen,  but  com;  lexity  in  mathematical  manipulation  makes  them  very 
undesirable.  F  r  this  study,  the  fuzzy  sets  chosen  to  describe  degrees  of 
damage  are  shown  in  Figure  21.  Typically,  the  shapes  and  their  locations  in 
the  universe  of  discourse  (abscissa  axis)  are  decided  by  consensus  of  expert 
opinion. 

In  particular,  the  reader  should  notice  the  fuzzy  set  representing  the 
linguistic  term  "unknown."  After  much  debate  and  many  trial  shapes  and 
calculations  with  these  shapes,  a  rectangle  (unit  membership  throughout)  was 
chosen.  Upon  reflection,  this  shape  appears  to  be  the  most  intuitive  choice. 
The  very  idea  of  an  "unknown"  damage  level  implies  that  it  could  be  anywhere 


59 


u(x)  •  Membership 


Damage  Grade 


Figure  21.  Fuzzy  sets  representing  linguistic  damage  levels. 

between  “none"  and  "very  extensive,*  The  implication  of  including  unknown 
into  a  calculation  it  to  extend  the  rtiulting  answer  in  both  directions,  j.e., 
toward  “none'  and  toward  “very  extensive.*  Thut,  tht  addition  of  “unknown* 
into  an  analytit  ltndt  tort  uncertainty  to  th»  antwtr,  juat  at  it  mould. 

Tht  inttrprttation  of  tht  calculated  fuzzy  sett  it  a  difficult  tatk. 
Ideally,  the  goal  it  to  find  a  linggiatic  tera  that  ccrretpondt  to  each 
calculated  fuzzy  tat.  Bocause  of  the  nature  of  the  weighted  average 
algorithm,  however,  the  actual  thapes  of  the  reaulting  fuzzy  tets  aay  not  he 
at  neat  and  tyaaetrical  as  the  fuzzy  eett  used  to  produce  thee  (l.e., 
triangles,  trapezoids,  or  delta-functional.  If  we  tuperiepote  a  typical 
calculated  fuzzy  set  onto  the  universe  of  fuzzy  sett  representing  the 
different  damage  levels  at  shown  in  Figure  22,  it  it  possible  to  make  a 
comparisun  between  the  shape  of  the  calculated  set  and  each  of  the  predefined 
damage  level  sett.  One  way  of  coeparinq  fuzzy  sett  it  to  use  a  difference 
measure,  denoted  D.  If  the  value  of  0  between  the  calculated  fuzzy  set  and 
one  of  the  predefined  sett  it  small,  then  it  it  possible  to  tay  the  sett  are 
timilar.  In  the  cate  of  good  similarity,  th«  linguistic  term  for  the 
predefined  fuzzy  set  aay  be  applied  to  the  calculated  fuzzy  tet. 
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The  difference  D  can  be- measured  by  a  number  of  techniques,  some  of  which 
were  examined  for  this  study.  The  difference  measure  chosen  is  a  modified 
Euclidean  approach  and  is  given  by  Equation  h.  This  measure  was  chosen 
because  it  was  most  easily  adapted  to  the  DSW  algorithm,  and  because  it 
provides  a  relatively  good  measure  of  similarity  between  a  predefined  fuzzy 
set  and  a  calculated  fuzzy  set.  In  this  context,  similarity  is  meant  to 
describe  both  "similarity  of  shape"  and  "distance  between"  the  predefined  and 
calculated  fuzzy  sets. 


[(X’ 


R  *i 


(6) 


where 

n  »  number  of  a-cuts 

i  “  1,2,  •  •  *10 

(XM.  »  left-hand  damage  level  value  at  the  ith  a-cut  level  of  B’ 

~  *  1L 

(X. 1.  =  left-hand  damage  level  value  at  the  il  a-cut  level  of  B. 

V  •  4  u  * 

(Xii ,  »  Right-hand  damage  level  value  at  the  i  a-cut  level  of  B' 

^  *  4  u 

< XR ) ^  •  Right-hand  damage  level  value  at  the  i  a-cut  level  of  Bi 

The  modification  to  the  Euclidean  approach  involves  the  use  of  the  a-cut 
interval  method,  which  forms  the  basis  of  the  OSN  algorithm.  As  we  see  in 
Figure  22,  a-cuts  are  simply  parallel  horizontal  lines  of  constant  membership 
that  cut  through  the  fuzzy  sets.  The  intersection  of  these  lines  with  a  fuzzy 
set  identify  points  on  the  universe  of  discourse  which  have  the  same 
membership  value. 

Using  Equation  6,  it  is  possible  to  calculate  the  difference  measure,  0, 
which  is  a  measure  of  the  difference  between  the  calculated  fuzzy  set  and  eacn 
predefined  fuzzy  set.  The  linguistic  term  chosen  for  the  calculated  fuzzy  set 
is  simply  the  linguistic  term  corresponding  to  the  predefined  fuzzy  set  with 
the  lowest  0  ti.e.,  highest  similarity).  Although  this  is  probably  not  the 
most  rigorous  or  elegant  method  available  for  determining  similarity,  for 
demonstration  purposes,  it  provides  a  good  indication  of  relative  similarity. 


62 


Explanation  of  External  Proarams  Used  in  Soft  Data  Analysis 


The  previous  discussion  Mas  -focused  on  describing  the  methodology  used  in 
combining  subjective  expert  opinion  with  observed  structural  damage  provided 
by  the  user  of  the  code.  In  this  section,  the  focus  will  be  on  presenting  the 
manner  in  which  this  information  was  coded  into  the  computer.  As  was  shown  in 
Figure  18-,  the  code  begins  by  determining  if  complete  failure  of  the  structure 
has  occurred.  If  true,  the  code  proceeds  as  discussed  previously.  However, 
if  the  answer  is  false,  a  series  of  calls  to  external  programs  is  made  by 
EXSYS  in  order  to  begin  processing  the  soft  data.  There  are  four  external 
programs  used  in  the  soft  data  manipulation  process,  as  shown  in  Figure  23. 

All  external  programs  were  written  in  either  FORTRAN  or  BASIC  computer 
languages.  Programs  written  in  BASIC  were  done  solely  because  of  the 
extensive  graphics  capabilities  available  in  the  language.  Each  of  the  four 
programs  will  be  discussed  next. 

Prooram  FUZSET.  The  program  FUZSET  is  the  first  program  to  be  called 
from  EXSYS.  This  portion  of  the  code  is  written  in  BASIC  and  is  composed  of 
two  major  sections.  The  first  section  of  the  program  presents  graphic 
pictures  of  the  fuzzy  sets  that  Mill  represent  the  linguistic  terms  none, 
slight,  moderate,  severe,  very  extensive,  and  unknown,  similar  to  Figure  21. 
These  fuzzy  sets  were  predetermined  by  the  authors  and  are  stored  in  the  data 
file  USER. OAT.  During  the  pr< an  tat  ion  of  these  fuzzy  sets,  the  user  js  given 
the  option  of  viewing  or  changing  their  shapes  if  he/she  is  uncomfortable  with 
the  current  defined  shapi 

The  second  section  of  the  program  deals  with  the  user  matrix  IU3  and  the 
expert  matrix  C E 3  described  previously.  It  is  at  this  point  in  the  code  that 
the  comouter  questions  the  user  about  the  ten  structural  attributes  listed  in 
Table  9  (assuming  they  have  access  to  the  information).  During  this  Question 
and  answer  period,  the  user  also  has  the  option  of  viewing  the  fuzzy 
relationships  in  the  expert  matrix  C E 3  dealing  with  the  attribute  about  which 
they  are  being  questioned.  The  purpose  of  this  is  that  it  may  aid  them  in 
their  decision  making  process.  As  an  added  option,  the  user  may  also  change 
the  relationships  in  the  C E 3  matrix  if  they  are  unsatisfied  with  expert 
opinion  on  a  given  mode/attribute  relationship.  All  fuzzy  sets  are  stored  in 
data  files  as  shown  in  Figure  23. 
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Program  CQMBIN.  The  program  COHBIN  is  a  program  written  in  FORTRAN  and 
has  no  direct  interaction  with  the  user.  COHBIN  is  called  immediately  upon 
completion  of  the  program  FUZSET,  and  its  major  purpose  is  the  combination  of 
information  derived  from  the  user  in  FUZSET.  This  combination  process  is 
accomplished  through  the  use  of  Equation  5  which  outputs  four  fuzzy  sets. 
Again,  these  four  fuzzy  sets  represent  the  level  of  damage  of  the  structure 
associated  with  the  four  modes  of  response.  After  the  DSW  algorithm  has  been 

ft 

completed,  a  similarity  analysis  is  performed  on  each  calculated  fuzzy  set 
using  Equation  b.  Thus,  a  linguistic  term  is  applied  to  each  mode  of 
response.  Of  these  four  linguistic  descriptions,  the  term  describing  the 
highest  level  of  damage  becomes  the  overall  damage  indicator  of  the  structure 
as  determined  by  subjective  data. 

Program  COLORRLT.  The  BASIC  program  CQLORPLT  is  run  by  EXSYS  after  all 
calculations  in  COHBIN  have  been  completed.  The  function  of  CQLORPLT  is  quite 
simple  —  to  present  to  the  user  a  graphical  representation  of  the  results 
obtained  in  the  program  COHBIN.  The  user  is  asked  if  he/she  would  like  to 
view  the  results.  If  yes,  the  program  is  rum  if  no,  EXSYS  continues  without 
executing  CQLORPLT. 

Prooram  RETURN.  As  seen  in  Figure  24,  the  sole  function  of  program 
RETURN  (written  in  FORTRAN)  is  to  transfer  to  EXSYS  answers  obtained  in  all 
external  programs.  The  values  obtained  in  COHBIN  are  placed  in  a  temporary 
data  file  called  TEHPl.DAT,  RETURN  simply  reads  these  values  and  returns  them 
to  EXSYS  via  file  RETURN. OAT,  as  explained  in  Appendix  B.  Five  values  are 
returned  to  EXSYS  in  this  manner,  the  first  four  representing  the  damage 
associated  with  each  mode,  and  the  fifth  representing  the  overall  damage  of 
the  structure. 

Once  the  values  have  been  passed  to  EXSYS,  there  are  several  rules  (rules 
13  to  42  in  Appendix  C)  which  provide  a  linguistic  interpretation  of  the  modes 
and  damage  levels  computed  by  the  program  COHBIN.  An  example  rule  follows. 
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Via  File 
Return. OAT 


Return  Values  Calculated 
In  TIM  and  IFPLOT 


From 


To  EXSYS 


Via  File 
I  Return .OAT 


Return  Values  Calculated 
In  DEFLECT 


Figure  24.  Organization  of  hard  data  external  functions. 
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RULE  NUMBER:  27 


IF: 

complete  Failure  oF  the 'structure  is  False  (i.e.,  the  rooF  slab  has  not 
been  completely  separated  From  one  or  both  oF  its  supports) 
and  Flexure  damage  is  moderate 

THEN: 

the  damage  to  the  structure  associated  with  a  Flexure  response  as 
determined  by  structural  integrity  analysis  From  visual  inFormation  is 
moderate  -  Probability*! 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  Fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
oF  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  oF 
Failure  to  these  same  parameters. 


This  inFormation  not  only  is  output  at  the  end  oF  an  EXSYS  run  but  can  be  used 
later  in  the  code  as  a  priori  knowledge  in  rules  dealing  with  the 
Functionality  and  repairability  oF  the  structure. 

St  rue  tun  1  _  integrity  ..An  a  1  v » t  s.  m  i  n  a  H  tr  d  _P  a  t  a 

Having  completed  the  analysis  oF  the  soFt  data,  the  code  proceeds  to  the 
hard  data  portion,  as  shown  in  Figure  10.  In  the  soFt  data  analysis  just 
discussed,  the  inFormation  used  was  derived  From  "extremely  late-time*  data, 
i.e.,  posttest  observations.  On  the  other  hand,  all  available  data  to  be  used 
in  the  hard  data  analysis  corresponds  to  structural  responses  occurring 
between  0  and  IS  ms  aFter  initial  loading.  Nevertheless,  the  hard  data 
analysis  described  in  this  section  will  be  assigned  to  two  segments,  one  based 
on  early-time  analysis  and  one  on  late-time  analysis. 

The  reasoning  For  these  separate  sections  is  based  on  the  previous 
discussion  oF  the  hard  data.  In  that  section,  it  was  pointed  out  that  not  aU 
data  are  oF  equal  value  in  determining  modes  qF  response  at  diFFerent  points 


67 


in  tine.  Interface  pressures,  for  example,  give  a  good  indication  of  early- 
time  (less  than  approximately  2  ms)  slab  response,  whereas  deflection  profiles 
give  a  good  indication  of  late-time  (greater  than  approximately  2  as)  slab 
response.  The  early-time  analysis  is  important  because  it  provides  an 
indication  of  trends  in  response.  Thus  by  knowing  the  early  response  first, 
late-time  analysis  (from  deflection  profiles)  and  overall  response  modes  may 
be  better  understood  and  interpreted. 

Therefore,  in  the  next  subsection  the  techniques  and  external  programs 
used  to  calculate  “early-time“  response  will  be  discussed.  In  the  following 
subsections  the  methodology  used  to  provide  "late-time"  analysis  is  discussed, 
as  well  as  the  rules  used  to  combine  both  analyses  into  an  overall  "hard" 
assessment  of  the  structure. 

Early  Time  Analysis 

The  calculation  of  the  early-time  response  of  the  structure  in  DAPS  is 
limited  to  two  programs  external  to  EXSYS.  These  two  programs,  TIM  and 
IFPLOT ,  are  illustrated  in  the  top  half  of  Figure  24,  and  an  explanation  of 
each  is  given  below. 

Prooram  TIM.  Program  TIM  is  a  FORTRAN  code  developed  by  Mickelsen  (Ref. 
24),  which  determines  the  early  time  response  of  a  fixed-fixed  beam  using  a 
Timoshenko  beam  model.  The  computer  program  first  solves  for  the  natural 
frequencies  of  the  Timoshenko  beam  model.  These  natural  frequencies  are  then 
used  in  a  modal  analysis  to  obtain  the  response  characteristics  of  the  bean. 
The  purpose  at  the  program  is  to  determine  the  controlling  early  time  response 
mode  of  the  structural  element.  This  is  accomplished  by  determining  the  times 
at  which  the  support  shear  and  support  bending  moment  reach  their  ultimate 
capacities.  Whichever  force  exceeds  its  capacity  first,  controls  the  early 
time  response. 

The  fixed-fixed  beam  model  is  used  because  it  represents  very  well  the 
one-way,  early-time  response  of  the  roof  slab  in  the  buried  box  structure 
(Ref.  7).  The  inputs  to  the  program  include  material,  geometrical,  and 
loading  parameters  of  the  box  structure  under  study.  These  parameters  are 
obtained  interactively  from  the  user  before  the  program  TIM  is  executed.  The 
output  of  the  program  is  a  value  equal  to  either  0  or  1,  representing  a  direct 
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shear  or  flexure  aode  respectively.  The  value  is  written  to  the  teaporary 
data  file  TEMP2.DAT  for  later  use,  as  shown  in  Figure  24. 

Prooraa  I F PLOT .  Prograa  IFPLQT  is  a  BASIC  prograa  whose  function  is  to 
deteraine  the  siailarity  of  interface  pressure  records  as  discussed 
previously.  The  prograa  begins  by  showing  the  user  the  superiaposed  interface 
pressure  records  and  corresponding  iapulse  plots  of  the  structure  under  study. 
The  user  is  then  presented  with  a  question  concerning  the  siailarity  of  the 
superiaposed  plots.  Because  this  question  involves  a  large  aaount  of 
subjectivity,  the  user  is  allowed  to  view  exaeple  plots  representing  two 
possible  responses.  The  user  is  also  allowed  to  toggle  between  the  question, 
exaeple  plots,  and  actual  plots  as  desired,  before  eaking  a  decision.  As 
discussed  earlier,  the  purpose  of  this  coaparison  is  to  deteraine  whether  the 
early-tiae  response  falls  into  the  category  of  shear,  flexure,  or  shear* 
flexure.  The  output  of  IFPLOT  is  a  value  between  1  and  3  (shear  *  1,  flexure 
*  2,  shear-flexure  *  3),  representing  each  of  the  possible  aodes,  and  is 
stored  in  the  teaporary  data  file  TENP3.DAT. 

With  prograas  TIM  and  IFPLOT  coapleted,  the  data  obtained  froa  these 
analytes  are  transferred  back  to  EXSYS  via  the  prograa  RETURN.  Within  EXSYS, 
the  data  froa  TIM  and  IFPLOT  are  linguistically  interpreted  in  rules  43-46, 
providing  an  assessaent  of  earlytlae  response.  An  exaeple  rule  follows. 

RULE  NUMBER t  43 
IF* 

coaplete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  coapletely  separated  froa  one  or  both  of  its  supports) 
and  early  tiae  response  as  deterained  by  aodified  Tiaoshenko  beat  analysis  is 
shear 

and  early  tiae  response  as  deterained  by  siailarity  analysis  of  interface 
pressure  plots  is  shear 

THENt 

early  tiae  response  as  deterained  by  a  coabination  of  Tiaoshenko  beaa 
analysis  and  interface  pressure  data  is  predoainantly  shear  - 
Probability*! 


In  the  case  of  conflicting  answers  from  the  two  programs  (e.g.,  shear  and 
flexure),  the  response  is  assumed  to  be  a  combination  of  both.  Therefore, 
unless  both  TIM  and  IFRLOT  produce  the  same  answer  (i.e.,  both  shear  or  both 
flexure),  the  early-time  response  is  termed  shear-f lexure.  This  conclusion 
simply  implies  that  there  is  no  overwhelming  evidence  to  substantiate  one 
dominating  mode.  This  information  is  used  later  to  help  determine  the  late- 
time  response  of  the  structure. 

Late  Time  Analysis 

The  calculation  of  late-time  response  of  the  structure  in  DAPS  is  done  in 
the  external  program  DEFLECT,  as  shown  in  the  bottom  half  of  Figure  24.  An 
explanation  of  this  program  follows. 

Program  DEFLECT.  DEFLECT  is  a  program  written  in  BASIC.  It  shows  the 
user  the  deflection  profiles  of  the  underside  of  the  roof  slab  at  3  ms  and  IS 
as.  The  usefulness  »?  deflection  data  in  determining  modes  of  response  was 
presented  earlier  and  is  implemented  in  DEFLECT.  Only  two  deflection  profiles 
were  used  to  illustrate  the  concept  of  late-time  analysis.  Other  deflection 
profiles  at  different  times  could  have  been  used  Just  as  easily. 

After  viewing  the  deflection  profllest  the  user  is  asked  three  Questions 
(Fig.  23)  concerning  the  magnitude  of  deflection  at  the  centerline  and 
deflection  ratios  calculated  at  the  wall,  Quarter  point,  and  centerline.  To 
aid  the  decision  process,  the  user  is  allowed  to  toggle  between  questions  and 
profiles  in  order  to  study  them  more  carefully.  Because  the  questions  are 
somewhat  subjective,  numerical  limits  on  linguistic  terms  <e.g,,  'very  small 
deflection*  is  approximately  0-3  in)  are  given  to  aid  the  user  in  his/her 
response. 

The  damage  level  of  the  structure  as  determined  by  objective  information 
is  also  calculated  within  the  program  DEFLECT.  Using  Equation  3,  the 
deflection  profile  for  15  ms  is  numerically  integrated  to  orovide  a  measure  of 
overall  damage.  The  value  calculated  by  Equation  3,  along  with  the  answers  to 
the  three  questions,  are  stored  in  a  temporary  data  file,  TENP4.DAT,  as  shown 
in  Figure  24. 


Question  1 


The  deflection  (at  15  ms)  of  the  roof  slab  near  the  *811  (2  1/2  In  from 
face)  Is  5.5  in.  Typically,  when  evaluating  modes  of  response  (especially 
shear),  a  near-wall  deflection  In  the  range  of 

■  <  0.5  In  Is  considered  very  small. 

■  0.5-1. 5  In  Is  considered  small. 

*  1. 5-3.0  In  Is  considered  moderate. 

•  >3.0  In  Is  considered  large. 

In  your  opinion,  given  this  Information  and  the  deflection  profile,  which 
category  do  you  think  the  given  deflection  of  5.5  In  belongs? 

1.  Very  small 

2.  Small 

3.  Moderate 

4.  Large 

Enter  a  value  of  1-4,  or  just  ENTER  to  view  plot? 

Question  2 


At  15  ms,  the  ratio  of  differential  deflection  (near-wall  deflection  minus 
center’ine  deflection)  to  slab  thickness  Is  0.2.9.  Typically,  when 
evaluating  modes  of  response  (especially  flexure  and  tension  membrane),  a 
ratio  In  the  range  of 

•  <  0.3  is  considered  In  relatively  flat. 

•  0.3-0.6  is  considered  moderately  curved. 

•  >  0.6  is  considered  highly  curved. 

In  your  opinion,  given  this  Information  and  the  deflection  profile,  which 
category  do  you  think  the  qlven  ratio  of  0.26  belongs? 

1.  Relatively  flat 

2.  Moderately  curved 

3.  Highly  curved 

Enter  a  value  of  1*3,  or  just  ENTER  to  view  plot? 

Question  3 

At  15  »s,  the  ratio  of  differential  deflection  (near-wall  deflection  minus 
centerline  deflection)  to  slab  thickness  Is  0.20.  At  15  «s,  the  ratio  of 
differential  deflection  (near-wall  deflectJon  minus  the  deflection  a 
distance  t  «  thickness  from  the  wall)  is  0,06.  The  difference  between  these 
two  ratios  is  0.22.  Typically,  whe"  evaluating  modes  of  response 
(nsoecially  diagonal  tension),  a  difference  of 

•  <  0.15  is  considered  to  be  small. 

•  0.15*0.30  is  considered  to  be  moderate. 

•  >  0.30  is  considered  to  be  large. 

In  your  opinion,  given  this  Information  and  the  deflection  profiles,  which 
category  do  you  believe  the  difference  of  0.22  belongs? 


1.  Small 

2.  federate 

3.  Largo 


Enter  a  value  of  1—3,  Or  Just  ENTER  to  view  plot.. 


Figure  25.  Questions  in  Program  DEFLECT. 
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Combination  of  Early-  and  Late-Time  Analyses 


Upon  completion  of  the  program  DEFLECT,  the  answers  are  transferred  back 
to  EXSYS  via  program  RETURN  for  interpretation.  Within  QSYS,  rules  30-73  are 
used  to  linguistically  interpret  the  answers  obtained  from  DEFLECT,  A  typical 
example  is  given  as  rule  30  in  Figure  26.  In  rule  SO,  the  s&cond  antecedent 
deals  with  the  early-time  response  of  the  structure.  The  answer  to  this 
antecedent  (i.e.,  shear  or  not  shear)  is  predetermined  by  the  early-time 
analysis  as  explained  previously.  The  third  and  fourth  antecedents  of  the 
rule  are  determined  by  the  external  program  DEFLECT.  Thus,  after  all  external 
programs  have  been  run,  if  each  of  these  antecedents  is  true,  then  the 
consequent  portion  of  the  rule  is  fired  (i.e.,  late  time...  predominantly 
shear ) . 

This  example  shows  that,  although  the  early-tiae  response  and  the  late¬ 
time  response  are  calculated  in  separate  external  programs,  the  late-time 
response  as  determined  by  rules  50  through  73  is  dependent  on  the  combined 
responses.  Thus,  rules  50  through  73  represent  ail  possible  combinations  of 
output  from  the  external  programs,  and  the  late-tiae  response  is  actually 
determined  by  looking  at  both  early-  and  late-time  responses  together. 

As  for  the  overall  damage,  the  value  calculated  within  DEFLECT  using 
Equation  3  is  given  a  corresponding  linguistic  interpretation  in  rules  74-78. 
Rule  74  is  shown  in  Figure  26  as  an  example. 

Summary 

In  this  section,  the  methodology  used  to  organiae  a  large  problem  in  the 
form  of  an  expert  system  was  described.  In  particular,  it  was  shown  now  rule- 
based  information,  fuiry  set  theory,  and  cal culati onal  programs  with  varied 
functions  can  be  implemented  efficiently  into  a  single  system  to  help  solve  a 
complex  problem.  As  an  example,  the  programming  techniques  used  to  codify 
this  information  for  structural  integrity  analysis  of  a  buried  box  were 
described. 

The  output  of  this  structural  integrity  module  are  the  modes  of  response 
and  levels  of  damage  as  calculated  by  subjective  and  objective  information. 

At  this  point  it  is  important  to  remind  the  reader  that  the  output  from  this 
module  is  only  a  partial  answer.  This  information,  in  turn,  could  be  used  as 
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RULE  NUMBER:  50 


IF: 

complete  Failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early-time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  shear  =  1 
and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  very  small  or  small 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 
centerline  deflection)  to  slab  thickness  implies  a  relatively  flat 
surface 


THEN: 

late-time  response  as  determined  by  deflection  profiles  was  predominantly 
shear  -  Probabi 1 i ty=l 
and  CLATE  RESPONSE]  IS  GIVEN  THE  VALUE  1 


RULE  NUMBER:  74 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  soparatad  from  one  or  both  of  its  supports) 
and  [DAMAGE  LEVEL]  <  0.2 

THEN: 

the  overall  damage  to  the  structure  as  determined  by  deflection 
information  is  none  -  Probability5! 


NOTE: 

The  value  of  the  variable  [DAMAGE  LEVEL]  is  determined  in  the  external 
program  DEFLECT. FOR  and  returned  via  thm  program  RETURN. FOR. 


Figure  26.  Example  rules  for  the  combination  of  early-  and 
late-trme  analyses. 


« 
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a  priori  knowledge  in  the  nodules  on  Functionality  and  Repairabi 1 ity,  as  well 
as  being  output  at  the  end.  Again,  the  eventual  goal  is  a  combination  of 
analyses  produced  From  each  of  the  modules  in  order  to  provide  an  overall 
assessment  of  the  damaged  structure.  How  the  information  from  each  of  the 
individual  modules  is  used,  is  limited  only  by  the  imagination  of  the 
developer. 


V.  AN  EXAMPLE  SESSION 


Having  detailed  the  background  and  structure  of  the  computer  code  DAPS  in 
Sections  I  through  IV,  an  example  session  is  presented  here  in  order  to  give 
the  reader  a  thorough  understanding  of  how  it  works.  In  a  typical  case  study, 
the  specimen  used  to  test  the  prototype  is  usually  different  than  those  used 
to  develop  it.  Because  of  data  limitations,  however,  the  code  will  be 
demonstrated  by  analyzing  test  element  DS3. 

The  organization  of  this  section  is  quite  simple.  Starting  with  the 
introduction  of  the  code,  the  chronology  of  events  from  program  start  to 
program  finish  will  be  explained  using  figures  taken  from  the  computer 
monitor.  The  several  figures  presented  throughout  this  section  are  actual 
“screen  dumps"  taken  directly  from  the  computer,  and  thus  represent  exactly 
what  the  user  would  see  when  sitting  at  the  terminal. 

The  program  is  initiated  by  typing  the  words  EXSYS  DAPS  at  the  DOS  (disk 

operating  system)  command  level.  This  set  of  commands  loads  the  expert  system 
shell  EXSYS  and  all  rules  and  external  programs  needed  for  execution.  Tne 
beginning  of  a  session  with  DAPS  is  marked  by  the  appearance  on  the  screen  of 
the  title/author  block  as  shown  in  Figure  27.  Upon  pressing  any  key  on  the 
keyboard,  an  introduction  to  the  code  is  presented  to  the  user,  as  shown  in 
Figure  28.  The  introduction  provides  a  general  background  on  what  is  to  come, 
and  also  what  information  the  user  must  have  available  in  order  to  use  the 
code. 

The  screen  following  the  introduction  (Fig.  29)  contains  the  first 
question  relative  to  the  integrity  of  the  structure.  In  particular,  it  is 
trying  to  determine  if  complete  failure  of  the  structure  occurred.  At  this 

point,  it  is  important  to  note  the  text  at  the  bottom  of  the  screen  below  the 

thin  line.  This  text  (which  appears  only  when  inside  EXSYS),  is  a  menu  of 
possible  actions  that  the  user  may  take.  The  absence  of  this  information  at 
the  bottom  of  the  screen  indicates  that  the  code  has  exited  EXSYS  temporarily 
to  an  external  program.  Using  the  photographs  in  Appendix  A  as  reference,  the 
answer  to  the  question  for  test  DS3  is  2}  complete  collapse  is  false. 

Because  the  answer  was  false,  the  program  immediately  exits  EXSYS,  and 
enters  the  program  FUZSET.  Once  inside  FUZSET,  an  explanatory  screen  appears 
as  shown  in  Figure  30.  At  this  point,  the  user  is  given  the  option  of 
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Structural  damage  assessment  of 
reinforced  concrete  buried  boxes 
subjected  to  impulsive  blast  loadings. 


Figure  27.  Title  block  screen. 


Tha  following  program  is  an  initial  attampt  at  combining  linguistic 
Information  in  tha  form  of  expert  opinion  on  damaga  assassamant,  with  crisp 
numerical  data  obtained  from  instrumentation  waveforms.  This  program  is  not  an 
end  product  as  such,  but  rather  a  demonstration  tool  of  what  can  possibly  be 
achaived  using  axpart  system  techniques  in  attacking  such  a  problem.  The 
program  makes  several  calls  to  external  programs  which  taka  from  a  few  seconds 
up  to  several  minutes,  so  please  be  patient  with  tha  system. 

The  first  portion  of  the  code  deals  with  the  structural  integrity  analysis 
of  the  structure  thru  the  use  of  visual  information,  therefore,  the  user  will 
need  have  observed  or  lncpected  the  structure,  or  have  visual  information  in 
tha  form  of  photographic  data  available. 

The  second  portion  of  the  code  deals  with  structural  integrity  analysis  of 
the  structure  via  numerical  data.  Host  of  the  hard  data  has  been  stored 
Internally  in  the  form  of  data  files,  but  some  information  must  be  derived 
from  the  user.  Please  have  the  following  Information  available  relating  to  the 
material,  geometrical,  and  loading  parameters  of  the  structure  i  beam  density, 
shear  modulus,  roof  slab  clear  span  length,  slab  thickness,  poisson’s  ratio, 
damping  coefficient,  concrete  compressive  strength,  steel  yield  stress, 
percentage  steel  ratio,  maximum  slab  overpressure,  load  rise  time,  and  load 
duratl on. 


Press  any  key  to  start! 


Figure  28.  Introduction  to  DAPS. 
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complete  failure  of  the  structure  is 

1  true  (i.e:,  the  roof  slab  has  been  completely  separated  from  one  or 
both  of  its  supports) 

2  false  (i.e.,  the  roof  slab  has  not  been  completely  separated  from  one 
or  both  of  its  supports) 


Enter  number (*)  of  value(s),  WHY  for  Information  on  the  rule, 
QUIT  to  save  data  entered,  <ESC>  far  edit  mode  or  <H>  for  help 


Figure  29.  Determination  of  failure. 


In  a  moment  you  Mill  be  asked  10  questions  pertaining  to  the  structural 
integrity  of  the  buried  box  structure.  These  questions  will  deal  with  the 
DAMAGE  LEVEL  of  different  features  obtained  thru  visual  inspection  of  the 
damaged  structure.  The  possible  responses  to  these  questions  will  bei 

1. )  none 

2. )  slight 

3.  )  moderate 

4. )  severe 

3.)  very  extensive 
6.)  unknown 


For  each  of  these  linguistic  terms,  there  is  stored  In  a  data  file  a 
corresponding  fuzzy  set  representation.  These  fuzzy  sets  will  be  used 
in  a  fuzzy-weighted-average  algorithm  to  determine  modes  of  damage  and 
associated  levels  of  damage  in  each  mode.  If  you  want  to  see  these  fuzzy 
set  representations  (with  the  choice  of  changing  them  if  you  disagree), 
simply  answer  yes  to  the  following  question. 


Do  you  want  to  see  the  fuzzy  set  representati ons 
<y  *  view/change,  n  ■  continue  with  program)  ?  y 


Figure  30.  Introduction  to  FUZZSET. 
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viewing/ch an.g i n g  the  fuzzy  sets  that  Hill  represent  the  linguistic  damage 
levels  described  in  Section  IV.  Having  chosen  yes  as  the  response,  Figure  31 
appears.  These  two  screens  are  used  as  a  preview  of  the  fuzzy  sets  to  come. 

As  described  in  the  figure,  these  examples  are  given  for  comparative  purposes 
only.  The  code  then  goes  on  to  show  the  fuzzy  sets  for  each  of  the  linguistic 
terms:  none,  slight,  moderate,  severe,  very  extensive,  and  unknown.  An 
example  showing  the  screen  for  "moderate  damage"  is  shown  in  Figure  32. 

The  major  goal  of  program  FUZSET  is  the  determination  of  structural 
damage  as  obtained  from  subjective  information.  In  order  to  complete  this 
task,  the  program  begins  the  interrogation  of  the  user  concerning  the  damage 
level  of  the  ten  structural  attributes  discussed  in  Section  III.  The 
introduction  to  that  process  is  given  in  Figure  33.  The  remainder  of  program 
FUZSET  is  essentially  a  question  and  answer  session,  in  which  the  user 
provides  the  necessary  level  of  damage  for  each  of  the  ten  structural 
attributes.  The  standardized  question  format  is  shown  in  Figure  34.  Note 
that  as  one  of  the  responses,  the  user  is  allowed  to  view  and/or  change  the 
fuzzy  relations  as  outlined  in  Section  IV. 

If  the  viewer  chooses  7  (view/change),  he/she  is  greeted  with  Figure  35, 
which  explains  how  the  view/change  option  works.  Upon  choosing  a  damage  level 
to  view,  the  program  retrieves  the  four  fuzzy  relationship  (one  for  each  mode) 
data  files  that  the  user  has  requested.  Sample  screens  relating  structural 
attribute  number  one  to  shear  and  flexure  are  shown  in  Figure  36.  Once  the 
user  has  finished  viewing/changing  a  given  set  of  relationships,  Figure  37 
appears,  and  the  user  has  the  option  to  view/change  fuzzy  relations  at 
different  damage  levels,  or  to  just  answer  the  original  question.  In  this 
case,  the  option  to  continue  was  picked,  and  a  value  of  severe  was  chosen  for 
structural  attribute  number  one  (Example  session  values:  4, 3, 2, 3, 2,3, 1 , 5, 1 , 2) . 

This  process  may  be  repeated  for  each  of  the  ten  attributes,  or  the  user 
may  simply  answer  the  questions  as  they  are  presented.  In  this  manner,  the 
other  nine  questions  were  answered  with  appropriate  responses  for  test  DS3. 
Upon  completion  of  all  questions,  the  program  FUZSET  is  exited,  and  the 
program  COMB  I N  is  run  in  order  to  analyze  the  information  using  the  DSW 
algorithm.  When  the  D5U  algorithm  is  completed,  the  code  returns  to  EXSYS 
where  it  asks  the  user  (Fig.  39)  if  he/she  would  like  to  view  the  results 
obtained  in  COMB  I N .  If  yes,  the  program  CQLQRPLT  is  executed,  and  the 
calculated  damage  levels  for  each  mode  are  displayed,  e.g.  see  Figure  39. 
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Comparative  examples  of  Linguistic 
terms  describing  Levels  of  Damage 

Screen  1 


NO 

DAMAGE 


SLIGHT 

DAMAGE 


MODERATE 

DAMAGE 


PRESS  AN¥  KEV 


Conparative  examples  of  Linguistic 
terms  describing  Levels  of  Damage 


SEVERE 

DAMAGE 


Screen  2 


Il:sisa5j0:ri 


n 


mm 


UERV 

EXTENSIVE 

DAMAGE 


PRESS  ANV  KEV 


TOTAL 

FAILURE 


Figure  31.  Graphical  representation  of  linguistic  terms 


>  DAMAGE  GRADE 


The  given  graph  is  the  fuzzg  set 
representation  of  the  linguistic 
tern  :  MODERATE  DAMAGE 

Enter  c  to  change  the  relation,  n  for 
the  next  damage  level,  or  e  to  exit:  ?  e 


Figure  32.  Fuzzy  set  representation  of  "Moderate  Damage." 


I  am  now  going  to  ask  you  10  questions  concerning  the  level 
of  damage  of  structure!  peremetert  related  to  structural 
Integrity  of  the  burled  box,  Vour  responses  will  be  combined 
with  expert  opinion  relating  these  same  structural  parameters 
to  modes  of  response.  This  information  on  expert  opinion  is 
stored  In  data  files  in  the  form  of  fuxxy  sets. 


Press  ENTER  to  continue  7 


Figure  33.  Introduction  to  subjective  questions. 
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The  level  or  amount  of  opal  ling  and  loss  of  concrete/rebar 
Interaction  on  the  bottom  side  of  the  main  slab  (central  1/3 
to  1/2  span)  is  needed. 

Depending  on  the  amount  or  level  of  damage  associated  with 
this  parameter,  a  fuzzy  relation  (as  determined  by  expert 
opinion)  between  the  given  parameter  and  the  different  modes 
of  deformation  will  be  added  to  a  data  file  for  later  manipul¬ 
ation.  You  may  either  choose  the  appropriate  answer  to  the 
question,  or  you  may  request  to  view  and/o"  change  the  relation¬ 
ships  as  outlined  above,  before  making  a  decision. 

CHOICES) 

1)  negligible 

2)  slight 

3)  moderate 

4)  severe 

3)  very  extensive  (concrete  core  in  chunks,  etc..) 

6)  unknown 

7)  View  and/or  change  fuzzy  relations 

I  Please  enter  the  number  of  your  choice  i  ?  7  | 


Figure  34.  Typical  format  of  structural  attribute  questions. 


The  current  structural  parameter  1st  spalling  and  loss  of  concrete/rebar 
interaction  on  the  bottom  side  of  the  main  slab  (central  1/3 
to  1/2  span) 

The  fuzzy  relation  between  the  given  structural  parameter  and 
the  different  modes  of  deformation  may  be  different  for  each 
level  of  damage.  Therefore,  before  viewing/changing  any  of 
these  relations,  you  must  choose  a  specific  level  of  damage 
to  have  displayed. 

The  choices  again  aret  1)  None 

2)  Slight 

3)  Moderate 

4)  Severe 

3)  Very  extensive 
b)  Unknown 

Which  damage  level  (vi,  given  structural  parameter)  would  you  like  to 
look  at  7  3 


Figure  35.  Preview  of  view/change  option. 
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none  1  2  3  4  5  6  7  8  9  total 

DEGREE  OF  ASSOCIAIION/RELAIION  BETWEEN  GIVEN  STRUCTURAL 
PARAMETER  A  SHEAR 

GIVEN  STRUCTURAL  PARAMETER:  spalling  and  loss  of  concrete/rebar 
interaction  on  the  bottom  side  of  the  Main  slab  (central  i/3 
to  1/2  span)  is!  MODERATE 

Enter  n  to  view  next  node  relation,  or  c  to  change  this  fuzzy  relation 
(NOTE!  If  you  change  this  relation  you  will  not  influence  similar  relations  for 
other  structural  parameters,  damage  modes,  and  damage  levels,  7ou  only  change 
the  relation  between  the  parameter,  damage  level,  and  damage  mode  shown  above)? 


DEGREE  OF  ASSOCIATION/RELATION  BETWEEN  GIVEN  STRUCTURAL 
PARAMETER  I  FLEXURE 

GIVEN  STRUCTURAL  PARAMETER!  spalling  and  loss  of  cgncrete/rebar 
interaction  on  the  bottom  side  of  the  main  slab  (central  1/3 
to  1/2  span)  is!  MODERATE 

Enter  n  to  view  next  mode  relation,  or  c  to  change  thi$  fuzzy  relation 
(NOTE:  If  you  change  this  relation  you  will  not  influence  similar  relations  for 
other  structural  parameters,  damage  modes,  and  damage  levels,  ¥qu  only  ' 
between  the  parameter,  damage  level,  and  damage  mode  s 


Figure  36.  Typical  fuzzy  set  relation  of  mode  versus 
structural  attribute. 
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You  may  now  view/changa  the  other  fussy  damage  relations  for 
the  current  structural  parameter  at  a  different  level  of  damage, 
or  continue  the- program  and  simply  choose  an  appropriate 
level  of  damage  to  answer  the  question. 

Your  choice  is  (v  =  view/change  c  «*  continue)  ?  c 


The  level  of  damage  of  spalling  and  loss  of  concr ete/rebar 
interaction  on  the  bottom  side  of  the  main  slab  (central  1/3 
to  1/2  span)  is  i 

1 )  None 

2)  Slight 

3)  Moderate 

4)  Severe 

5)  Very  extensive 

6)  Unknown 

Enter  the  number  of  your  choicei  ?  4 


Figure  37,  Exiting  fro©  the  view/change  option. 


the  user  would  like  to  see  the  results  of  the  OSW  combination  process  lust 
completed,  i.e.,  the  calculated  (fussy  set)  levels  of  damage  few  each  mode  is 

t 

1  yes  -  please  show  them 

2  no  -  continue  program  execution 
l 


Enter  number (S)  of  value(s),  WHY  for  Information  on  the  rule, 
QUIT  to  save  data  entered,  <£SC>  for  edit  mode  or  <H>  for  help 


Figure  38.  Typical  format  of  an  EXSYS  question. 
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0.8 


4  5  6  7  8 

DAMAGE  GRADE 


Category  =  Structural  Integrity 
THIS  GRAPH  REPRESENTS  THE  CALCULATED 
FUZZY  DAMAGE  LEUEL  OF  THE  STRUCTURE 
INDUCED  BY  SHEAR  RELATED  PARAMETERS . 

The  linguistic  tern  corresponding  to 
this  fuzzy  set  is  Host  nearly: 

SLIGHT  DAMAGE 

Enter  c  to  continue,  or  e  to  exit  ?  ■ 


01234  54789  10) 

- DAMAGE  GRADE 


INDUCED  BY  FLEXURE  RELATED 
PARAMETERS. 

The  linguistic  tern  corresponding  to 
this  ruzzu  set  is  Most  nearly: 
MODERATE  DAMAGE 

Enter  c  to  continue,  or  e  to  exit  ?  I 


Figure  39.  Results  of  Program  CQMBIN  as  shown  by  COLORPLT. 
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The  completion  of  COLORPLT  signals  an  end  to  the  subjective  information 
analysis,  and  the  code  returns  to  EXSYS.  Inside  EXSYS,  the  answers  obtained 
in  C0J181N  are  linguistically  interpreted  and  presented  to  the  user  in  the  form 
of  rules.  In  test  DS3,  the  dasage  sustained  in'a  shear  mode  was  slight,  and 
the  overall  damage  to  the  structure  is  severe,  as  shown  in  Figures  40  and  41 
respectively. 

The  next  step  in  the  analysis  process  is  initiated  when  EXSYS  calls  the 
external  program  TIH.  As  discussed  in  Section  10,  the  program  TIM  evaluates 
the  early-tiee  response  of  the  structure  using  a  Timoshenko  beam  model.  Input 
to  the  program  is  interactive  and  includes  the  material,  geometrical,  and 
loading  parameters  shown  in  Figure  42.  All  input  for  test  DS3  can  be  found  in 
Section  U.  The  analysis  of  the  data  in  program  TIM  takes  approximately  i 
min,  after  which  it  returns  to  EXSYS, 

The  second  step  in  the  early-tiee  analysis  is  accomplished  through  the 
use  of  program  IFPLQT.  Figure  43  shows  the  introduction  provided  to  the  user. 
After  viewing  this  introduction,  the  pressure/ impulse  versus  time  plots  for 
the  given  structure  are  presented  to  the  user  as  shown  in  Figures  44  and  43, 
respectively.  In  Figure  46,  the  user  is  Questioned  as  to  the  similarity  of 
the  plots  just  viewed.  As  options,  the  user  is  allowed  to  look  at  example 
piots  (Figs.  9  and  10)  representing  case  1  and  case  2,  or  to  return  to  the 
plots  to  study,  them  more  carefully  before  answering.  As  can  be  seen  in 
Figures  44  and  43,  the  best  answer  for  test  DS3  is  3,  l.e,,  no  significant 
similarity. 

Following  the  early-time  analysis,  EXSYS  begins  the  late-time  analysis  by 
executing  program  DEFLECT.  The  introduction  to  this  program  is  given  in 
Figure  47.  In  Figure  40  the  user  is  presentod  with  the  deflection  profiles  of 

the  damaged  structure  at  3  ms  and  IS  me.  After  viewing  this  plot,  the  user  is 

asked  three  questions  pertaining  to  the  deflection  profile,  one  of  which  is 

shown  in  Figure  49.  Although  the  answers  to  these  questions  would  at  first 

seem  obvious,  in  the  event  a  value  is  borderline  between  two  linguistic 
classes,  the  intuitive  opinion  of  the  user  is  quite  useful. 

Once  DEFLECT  is  completed,  control  is  returned  to  EXSY3.  Since  DEFLECT 
is  the  final  external  program  to  be  executed,  the  only  task  remaining  is  for 
EXSYS  to  analyze  the  information  obtained  from  the  hard  data  external 
functions.  This  task,  as  discussed  in  Section  IV,  is  accomplished  through  the 
use  of  several  rules.  As  an  example,  Figures  SO  and  31  show  rules  64  and  78 
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NUM3ER:  14 


(1)  complete  fatlun  of  the  structure  is  false  (i.e.,  the  roof  slab  has 
not  been  completely  separated  from  one  or  both  of  its  supports) 
and  (2)  shear  damage  is  slight 

THENt 

the  damage  to  the  structure  associated  with  shear  response  as 
determined  by  structural  integrity  analysis  from  visual 
information  is  slight  -  Probabi 1 i ty“l 

NOTE)  The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 

algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels  of 
the  10  structural  integrity  parameters  supplied  by  the  user  are  combined 
by  weighting  them  with  expert  knowledge  relating  modes  of  failure  to 

these  same  parameters. 


Figure  41.  Typical  rule  explaining  overall  damage 


I  -?<!i  now  going  to  evaluate  the  early  time  resoonse  of  the  structure 
usinq  a  modified  Timoshenko  beam  analysis.  Before  I  can  do  this,  however, 

1  will  need  the  values  tor  the,  following  material,  geometrical,  and  loading 
properties: 

Beam  density 
Shear  modulus 
Slab  length 
Slab  thickness 

Shear  deformation  coefficient 
F'oisson  ratio 
Damping  coefficient 
Concrete  compresssive  strength 
Steel  yield  stress 
Reinforcing  ratio 
Maximum  slab  overpressure 
Load  rise  time 
Load  duration  time 

What  is  the  value  of  the  beam  density  Clb  s**2  /  in**43 
, 000223 

Is  this  value  correct  <y  or  n>  ? 


Figure  42.  Introduction  to  Program  TIM. 


In  a  moment  you  will  be  shown  a  pressure  vs.  time  plot 
and  an  Impulse  vs.  time  plot  of  a  test  specified  by  you. 

The  Interface  pressure  and  correspond  ng  impulse  at  the 
center  line,  near  the  support,  and  over  the  support  will  be 
superimposed  on  each  other.  Th*  purpose  of  this  is  to 
determine  the  similarity  of  the  response  at  different 
locations.  The  similarity  of  pressure/impulse  at  different 
locations  is  helpful  in  determining  the  initial  response 
mode. 

After  viewing  the  pressure  vs.  time  plot  and  the 
impulse  vs.  time  plots  you  will  t>«  asked  a  question 
pertaining  to  the  above  mentioned  similarity.  Vou  may  either 
answer  the  question  promptly,  or  view  the  plots  again  in 
order  to  study  them  more  carefully  before  producing  an 
answer. 

STRIKE  ANV  KEV  TO  CONTINUE 


Figure  43.  Introduction  to  Program  IFPLOT. 
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Figure  44.  Interface  pressure  versus  time  at  three 
locations . 
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QUESTION: 

In  your  opinion,  which  of  the  cases  below  best  describes  the  given  plots  ? 

Case  1  -  Pressure/impulse  plots  at:  location  1  closely  resembles 

location  3,  AND  location  2  is  significantly  different 
than  both  1  and  3. 

Case  2  -  Pressure/Impulse  plots  at:  location  2  closely  resembles 

location  3,  AND  location  1  is  significantly  different 
than  both  2  and  3. 

Case  3  -  The  plots  suggest  that  the  similarity  of  pressure/ impulse 

is  somewhat  between  conditions  1  and  2  above  (i.e.,  no 
significant  similarity  exists  as  stated  above). 

NOTE:  As  this  is  a  very  subjective  question,  you  may 

view  two  example  plots  that  in  the  authors  opinion, 
represent  cases  1  and  2.  This  may  aid  you  in  making  your 
decision.  If  you  would  like  to  view  these  examples 
choose  'e'  below. 

YOUR  CHOICE  (1,2,3,  e  for  example,  or  c  to  view  the  plots  again):  3 


Figure  46.  IFPLOT  question  concerning  interface  pressure 
plot  similarity. 


In  a  moment,  you  will  be  shown  the  deflection  profiles  of  the 
underside  of  the  roof  slab  at  3  and  IS  msec  after  initial  loading. 
You  will  be  asked  for  the  particular  test  number  you’re  assessing. 
Along  with  the  deflection  profiles,  you  will  be  asked  a  few 
questions  pertaining  to  this  plot.  As  you  answer  the  questions, 
you  may  want  to  view  the  deflection  profiles  first}  therefore, 
you  may  use  the  ENTER  key  to  toggle  back  and  forth  between  the 
questions  and  the  plot. 


Press  ENTER  to  continue? 


Figure  47.  Introduction  to  Program  DEFLECT. 
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At  13  msec,  the  ratio  of  differential  deflection  (near  wall  deflection 
minus  centerline  deflection)  to  slab  thickness  Is  .28  . 

Typically,  when  evaluating  modes  of  response  (especially  flexure 
and  tension  membrane),  a  ratio  in  the  range  ofs 

s  <  0.3  Is  considered  relatively  flat 
s  0.3  -  0.6  is  considered  moderately  curved 
s  >  0.6  Is  considered  highly  curved 


In  your  opinion,  given  this  information  and  the  deflection  profile, 
which  category  do  you  think  the  given  ratio  of  .28  belongs  i 

1. )  Relatively  flat 

2. )  Moderately  curved 

3. )  Highly  curved 


Enter  a  value  of  1-3,  or  just  ENTER  to  view  plot  ?  I 


Figure  49.  Question  2  in  Program  DEFLECT. 
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RULE  NUMBER i  64 
IF: 


and  (2) 


and  (3) 


and  (4) 


complete  failure  of  the  structure  is  false  (i.#.,  the  roof  slab  has 
not  been  completely  separated  from  one  or  both  of  its  supports) 
early  time  response  as  determined  by  a  combination  of  Timoshenko 
beam  analysis  and  Interface  pressure  data  falls  into  the  vague  area 
of  f lexure-shear  «  1 

the  near  Mall  deflection  as  determined  by  the  external  program 
OEFLECT.BAS  is  moderate  or  large 

the  ratio  of  differential  deflection  (near  Mall  deflection  minus 
centerline  deflection)  to  slab  thickness  implies  a  relatively  flat 
surface 


late  time  response  as  determined  by  deflection  profiles  was 
predominantly  shear  (even  though  early  response  did  not  shoM 
overwhelming  evidence  to  this)  with  some  residual  flexure 
Pr ob abi llty“l 

(LATE  RESPONSE  1  IS  GIVEN  THE  VALUE  1 


IF  line  #  for  derivation,  <K>-known  data,  <C>-cholcee 
or  -  prey,  or  next  rule,  <J>-jump,  <H>-help  or  <ENTER>  to  continue: 


Figure  50,  Typical  rule  interpreting  late-time  response. 


RULE  NUMBER!  78 
IF: 


and  t2> 


complete  failure  of  the  structure  is  false  U.e.,  the  roof  slab  has 
not  besn  completely  separated  from  one  or  both  of  its  suooorts> 

C DAMAGE  LEVEL!  >  1.0 


the  overall  damage  to  the  structure  as  determined  by  deflection 
information  is  very  extansive  -  Probabi l l ty»l 

NOTE:  the  value  of  the  variable  (DAMAGE  LEVEL!  is  determined  in  the  external 

program  DEFLECT. FOR  and  returned  via  the  program  RETURN. FOR.  The  value  of 
the  variable  is  computed  as  the  Integral  of  the  deflection  (over  the 
lenQth  of  the  slab)  at  IS  mssc.  divided  by  length  times  thickness  (to 
produce  a  dimensionless  parameter). 


IF  line  S  for  derivation,  <K>-knOwn  data,  <C>-cholcet 
or  -  prev.  or  next  rule,  <J>-ju*p,  <H>-help  or  < ENTER)  to  continue: 


Figure  51.  Typical  rule  interpreting  overall  damage  to  the 
structure  as  determined  by  the  hard  data. 
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which  describe  the  nodes  of  response  and  overall  damage  of  the  structure, 
respectively.  Finally,  all  inportant  information  derived  throughout  the 
session  is  given  as  output  at  the  end  of  the  run  as  shown  in  Figure  52.  At 
this  point,  the  user  may  either  end  the  session  or  use  the  menu  at  the  bottom 
of  the  page  to  explore  the  answers.  By  using  the  facilities  within  EXSYS,  the 
user  is  able  to  back-chain  through  the  rules  in  order  to  see  how  a  particular 
piece  of  information  was  derived. 
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Values  based-  on  0/1  system 


VALUE 


1  the  external  program  FUZSET.BAS  has  been  run  1 

2  the  external  program  COMBIN.FOR  has  been  run  1 

Z  the  external  program  PLOT. BAS  has  been  run  1 

4  the  damage  to  the  structure  associated  with  shear  response  as 

determined  by  structural  integrity  analysis  from  visual 
information  is  slight  1 

5  the  damage  to  the  structure  associated  with  diagonal  tension 

response  as  determined  by  structural  integrity  analysis  from 
visual  information  is  slight  1 

6  the  damage  to  the  structure  associated  with  a  flexure  response 

as  determined  by  structural  Integrity  analysis  from  visual 
information  is  moderate  1 

V  the  damage  to  the  structure  associated  with  tensile  membrane 
response  as  determined  by  structural  integrity  analysis  from 
visual  information  is  severe  i 

8  the  overall  level  of  damage  to  the  structure  as  determined  by 
structural  Integrity  analysis  from  visual  information  is 

severe  1 

9  the  external  program  TIM. FOR  has  been  run  1 

10  the  external  program  1FPL0T.BAS  has  been  run  1 


Press  any  key  for  morel 


Values  based  on  0/1  system  value 

11  early  time  response  as  determined  by  a  combination  of 
Timoshenko  beam  analysis  snd  Interface  pressure  data  falla 

Into  the  vague  area  of  f lexure*shear  l 

12  the  external  program  DEFLECT. has  nas  been  run  1 

13  lete  time  response  es  determined  by  deflection  profiles  wee 

predominantly  shear  (even  though  eerly  response  old  not  Show 
overwhelming  evidence  to  this)  with  some  residual  flexure  1 

14  the  overall  damage  to  the  structure  as  determined  by 

deflection  Information  la  very  extensive  I 


Alt  choleee  <A>,  only  If  valueil  <6>,  Print  <P>,  Change  and  rerun  <C>* 
rulee  used  time  number >,  Qult/eeve  <0>>  Help  <M>,  Done  <0>s 


Figure  52.  Output  of  DAPS  for  Test  DS3. 
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VI.  REVIEW,  CONCLUSIONS,  AND  RECOMMENDATIONS 


REVIEW  ‘ 

Structural  damage  assessment  of  buried  protective  structures  is  a  very 
complex  subject  which  is  imbued  with  a  large  amount  of  uncertainty  and 
vaqueness.  This  is  due  to  the  Tact  that  much  of  the  information  used  in  the 
analysis  process  is  derived  from  expert  opinion  and  uncertain  numerical  data. 
Because  of  this,  a  new  approach  for  combining  this  type  of  information  was 
studied.  Of  particular  interest  was  the  feasibility  of  incorporating  both 
soft  data  in  the  form  of  expert  opinion  and  hard  data  produced  from 
instrumentation  waveforms  into  a  structural  damage  assessment  code. 

The  data  base  used  in  this  study  is  derived  from  a  series  of  experimental 
tests  conducted  in  1981  and  1982  on  eleven  buried  reinforced  concrete  boxes. 
The  boxes  were  subjected  to  extremely  high  impulsive  blast  pressures  and  thus 
sustained  levels  of  damage  ranging  between  slight  degradation  and  complete 
collapse.  Major  modes  of  response  observed  in  the  test  structures  included 
direct  shear  at  the  supports,  flexure  of  the  main  slab,  and  tension  membrane 
at  both  supports  and  midspan  induced  by  large  deflections.  Evidence  of  all 
these  modes  of  response  were  obtained  from  experimental  data  comprised  of 
interface  pressure  data,  strain  data,  and  deflection  information  obtained  from 
high-  speed  photography. 

The  framework  developed  in  this  report  is  incorporated  in  a  rule-based 
expert  system  approach.  In  an  expert  system  scheme,  difficult  problems  are 
subdivided  into  smaller  problems,  which  in  turn  are  represented  in  antecedent- 
consequent  pairs  as  rules.  These  rules  are  combined  with  othfc'  data  and 
information  to  form  what  is  called  the  knowledge  base.  The  processing  and 
analysis  of  this  information  is  controlled  through  an  inference  mechanism, 
which  retrieves  necessary  information  from  the  user  or  the  knowledge  base 
using  either  backward  or  forward  chaining. 

In  this  study,  an  analysis  of  structural  integrity  of  the  buried  box 
element  was  accomplished  using  the  expert  system  shell  EXSVS  as  the  control 
mechanism  to  combine  subjective  and  objective  information  and  as  a  mechanism 
to  chain  through  the  rule  base.  Expert  opinion  on  damage  assessment  obtained 
from  questionnaires  was  used  as  the  basis  for  the  subjective  portion  of  the 
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code.  Because  this  type  of  information  is  vague  and  uncertain,  fuzzy  set 
theory  Mas  used  to  quantify  linguistic  variables.  Also,  because  of  the  size 
and  complexity  of  the  problem,  a  numerical  method  in  the  form  of  a  fuzzy- 
Mei ghted-average  algorithm  Mas  used  instead  of  rules  to  synthesize  the 
subjective  information. 

Objective  data,  on  the  other  hand,  Mas  obtained  through  the  use  of 
external  programs  Mhich  “reach  out“  and  “hook”  numerical  information  from 
digitized  Maveforms.  This  information  is  then  passed  back  to  the  rule-based 
inference  module  Mhere  it  is  analyzed  or  interpreted  through  the  use  of  rules. 

CONCLUSIONS 

The  basic  purpose  of  this  report  mas  to  study  the  feasibility  of  using  an 
expert  system  approach  in  the  area  of  survivability  and  vulnerability  analyses 
of  buried  protective  structures.  As  such,  several  interesting  conclusions  nay 
be  dramn  from  this  mork. 

First,  the  assessment  of  damage  to  protective  facilities  is  a  highly 
complex  problem  requiring  innovative  analysis  techniques  for  three  reasons. 
One,  due  to  the  extreme  nature  of  the  short-lived  transient  loading  condition, 
many  factors  affect  the  dynaeic  response  of  the  structure.  These  factors 
include  loading  parameters  such  as  peak  pressure  and  rise  tlee,  material 
properties  like  concrete  strength  and  soil  type,  and  geometrical  properties 
like  1 ength-to-oepth  ratio,  end  restraint,  and  depth  of  soil  overburden. 
Because  of  the  interaction  of  these  various  paraaeters,  modes  of  response 
Mhich  art  unique  to  this  problea  (direct  shear  for  example)  can  be  induced  in 
the  structure.  Tmo,  damage  assessaent  is  a  fundamentally  subjective  concept 
and  thus  requires  the  use  of  linguistic  interpretation  and  engineering 
judgment  to  analyte  the  problem.  This  adds  to  the  complexity  of  the  situation 
because  both  linguistic  information  and  engineering  judgment  are  inherently 
difficult  to  quantify.  Three,  an  overall  assessaent  of  the  structure  cannot 
rely  on  structural  integrity  analysis  aione  but  eust  also  include  the  concepts 
of  functionality  and  repairability  for  coepleteness.  Analysis  is  complicated 
further  Mhen  it  is  alto  realized  that  each  of  these  concepts  may  be 
interrelated. 

Second,  although  the  information  obtained  from  expert  opinion  can  be 
quite  useful,  the  decomposition  of  this  information  is  no  trivial  task. 
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Problems  encountered  include]  interpreting  linguistic  teres;  uncertainty  in 
expert  opinion;  deciphering  different  terminology  used  to  describe  the  sane 
concept;  and  subdividing  and  organizing  the  information  into  a  meaningful 
structure,.  Also,  if  the  problem  (to  which  the  information  will  be  applied)  is 
sufficiently  complex,  difficulties  such  as  the  combinatorial  explosion 
described  in  Section  IV  may  arise.  In  this  particular  case,  the  problem  was 
circumvented  through  the  use  of  a  fuzzy-weighted-average  algorithm,  rather 
than  the  conventional  rule-based  format. 

Third,  although  the  term  "objective  data"  was  used  in  this  work  to 
describe  numerical  information  obtained  from  instrumentation  waveforms,  it  was 
found  that  these  data  do  contain  some  inherent  subjectivity.  The  subjectivity 
arises  when  we  try  to  use  discrete  numerical  quantities  derived  from 
continuous  functions.  For  example,  discrete  values  of  deflection  are  used  in 
Section  IV  to  describe  various  nodes  of  damage.  Deflection  of  the  roof  slab 
at  any  given  point,  however,  is  a  function  of  both  time  and  position. 
Therefore,  tfr  choice  of  a  discrete  deflection  value  is  a  subjective  decision. 
The  approach  used  in  this  report  was  to  use  engineering  judqment  to  classify 
the  various  numerical  Quantities  into  different  linguistic  levels,  then 
proceed  with  the  analysis  using  these  linguistic  values. 

Fourth,  a  careful  analysis  of  the  problem  and  a  review  of  the  literature 
reveals  that  typical  analytical  methods  (deflection  ratios,  dissipated 
hysteretic  energy,  etc.)  employed  in  describing  levels  of  damage,  have  limited 
use  in  the  case  under  study.  Calculation  of  different  damage  measures 
revealed  that,  at  the  present  time,  not  enough  is  known  about  late-time 
effects  on  structural  response  produced  by  changes  in  loading,  material,  and 
geometrical  parameters.  For  example,  a  simple  damage  measure  using  energy  and 
impulse  was  calculated  at  two  points  along  the  roof  span.  Although  the 
measure  predicted  damage  well  within  a  group  containing  similar  geometrical 
and  material  parameters,  it  was  not  able  to  predict  daeage  when  these 
parameters  were  changed.  Therefore,  without  an  understanding  of  parameter 
effects  on  response,  developing  a  good  measure  of  structural  damage  is 
difficult. 

Fifth,  the  theory  of  fuzzy  sets  can  be  quite  effective  in  interpreting 
linguistic  information  that  is  mither  vague,  imprecise,  or  uncertain.  The 
assignment  of  fuzzy  sets  to  linguistic  tens  and  the  combination  of  these 
quantities  is  a  straightforward  procedure.  The  major  difficulty  in 
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confronting  this  situation  is  that  of  providing  linguistic  teres  for 
calculated  fuzzy  sets.  Equation  6  was  used  in  this  study  to  provide  a 
relative  eeasure  of  siailarity  between  sets.  In  the  event  that  a  calculated 
fuzzy  set  falls  between  two  predefined  fuzzy  sets,  the  equation  siaply  assigns 
the  linguistic  tera  corresponding  to  the  aore  siailar  of  the  two.  In 
practice,  however,  a  better  approach  would  be  to  assign  a  linguistic  tera 
which  is  a  coaproaise  of  the  two.  Thus,  although  Equation  &  is  adequate  for 
deaonstration  purposes,  other  techniques  could  be  developed  which  would  better 
interpret  results  obtained  froa  fuzzy  operations. 

Sixth,  although  quite  useful  for  soae  parts  of  a  problea,  a  rule-based 
foraat  is  not  feasible  in  certain  cases  where  the  i nterrel ati onshi p  of  a  large 
nuaber  of  variables  aakes  the  problea  too  coaplex.  Such  a  problea  was 
described  in  Section  IV  concerning  the  experts'  opinion.  However,  the  use  of 
expert  systee  techniques  allows  the  integration  of  such  situations  into  the 
problea.  In  this  case,  the  problea  was  circuavented  through  the  use  of  an 
external  proqraa  which  aanipulated  and  contained  the  data.  Similarly,  other 
types  of  strategies  can  be  caployed  in  the  solution  of  a  problea.  Besides 
those  discussed  in  this  report,  other  tools  such  as  finite  eieaent  codes  or 
pattern  recognition  techniques  (Refs.  8,9)  could  also  be  used. 

RECOMMENDATIONS 

Based  on  the  inforaation  learned  froa  this  study,  the  following 
recoaaendations  are  presented  for  possible  future  work. 

1.  Because  of  tiae  constraints,  it  was  not  possible  to  test  the  code  on 
structures  other  than  those  in  the  N£$  dynaaic  test  series.  Therefore,  the 
code  should  be  tested  against  another  test  series  in  order  to  deteraint  its 
strengths  and  weaknesses. 

2.  As  was  shown  in  Figure  17,  the  work  undertaken  in  this  report 
coapnses  only  the  structural  integrity  aodule.  In  order  to  complete  the 
code,  aodule*  on  functionality  and  rc^.rability  should  also  be  developed. 

3.  The  only  objective  data  used  in  this  study  included  interface 
pressures  and  deflection  profiles.  As  discussed  in  Section  III,  other 
objective  data  are  available  and  should  be  studied  for  possible  inclusion  into 
the  hard  data  saction  in  order  to  aake  it  aore  coeplete.  The  possibility  of 
using  pattern  recognition  procedures  on  the  data  should  be  explored. 


97 


4.  The  expert  aatrix  [El  described  in  Section  IV  was  developed  using  the 
authors’  opinions  as  its  basis.  A  final  version  of  the  eatrix,  however, 
should  be  coaposed  of  the  opinions  of  several  experts.  The  inforaation  used 
to  fill  this  aatrix  was  collected  through  the  use  of  a  second  questionnaire 
sent  to  seven  experts,  but  tiae  constraints  did  not  persit  its  inclusion  in 
this  report.  Therefore,  analysis  of  the  inforaation  and  integration  into  the 
code  should  be  accoaplished. 

5.  A  detailed  paraaetric  study  should  be  conducted  to  deteraine  the 
effects  of  different  paraaeters  (load  intensity,  concrete  strength,  length-to- 
depth  ratio,  and  reinforcing  ratio)  on  the  late-tiae  response  of  the 
structure.  An  early-tiae  analysis  capability  currently  exits  in  DAPS.  With 
the  incorporation  of  a  late-tiae  analysis  to  coapleaent  the  existing  early- 
tiae  analysis,  a  study  of  suitable  daaage  level  aeasures  can  be  undertaken. 

6.  More  work  needs  to  be  done  in  the  area  of  siailarity  and  difference 
aeasures.  Specifically,  better  techniques  are  needed  for  converting 
calculated  fuzzy  sets  into  aeaningful  linguistic  interpretations. 
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APPENDIX  A 

PHOTOGRAPHS  OF  DAMAGED  STRUCTURES 
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1C2 


Figure  A-l.  Posttest  photograph  of  Test  Element  DS1. 


Figure  A-2.  Posttest  photograph  of  Test  Element  DS2 . 


Figure  A- 3 .  Posttest  photograph  of  Test  Element  DS3 . 


Figure  A-4.  Posttest  photograph  of  Test  Element  DS4. 


Figure  A-5.  Posttest  photograph  of  Test  Element  DS5. 
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Figure  A-6 .  Posttest  photograph  of  Test  Element  DS5. 


Figure  A-10.  Posttest  photograph  of  Test  Elo 
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-Figure  A-ll. 


DS2-4 


:test  photograph  of  Test  Element  DS2-4. 
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Figure  A-12.  Posttest  photograph  of  Test  Element  DS2-5. 


113 


APPENDIX  B 

EXPLANATION  OF  EXSYS  FEATURES 


EXSYS  is  a  generalized  expert  system  development  package.  The  package 
runs  on  IBM  PC  or  compatible-type  computers  with  256K  RAM  or  greater.  EXSYS 
can  create  about  700  rules,  with  an  average  of  6  or  7  conditions,  per  64K  of 
memory  over  192K.  This  is  roughly  5000  rules  in  a  PC  with  640K  RAM.  Some  of 
the  more  important  features  of  EXSYS  are  detailed  here. 

EDITOR 

Expert  systems  are  generated  in  EXSYS  using  a  rule  editor  program  called 
EDITXS.  This  editor  enables  you  to  easily  edit  rules  along  with  providing  a 
simple  and  rapid  way  of  putting  rules  into  the  computer  and  testing  the  rules 
as  they  are  generated.  Rules  are  input  into  the  knowledge  base  through  the 
use  of  several  menus  and  templates.  All  linguistic  input  is  in  the  form  of 
normal  human  language,  and  thus  no  computer  languages  or  special  procedures 
need  be  learned  before  using  it.  This  feature  makes  the  use  of  expert  system 
shells  very  appealing  to  both  novices  and  experienced  computer  programmers. 

t  The  format  of  a  rule  entered  into  EXSYS  is  as  follows: 


IF: 

Conditions 

THEN: 

Conditions 

and 

Choices 

ELSE: 

Conditions 

and 

Choices 

NOTE:  • 

REFERENCES:  . 

"Condition"  -  A  condition  is  simply  a  statement  of  fact.  In  EXSYS,  there  are 
two  main  types  of  conditions:  text  and  mathematical.  A  text  condition  is  a 
sen!;  a  that  may  be  true  or  false.  The  condition  is  made  up  of  two  parts,  a 
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Qualifier  and  one  or  more  Values.  The  Qualifier  is  usually  the  part  of  the 
condition  up  to  and  including  the  verb.  The  Values  are  the  possible 
completions  of  the  sentence  started  by  the  Qualifier. 

"Choices"  -  Choices  are  all  the  possible  solutions  to  the  problem  among  which 
the  expert  system  will  decide.  The  goal  of  EXSYS  is  to  select  the  most  likely 
choice  based  on  input  data  or  to  provide  a  list  of  possible  choices  arranged 
in  order  of  likelihood. 

CALCULATIQNftL  SYSTEMS 

EXSYS  has  three  built  in  cal cul ational  systems  for  handling  uncertainty 
in  Qualifier  and  Choice  decisions. 

0  or  1  system  -  this  system  should  be  used  when  the  choice  or  qualifier  is 
perfectly  true  or  perfectly  false.  There  is  no  ranking  of  acceptable  choices 
in  this  system,  it  is  a  simple  yes  or  no  situation. 

0-10  system  -  this  is  the  most  generally  useful  system,  especially  if 
probabilistic  rules  are  used  (e.g.,  if  a  condition  is  true,  there  is  an  80'/. 
likelihood  that  choice  X  is  appropriate).  Values  of  0  and  10  allow  you  to 
completely  eliminate  (0)  or  definitely  include  (10)  a  choice.  The  remaining 
values  of  1-9  are  averaged  over  the  rules  used. 

-100  to  +100  -  this  system  is  similar  to  the  0-10  system,  except  it  allows 
data  to  be  combined  as  independent  or  dependent  probabilities.  If  the  rule 
base  is  such  that  precise  statistical  data  are  available,  then  use  this 
system. 

MATHEMATICAL  EXPRESSIONS  AND  VARIABLES 

While  many  expert  systems  can  be  developed  without  the  use  of  numeric 
variables,  the  ability  to  calculate  numeric  variables  is  a  powerful  tool. 

Even  if  mathematical  calculations  are  not  needed  in  the  expert  system,  numeric 
variables  can  be  used  to  flag  text  for  display  at  the  end  of  a  run. 

Mathematical  expressions  containing  variables  can  be  used  to  derive 
information  in  the  IF  and  THEN  portions  of  the  rule  just  like  any  other 
qualifier.  Variables  may  be  either  string  or  numeric.  Mathematical 
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expressions  involving  numeric  variables  can  be  combined  using  the  basic 
algebraic  operations  I),  trigonometric  operations,  exponentiation,  and 

others. 

EXTERNAL  PROGRAMS 

External  programs  can  be  called  from  within  EXSYS  for  data  acquisition 
and  calculation,  and  then  can  be  passed  back  to  EXSYS  for  analysis.  EXSYS  can 
directly  receive  data  from  automatic  testing  equipment,  data  bases,  spread 
sheets,  and  dedicated  programs.  The  process  involved  in  calling  external 
programs  to  get  the  value  for  a  specific  qualifiers  or  numeric  variables  is 
quite  simple.  Each  qualifier  or  variable  has  an  explanatory  test  associated 
with  it.  In  order  to  run  an  external  program,  it  is  necessary  only  to  start 
the  text  associated  with  the  qualifier/variable  with: 

RUN (f i 1  enamel 

where  “filename"  is  the  name  of  the  external  program  to  run.  Data  are  passed 
to  external  programs  via  the  file  PASS.DAT  and  returned  via  the  file 
RETURN.DAT.  Therefore,  the  program  must  only  be  able  to  read  and  write  to 
these  disk  files. 

INFERENCE  MECHANISM 

The  different  types  of  inference  mechanisms  (forward  and  backward 
chaining)  were  described  in  Section  IV.  A  newly  added  feature  of  EXSYS  that 
makes  it  very  powerful  is  the  ability  to  use  ei-ther  forward  or  backward 
chaining  or  combinations  of  both.  The  choice  of  inferencing  schemes  is  chosen 
before  running  the  expert  system  as  part  of  the  command  line  options. 

EXPLANATION  FEATURES 

At  the  end  of  an  EXSYS  run,  or  when  a  rule  is  displayed,  or  information 
is  asked  of  the  user,  the  user  has  the  option  of  asking  the  computer  one  of 
several  things: 

1.  The  computer  may  display  the  rule  or  rules  that  allowed  it  to  derive 
the  information.  A  rule  used  for  derivation  will  have  information  about  the 
condition  the  user  is  asking  for  in  its  THEN  part.  The  user  may  then  continue 
asking  how  the  computer  knew  that  the  rule’s  IF  conditions  were  true  and  so 
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on.  If  asked  about  a  condition  that  is  an  alQebraic  expression,  the  values  of 
the  variables  in  the  expression  will  be  displayed.  The  user  may  then  ask  how 
these  values  were  derived  by  entering  the  number  of  the  variable. 

2.  The  computer  may  respond  that  you  told  it  the  information  was  true. 

3.  The  computer  may  respond  that  it  does  not  yet  know  if  the  condition 
is  true. 

4.  If  the  information  cane  from  an  external  program,  EXSYS  will  give  the 
program  name  from  which  it  came. 

The  features  outlined  here  are  but  a  few  of  the  capabilities  of  EXSYS.  For 
further  information,  see  the  owner’s  manual,  or  contact  EXSYS  at:  EXSYS,  Inc., 
P.Q.  Box  7S1S8,  Contr.  Sta.  14,  Albuquerque,  NM  87194,  (505)  836-6676. 
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APPENDIX  C 
RULES 

Subject:  Structural  damage  assessment  of  reinforced  concrete  buried  boxes 
subjected  to  impulsive  blast  loadings. 

Author:  Steve  J.  Savage 

Starting  text: 

The  follOMing  program  is  an  initial  attempt  at  combining  linguistic, 
information  in  the  form  of  expert  opinion  on  damage  assessment  with  crisp 
numerical  data  obtained  from  instrumentation  waveforms.  This  program  is  not 
an  end  product  as  such,  but  rather  a  demonstration  tool  of  what  can  possibly 
be  achieved  using  expert  system  techniques  in  attacking  such  a  problem.  The 
program  makes  several  calls  to  external  programs  which  take  from  a  few  seconds 
up  to  several  minutes,  so  please  be  patient  with  the  system. 

The  first  portion  of  the  code  deals  with  the  structural  integrity 
analysis  of  the  structure  through  the  use  of  visual  information}  therefore, 
the  user  will  need  to  have  observed  or  inspected  the  structure  or  have  visual 
information  in  the  form  of  photographic  data  available. 

The  second  portion  of  the  code  deals  with  structural  integrity  analysis 
of  the  structure  via  numerical  data.  Host  of  the  hard  data  have  been  stored 
internally  in  the  form  of  data  files,  but  some  information  must  be  derived 
from  the  user.  Please  have  the  following  information  available  relating  to 
the  material,  geometrical,  and  loading  parameters  of  the  structure  :  beam 
density,  shear  modulus,  roof  slab  clear  span  length,  slab  thickness,  Poisson’s 
ratio,  damping  coefficient,  concrete  compressive  strength,  steel  yield  stress, 
percentage  steel  ratio,  maximum  slab  overpressure,  load  rise  time,  and  load 
duration. 

Ending  text: 

The  information  to  follow  is  the  output  from  the  program  DAPS.  At  the 
present  time,  the  only  information  to  be  output  is  the  damage  level  to  the 
structure  and  most  likely  modes  of  deformation  as  determined  by:  1)  structural 
integrity  analysis  via  soft  data  ,  and  2)  structural  integrity  analysis  via 
hard  data.  All  applicable  rules  are  used  in  data  derivations. 
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RULES! 


RULE  NUHBER:  1 
IFi 

complete  Failure  of  the  structure  is  true  (i.e.,  the  root  slab  has  been 
completely  separated  From  one  or  both ‘oF  its  supports) 


THEN: 

[FAILURE  MODES!  IS  GIVEN  THE  VALUE  “NOT  KNOWN" 


RULE  NUMBER:  2 
IF: 

complete  Failure  oF  the  structure  is  true  (i.e.,  the  rooF  slab  has  been 
completely  separated  From  one  or  both  oF  its  supports) 
and  the  Failure  has  occurred  at  the  slab-wall  connection 
and  the  separated  slab  is  lying  Flat  on  the  Floor  (implying  a  symmetric 
Fai lure) 

and  the  Failure  surFace(s)  at  both  supports  are  relatively  “clean"  with  a 
single  surFace  Failure 

and  the  inclination  oF  the  Failure  surface(s)  at  both  sides  are  relatively 
vertical 

and  inspection  oF  the  Failure  region  indicates  the  main  reinForcing  bars 
have  been  mostly  severed  or  sheared  oFF 

THEN] 

damage  level  as  determined  by  structural  integrity  analysis  From  visual 
inFormation  is  total  Failure  -  Probability*! 
and  the  models!  involved  in  the  deFormation  process  as  determined  by  visual 
(subjective)  inFormation  was  predominantly  direct  shear  -  Probability^ 
and  [FAILURE  NODES!  IS  GIVEN  THE  VALUE  "KNOWN" 


RULE  NUMBER]  3 
IFi 

complete  Failure  oF  the  structure  is  true  (i.e.,  the  rooF  slab  has  been 
completely  separated  From  one  or  both  oF  its  supports) 
and  the  Failure  has  occurred  at  the  slab-wall  connection 
and  the  separated  slab  is  lying  Flat  on  the  Floor  (implying  a  symmetric 
Fai lure) 

and  the  Failure  surFace(s)  at  both  supports  are  relatively  “clean"  with  a 
single  surFace  Failure 

and  the  inclination  oF  the  Failure  surFace(s)  at  both  sides  are  relatively 
inclined  or  unclear 

and  inspection  oF  the  Failure  region  indicates  the  main  reinForcing  bars 
exhibited  rupture  aFter  significant  deformation  (note  "necking"  or 
stretching) 
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THEN: 


damage  level  as  determined  by  structural  integrity  analysis  from  visual 
information  is  total  failure  -  Probabi  1  i.ty=l 
and  the  mode(s)  involved  in  the  deformation  process  as  determined  by  visual 
(subjective)  information  was  predominantly  diagonal  tension  - 
Probabi 1 ity=l 

and  [FAILURE  MODES 3  IS  GIVEN  THE  VALUE  "KNOWN" 


RULE  NUMBER:  4 
IF: 

complete  failure  of  the  structure  is  true  (i.e.,  the  roof  slab  has  been 
completely  separated  from  one  or  both  of  its  supports) 
and  the  failure  has  occurred  at  the  slab-wall  connection 
and  the  separated  slab  is  lying  flat  on  the  floor  (implying  a  symmetric 
failure) 

and  the  failure  surface(s)  at  both  supports  are  relatively  rough  with  many 
cracks  and  “concrete  teeth* 

and  the  inclination  of  the  failure  surface(s)  at  both  sides  are  relatively 
inclined  or  unclear 

and  inspection  of  the  failure  region  indicates  concrete  crushed  in  the 
compression  zone  with  possibly  ripped  out  or  protruding  rebar 


THEN: 

damage  level  as  determined  by  structural  integrity  analysis  from  visual 
information  is  total  (ailure  -  Probabi 1 i t y *  1 
and  the  mode(s)  involved  in  the  deformation  process  as  determined  by  visual 
(subjective)  information  was  predominantly  shear-compression  - 
Probabi 1 i t  y« 1 

and  [FAILURE  NODES!  IS  GIVEN  THE  VALUE  “KNOWN* 


RULE  NUMBER:  3 
IF: 

complete  failure  of  the  structure  is  true  (i.e.,  the  roof  slab  has  been 
completely  separated  from  one  or  both  of  its  supports) 
and  the  failure  has  occurred  at  the  slab-wall  connection 
and  the  separated  slab  is  lying  flat  on  the  floor  (implying  a  symmetric 
f al lure) 

and  the  failure  surface(s)  at  both  supports  are  relatively  rough  with  many 
cracks  and  “concrete  teeth* 

and  inspection  of  the  failure  region  indicates  the  main  reinforcing  bars 
exhibited  rupture  after  significant  deformation  (note  “necking*  or 
stretching! 


THEN: 

damage  level  as  determined  by  structural  integrity  analysis  from  visual 
i  formation  is  total  failure  *  Probabillty«l 
and  the  models)  involved  in  the  deformation  process  as  determined  by  visual 
(subjective)  information  were  predominantly  shear  and  diagonal  tension, 
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causing  rupture  of  reinforcement  after  significant  rotation  and 
’  deformation.  It  is  possible  that  the  failure  region  was  underreinforced 
-  Probabi 1 i ty=l 

and  C  FA  I  LURE  MODES!  IS  GIVEN  THE  VALUE  "KNOWN” 


RULE  NUMBER:  6 
IF: 

conplete  failure  of  the  structure  is  true  (i.e.,  the  roof  slab  has  been 
completely  separated  from  one  or  both  of  its  supports) 
and  the  failure  has  occurred  at  the  slab-wall  connection 
and  the  separated  slab  is  lying  flat  an  the  floor  (implying  a  symmetric 
f ai lure) 

and  the  failure  surface(s)  at  both  supports  are  relatively  rough  with  many 
cracks  and  “concrete  teeth* 

and  inspection  of  the  failure  region  indicates  concrete  crushed  in  the 
compression  zone  with  possibly  ripped  out  or  protruding  rebar 
and  the  inclination  of  the  failure  surface(s)  at  NOT  both  sides  are 
relatively  inclined  or  unclear 


THEN: 

damage  level  as  determined  by  structural  integrity  analysis  from  visual 
information  is  total  failure  -  Probability=l 
and  the  mode(s)  involved  in  the  deformation  process  as  determined  by  visual 
(subjective)  information  were  predominantly  shear  and  diagonal  tension, 
causing  crushing  of  concrete  after  significant  rotation.  It  is  possible 
that  the  failure  region  was  overreinforced  or  the  concrete  was  weak  in 
strength  -  Pr obabi 1 i ty« 1 
and  (FAILURE  MODES!  IS  GIVEN  THE  VALUE  “KNOWN* 


RULE  NUMBER:  7 

IF: 

complete  failure  of  the  structure  is  true  (i.e.,  the  roof  slab  has  been 
completely  separated  from  one  or  both  of  its  supports) 

and  the  failure  has  occurred  well  in  the  slab 

and  the  separated  slab  it  lying  flat  on  the  floor  (implying  a  symmetric 
f ai lure) 

and  the  failure  surface(s)  at  both  supports  are  relatively  "clean"  with  a 
tingle  surface  failure 

and  the  inclination  of  the  failure  surface(s)  at  both  sides  are  relatively 
inclined  or  unclear 

and  inspection  of  the  failure  region  indicates  the  main  reinforcing  bars 
exhibited  rupture  after  significant  deformation  (note  "necking"  or 
stretching) 


THEN: 

damage  level  as  determined  by  structural  integrity  analysis  from  visual 
information  is  total  failure  -  Probabi 1 i ty*l 
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and  the  mode(s)  involved  in  the  deformation  process  as  determined  by  visual 
(subjective)  information  was  predominantly  punching  shear.  This  -failure 
mode  is  slower -than  direct  shear,  where  failure  surfaces  are  vertical  - 
Probability=l 

and  [FAILURE  MODES]  IS  GIVEN  THE  VALUE  "KNOWN" 


RULE  NUMBER:  8 
IF: 

complete  failure  of  the  structure  is  true  (i.e.,  the  roof  slab  has  been 
completely  separated  from  one  or  both  of  its  supports) 
and  the  failure  has  occurred  at  the  slab-wall  connection 
and  the  separated  slab  is  leaning  on/towards  one  side  of  the  structure 
(implying  an  unsymmetrical  failure) 

and  the  failure  surface(s)  at  one  support  is  defined  by  a  single  surface 
failure  which  is  relatively  "clean",  and  the  other  is  still  partially 
attached  or  relatively  "rough"  with  many  cracks  and  “concrete  teeth" 
and  the  inclination  of  the  failure  surface(s)  at  one  side  is  relatively 
vertical  and  the  other  side  is  relatively  inclined  or  unclear 


THEN: 

damage  level  as  determined  by  structural  integrity  analysis  from  visual 
information  is  total  failure  -  Probabi 1 i ty® l 
and  the  models)  involved  in  the  deformation  process  as  determined  by  visual 
(subjective)  information  was  predominantly  direct  shear,  followed  by 
concrete  crushing  and  rebar  pullout  of  the  other  side  (possibly  due  to 
extremely  high  loading  and/or  weak  concrete)  -  Probability®! 
and  [FAILURE  MODES]  IS  GIVEN  THE  VALUE  "KNOWN" 


RULE  NUMBER:  9 
IF: 

[FAILURE  MOOES!  =>  "NOT  KNOWN" 

and  complete  failure  of  the  structure  is  true  (i.e.,  the  roof  slab  has  been 
completely  separated  from  one  or  both  of  its  supports: 


THEN: 

damage  level  as  determined  by  structural  integrity  analysis  from  visual 
information  is  total  failure  -  Probab i 1 i ty® l 
and  the  models)  involved  in  the  deformation  process  as  determined  by  visual 
(subjective)  information  is  unknown  due  to  insufficient  or  inadeOuate 
data  -  Probability8! 


NOTE: 

The  information  provided  from  the  visual  data  was  either  inadequate  or 
inconsistent  with  information  contained  in  previous  rules  to  determine  a 
mode(s)  of  failure  with  any  considerable  amount  of  certainty. 
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RULE  NUMBER:  10 


IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 


THEN: 

RUNIC: \DAPS\HAIN\FUZSET) 

and  the  external  program  FUZSET.BAS  has  been  run  -  Probability3! 

NOTE: 

The  THEN  portion  of  this  rule  simply  calls  an  external  program  which  asks 
questions  of  the  user  pertaining  to  the  structural  integrity  of  the 
structure,  as  ascertained  from  visual  information. 


RULE  NUHBER:  11 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 


THEN: 

RUN (C: \DAPS\MAIN\CQHBIN) 

and  the  external  program  CORBIN. FOR  has  been  run  -  Probabi 1 i ty®l 
NOTE: 

The  external  program  CORBIN. FOR  uses  fully  logic  in  the  fora  of  the  OSH 
(Dong,  Shah,  Mong)  algorithm  to 'combine  the  linguistic  information 
obtained  from  the  user,  with  expert  opinion  relating  this  information  to 
various  modes  of  failure.  The  output  is  a  group  of  4  futsy  sets  relating 
the  level  of  damage  associated  with  a  given  mode. 


RULE  NURBER:  12 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  tus  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  the  user  would  like  to  see  the  results  of  the  OSH  combination  process 
just  completed,  i.e.,  the  calculated  (funy  set)  levels  of  damage  for 
each  aodt  is  :  yes  -  please  show  them 


THEN: 

run (C:\0APS\RA1N\C0L0RPLT) 

and  the  external  program  PLOT. BAS  has  been  run  -  Probability3! 

NOTE: 

Program  C0L0RPIT. BAS  takes  the  results  of  the  OSH  combination  process 
obtained  from  the  program  CORBIN. FOR  and  simply  plots  the  resulting  fu*:v 
sets  (damage  level  of  a  given  mode). 
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•  RULE  NUMBER:  13 


IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  shear  damage  is  none 

THEN: 

the  damage  to  the  structure  associated  with  a  shear  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
none  (i.e.,  no  appreciable  daBage)  -  Probabi 1 i t y  =  l 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  14 
IF: 

complete  failure  of  the  structure  is  false  ti.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  shear  damage  is  slight 

THEN: 

the  damage  to  the  structure  associated  with  shear  response  as  determined 
by  structural  integrity  analysis  from  visual  information  is  siignt  - 
Probability®! 


NOTE- 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  IS 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supoorts* 
and  shear  damage  is  moderate 
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THEN: 


the  damage  to  the  structure  associated  with  shear  response  as  determined 
by  structural  integrity  analysis  free  visual  indorsation  is  moderate  - 
Probabi 1 ity=l 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  duzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
od  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  od 
dailure  to  these  same  parameters. 


RULE  NUMBER:  16 
IF: 

complete  dailure  od  the  structure  is  dalse  (i.e.,  the  rood  slab  has  not 
been  completely  separated  drom  one  or  both  od  its  supports! 
and  shear  damage  is  severe 

THEN: 

the  damage  to  the  structure  associated  with  shear  response  as  determined 
by  structural  integrity  analysis  drom  visual  indormation  is  severe  - 
Probabi 1 ity«l 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
od  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  od 
dailure  to  these  same  parameters. 


RULE  NUMBER:  17 
IF: 

complete  dailure  od  the  structure  is  dalse  (i.e.,  the  rood  slab  has  not 
been  completely  separated  drom  one  or  both  od  its  supports) 
and  shear  damage  is  very  extensive 

THEN: 

the  damage  to  the  structure  associated  with  shear  response  as  determined 
by  structural  integrity  analysis  drom  visual  information  is  very 
extensive  -  Probability*! 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  duzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
od  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  od 
dailure  to  these  same  parameters. 
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RULE  NUMBER:  18 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  froa  one  or  both  of  its  supports! 
and  shear  damage  is  unknown 

THEN: 

the  damage  to  the  structure  associated  with  shear  response  as  determined 
by  structural  integrity  analysis  from  visual  information  is  unknown  - 
Probabi ii ty=l 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  19 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  diagona:  tension  damage  is  none 

THEN: 

the  damage  to  the  structure  associated  with  diagonal  tension  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
none  (i.e., no  appreciable  damage)  -  Probab i 1 i ty= 1 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  .fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  20 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  diagonal  tension  damage  is  slight 
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TKENi 


the  damage  to  the  structure  associated  with  diagonal  tension  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
slight  -  Probability5! 


MOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by.  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  21 


IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  diagonal  tension  damage  is  moderate 

THEN: 

the  damage  to  the  structure  associated  with  diagonal  tension  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
moderate  -  Probabi 1 i ty=l 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  f uzzy-weightea  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  bv  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  22 
IF: 

complete  failure  of  the  structure  is  false  (i.e,,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  ot  its  supports) 
and  diagonal  tension  damage  is  severe 

THEN: 

the  damage  to  the  structure  associated  with  diagonal  tension  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
severe  -  Probabi 1 i ty= 1 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 
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RULE  NUMBER:  23 


IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not. 
been  completely  separated  from  one  or  both  of  its  supports) 
and  diagonal  tension  damage  is  very  extensive 

THEN: 

the  damage  to  the  structure  associated  with  diagonal  tension  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
very  extensive  -  Prebabi 1 ity=l 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  24 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  diagonal  tension  damage  is  unknown 

THEN: 

the  damage  to  the  structure  associated  with  diagonal  tension  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
unknown  -  Probability*! 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to-  these  same  parameters. 


RULE  NUMBER:  23 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  net 
been  completely  separated  from  one  or  both  of  its  supports) 
and  flexure  damage  is  none 


THEN: 


the  damage  to  the  structure  associated  with  a  flexure  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
none  (i.e.,  no  appreciable  damage)  -  Probabi 1 i t y  =  1 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm'  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  26 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  flexure  damage  is  slight 

THEN: 

the  damage  to  the  structure  associated  with  a  flexure  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
slight  -  Probability*! 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  27 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  flexure  damage  is  moderate 

THEN: 

the  damage  to  the  structure  associated  with  a  flexure  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
moderate  -  Probabi 1 ity=t 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 
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RULE  NUMBER:  28 


IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  flexure  damage  is  severe 

THEN: 

the  damage  to  the  structure  associated  with  a  flexure  response  as 
determined  by  structural  integrity,  analysis  from  visual  information  is 
severe  -  Probabil ity=l 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  29 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  flexure  damage  is  very  extensive 

THEN: 

the  damage  to  the  structure  associated  with  a  flexure  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
very  extensive  -  Probability®! 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  30 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  flexure  damage  is  unknown 
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THEN: 


the  damage  to  the  structure  associated  with  a  flexure  response  as 
determined  by  structural  integrity  .analysis  from  visual  information  is 
unknown  -  Probability5! 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  31 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roo.f  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports)  and  tension 
membrane  damage  is  none 


THEN: 

the  damage  to  the  structure  associated  with  tensile  membrane  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
none  (i.e.,  no  appreciable  damage)  -  Pr obabi 1 i ty=  1 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  32 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  tension  membrane  damage  is  slight 

THEN: 

the  damage  to  the  structure  associated  with  tensile  membrane  response  as 
determined  by  structural  integrity  analysis  from  visual  Information  is 
slight  -  Probability*! 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 
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RULE  NUHBERs  33 


IF  i 

complete  Failure  oF  the  structure  is  False  (i.e.,  the  rooF  slab  has  not 
been  completely  separated  From  one  or  both  oF  its  supports) 
and  tension  membrane  damage  is  moderate 

THEN! 

the  damage  to  the  structure  associated  with  tensile  ir.*iiibr  ane  response  as 
determined  by  structural  integrity  analysis  From  visual  inFormation  is 
moderate  -  Probabi 1 ity=l 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  Fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
oF  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledg*  relating  modes  oF 
Failure  to  these  same  parameters. 


RULE  NUMBER!  34 
IF! 

complete  Failure  oF  the  structure  is  False  (i.e.,  the  rooF  slab  has  not 
been  completely  separated  From  one  or  both  oF  its  supports) 
and  tension  membrane  damage  is  severe 

THEN! 

the  damage  to  the  structure  associated  with  tensile  membrane  response  as 
determined  by  structural  integrity  analysis  From  visual  inFormation  is 
severe  -  Probabi 1 i ty  =  1 


NOTE! 

The  actual  analysis  procedure  stated  above  is  a  Fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
oF  the  10  structural  Integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  uith  expert  knowledge  relating  modes  oF  - 
Failure  to  these  same  parameters. 


RULE  NUMBER:  35 
IFi 

complete  Fallire  oF  the  structure  is  False  (i.e.,  the  rooF  slab  has  not 
been  completely  separated  From  one  or  both  oF  its  supports) 
and  tension  membrane  damage  is  very  extensive 
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THEN: 


the  damage  to  the  structure  associated  with  tensile  Membrane  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
very  extensive  -  Probabi 1 ity=l 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  Zb 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  tension  membrane  damage  is  unknown 

THEN: 

the  damage  to  the  structure  associated  with  tensile  membrane  response  as 
determined  by  structural  integrity  analysis  from  visual  information  is 
unknown  -  Probability*! 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  37 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  overall  damage  is  none 

THEN: 

the  overall  level  of  damage  to  the  structure  as  determined  by  structural 
Integrity  analysis  from  visual  information  is  none  (i.e.,  no  appreciable 
damage)  -  Probability*! 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  mooes  of 
failure  to  these  same  parameters.  The  overall  damage  is  assumed  to  be 
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the  largest  damage  level  from  each  of  the  4  modes  -  shear,  flexure, 
diagonal  tension,  tension  membrane. 


RULE  NUMBER:  38 
IF: 

complete  failure  of  the  structure  is  false  (i.e,,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  overall  damage  is  slight 

THEN: 

the  overall  level  of  damage  to  the  structure  as  determined  by  structural 
integrity  analysis  from  visual  information  is  slight  -  Probabi 1 ity=l 


NOTE: 

The  actual  analysis  procedure  stated  above  is  a  f uzzy-wei ghted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 

The  overall  damage  is  assumed  to  be  the  largest  damage  level  from  each 
of  the  4  modes  -  shear,  flexure,  diagonal  tension,  tension  membrane. 


RULE  NUMBER:  39 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  Its  supports) 
and  overall  damage  is  moderate 

THEN: 

the  overall  level  of  damage  to  the  structure  at  determined  by  structural 
integrity  analysis  from  visual  information  is  moderate  -  Probability*! 


NOTE: 

The  overall  damage  is  assumed  to  be  the  largest  damage  level  from  each  of 
the  4  modes  -  shear,  flexure,  diagonal  tension,  and  tension  membrane. 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 


RULE  NUMBER:  40 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
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and  overall  daaage  is  severe 


THEN: 

the  overall  level  of  daaage  to  the  structure  as  determined  by  structural 
integrity  analysis  from  visual  information  is  severe  -  Probability-: 

NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program,  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 

The  overall  damage  is  assumed  to  be  the  largest  damage  level  from  each 
t^f  the  4  modes  -  shear,  flexure,  diagonal  tension,  tension  membrane. 


RULE  NUMBER:  41 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  overall  damage  is  very  extensive 

THEN: 

the  overall  level  of  damage  to  the  structure  as  determined  by  structural 
integrity  analysis  from  visual  information  is  very  extensive  - 
Probability*! 

NOTE: 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters. 

The  overall  damage  is  assumed  to  be  the  largest  damage  level  from  each 
of  the  4  modes  -  shear,  flexure,  diagonal  tension,  tension  membrane. 


RULE  NUMBER!  42 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  tha  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  overall  damage  is  unknown 

THEN: 

the  overall  level  of  daaage  to  the  structure  as  determined  by  structural 
integrity  analysis  from  visual  information  is  unknown  -  Probability5! 

NOTE* 

The  actual  analysis  procedure  stated  above  is  a  fuzzy-weighted  average 
algorithm  implemented  in  a  FORTRAN  computer  program.  The  damage  levels 
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of  the  10  structural  integrity  parameters  supplied  by  the  user  are 
combined  by  weighting  them  with  expert  knowledge  relating  modes  of 
failure  to  these  same  parameters.  . 

The  overall  damage  is  assumed  to  be  the  largest  damage  level  from  each 
of  the  4  modes  -  shear,  flexure,  diagonal  tension,  tension  membrane. 


RULE  NUMBER!  43 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 


THEN: 

RUN(C; \DAPS\NAIN\TIH) 

and  the  external  program  TIM. FOR  has  been  run  -  Probabi 1 i t y  =  1 
NOTE! 

Program  TIM. FOR  is  used  to  calculate  the  early  time  response  of  the 
structure  (i.e.,  at  less  than  2  ms)  by  using  a  modified  Timoshenko  beam 
analysis.  Although  the  program  is  highly  complex,  the  output  is  simple; 
early  time  response  is  categorizes  as  either  shear  or  flexure  depending 
on  which  value  reaches  its  critical  failure  level  first. 


RULE  NUHBERi  44 
IF: 

complete  failure  of  the  structure  is  false  (i.e,,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 


THEN! 

RUN (Ct \DAPS\KA1N\ IFPLQT) 

and  the  external  program  IFPLQT. BAS  has  been  run  -  Probability*! 

NQTEi 

Program  IFPL0T.8AS  is  used  to  help  determine  the  early  time  response  of 
the  structure  (i.e.,  at  less  than  2  ms)  by  making  a  similarity  comparison 
of  pressure  versus  time  plots  at  three  locations  of  the  roof  slab.  The 
three  locations  are  near  the  wall,  at  the  slab  centerline,  and  over  the 
wall.  Similarity  of  different  combinations  of  interface  pressures  at  the 
roof/soil  interface  help  determine  modes  of  response. 


RULE  NUMBER!  45 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  tiee  response  as  detereined  by  modified  Timoshenko  beam  analysis  is 
shear 
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and  early  tine  response  as  detemined  by  similarity  analysis  of  interface 
pressure  plots  is  shear 


THEN: 

early  tine  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  shear  - 
Probability*! 


RULE  NUMBER:  4 b 
IF: 

complete  failure  of  the  structure  is  false  li.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports! 
and  early  time  response  as  determined  by  modified  Timoshenko  beam  analysis  is 
flexure 

and  early  time  response  as  determined  by  similarity  analysis  of  interface 
pressure  plots  is  flexure 

THEN: 

early  time  response  as  determined  by  a  combination  of  Timosnenko  Dean 
analysis  and  interface  pressure  data  is  predominantly  flexure  - 
Probabi 1 ity°l 


RULE  NUMBER:  47 
IF; 

complete  failure  of  the  structure  is  false  U.e.,  the  ro«»f  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  modified  Timoshenko  beam  analysis  is 
shear 

and  early  time  response  as  detereined  by  similarity  analysis  of  interface 
pressure  plots  is  flexure  or  flexure-shear 


THEN: 

early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  falls  into  the  vague  area  of 
flexure-shear  -  Probabil ity»l 


RULE  NUMBER:  48 
IF: 

complete  failure  of  the  structure  is  false  ii.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  modified  Timoshenko  beam  analysis  * 
flexure 

and  early  time  response  as  determined  by  similarity  analysis  of  interface 
pressure  plots  is  f lexure-shear 
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THEN: 


early  tine  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  falls  into  the  vague  area  of 
f lexure-shear  -  Probabi 1 i ty=l 


RULE  NUMBER:  49 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 


THEN: 

RUNIC: \DAPS\MA I N\ DEFLECT) 

and  the  external  program  DEFLECT. BAS  has  been  run  -  Probabi 1 i tyal 
NOTE: 

Program  DEFLECT. BAS  shows  the  user  the  roof  slab  deflection  profiles  at  3 
and  IS  ms,  and  asks  three  questions  of  the  user  that  are  used  in 
successive  rules  to  help  determine  late  time  modes  of  response. 


RULE  NUMBER;  50 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  shear  *  l 
and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  very  small  or  small 

and  the  ratio  of  differential  deflection  (near  wail  deflection  minus 
centerline  deflection)  to  slab  thickness  implies  a  relatively  flat 
surface 


THENt 

late  time  response  as  determined  by  deflection  profiles  was  predominantly 
shear  -  Probability"! 
and  (LATE  RESPONSE!  15  GIVEN  THE  VALUE  l 


RULE  NUMBER:  51 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  r.ot 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  shear  a  i 
and  the  near  wall  deflection  as  deterained  by  the  external  program 
DEFLECT. BAS  is  very  small  or  small 
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and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  moderately  curved 
surface 


THEN: 

late  time  response  as  determined  by  deflection  profiles  was  initially 
shear,  becoming  substantial  flexure  as  time  progressed  -  Probability=l 
and  CLATE  RESPONSE]  IS  GIVEN  THE  VALUE  1 


RULE  NUMBER:  52 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  thp  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  shear  =  1 
and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  small 

and  the  ratio  of  differential  deflection  (near  mall  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  highly  curved  surface 

THEN: 

late  time  response  as  determined  by  deflection  profiles  mas  initially 
shear,  but  residual  strength  forced  it  into  a  tension  membrane  mode, 
possibly  due  to  weak  concrete  -  Probabi l i ty»l 
and  CLATE  RESPONSE!  IS  GIVEN  THE  VALUE  1 


RULE  NUMBER:  53 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports! 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  apd  interface  pressure  data  is  predominantly  shear  *  i 
and  the  near  mil  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  moderate 

and  the  ratio  of  differentia!  deflection  (near  nail  deflection  minus 
centerline  deflection)  to  slab  thickness  implies  a  relatively  fl^t 
surface 


THEN: 

late  time  response  as  determined  by  deflection  profiles  mas  predoainantiy 
shear  -  Probabi 1 i ty=l 
and  (LATE  RESPONSE]  IS  GIVEN  THE  VALUE  l 
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RULE  NUMBER:  54 


IF: 

complete  failure  of  the  structure  is  *a!se  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  shear  =  1 
and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  moderate 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  moderately  curved 
surface 


THEN: 

late  time  response  as  determined  by  deflection  profiles  was  a  combination 
of  shear  and  flexure  -  Probabi 1 i ty=l 
arid  CLATE  RESPONSE!  IS  GIVEN  THE  VALUE  1 


RULE  NUMBER:  55 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  shear  =  1 
and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  moderate 

ar.d  the  ratio  of  differential  deflection  (near  wall  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  highly  curved  surface 


THEN: 

late  time  response  as  determined  by  deflection  profiles  was  substantial 
shear,  followed  by  membrane  response  throughout  -  Probabi 1 ity=l 
and  CLATE  RESPONSE!  IS  GIVEN  THE  VALUE  1 


RULE  NUMBER:  56 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  shear  =  1 
and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  large 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 
centerline  deflection)  to  slab  thickness  implies  a  relatively  flat 
surface 


THEN: 


late  time  response  as  determined  by  deflection  profiles  was  almost 
exclusively  shear  -  Frobabi 1 i ty-1 
and  [LATE  RESPONSE!  IS  GIVEN  THE  VALUE  1 


RULE  NUMBER:  57 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  shear  =  1 
and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  large 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  moderately  curved 
surface 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  moderately  curved 
surface 


THEN: 

late  time  response  as  determined  by  deflection  profiles  was  mostly  shear 
with  some  residual  flexure  response  -  Probabi 1 i ty=l 
and  CLATE  RESPONSE!  IS  GIVEN  THE  VALUE  1 


RULE  NUMBER:  58 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  shear  =  1 
and  che  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  Urge 

and  the  ratio  of  differential  deflection  (near  wail  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  highly  curved  surface 


THEN: 

late  time  response  as  determined  by  deflection  profiles  was  a  combination 
of  shear  and  tension  memurane  response  (possibly  due  to  weak  concrete 
mix)  -  Probabi 1 i ty® l 

and  (LATE  RESPONSES  IS  GIVEN  THE  VALUE  1 
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RULE  NUMBER:  59 


IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  flexure  =  1 
and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  very  small 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 
centerline  deflection)  to  slab  thickness  implies  a  relatively  flat 
surface  or  a  moderately  curved  surface 


THEN: 

late  time  response  as  determined  by  deflection  profiles  was  predominantly 
flexure  -  Probability=l 
and  CLATE  RESPONSE!  IS  GIVEN  THE  VALUE  1 


RULE  NUMBER:  60 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  flexure  =  1 
and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  very  small 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  highly  curved  surface 


THEN: 

Ute  time  response  as  determined  by  deflection  profiles  is  unknown 
because  of  inconsistent  information  and/or  non-intuitive  combination  of 
antecedents  ~  Probabii  ityM 
and  CLATE  RESPONSE!  IS  GIVEN  THE  VALUE  i 


RULE  NUMBER:  61 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  bot:  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  falls  into  the  vague  area  of 
flexure-shear  »  l 

and  the  near  wait  deflection  as  determined  by  the  external  program 
DEFLECT, BAS  is  very  small  or  small 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 
centerline  deflection)  to  slab  thickness  implies  a  relatively  flat 
surface 
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THEN: 


late  tine  response  as  determined  by  deflection  profiles  was  a  combination 
of  shear  and  flexure  -  Probahility=l  . 
and  [LATE  RESPONSE]  IS  GIVEN  THE  VALUE  1 


RULE  NUMBER:  62 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  bath  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  falls  into  the  vague  area  of 
f i exure-shear  =  1 

and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  small 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  moderately  curved 
surface 


THEN: 

late  time  response  as  determined  by  deflection  profiles  initiated  as 
shear  and/or  diagonal  tension,  but  flexure  dominated  thereafter  - 
Probabi  1  i  ty=l 

and  [LATE  RESPONSE!  IS  GIVEN  THE  VALUE  l 


RULE  NUMBER:  63 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  Interface  pressure  data  is  predominantly  flexure  =  1 
and  the  near  Nall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  small 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  highly  curved  surface 


THEN: 

late  time  response  as  determined  by  deflection  profiles  was  predominantly 
tension  membrane  *  Probability®! 
and  [LATE  RESPONSE]  IS  GIVEN  THE  VALUE  1 


RULE  NUMBER:  64 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
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and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  falls  into  the  vague  area  of 
flexure-shear  =1 

and  the  near  Hall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  moderate  or  large 

and  the  ratio  of  differential  deflection  (near  nail  deflection  minus 
centerline  deflection)  to  slab  thickness  implies  a  relatively  flat 
surface 


THEN: 

late  time  response  as  determined  by  deflection  profiles  was  predominantly 
shear  (even  though  early  response  did  not  show  overwhelming  evidence  of 
this)  with  some  residual  flexure  -  Probabi 1 i ty=i 
and  CLATE  RESPONSE!  IS  GIVEN  THE  VALUE  1 


RULE  NUMBER:  65 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  falls  into  the  vague  area  of 
flexure-shear  =  1 

and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  moderate 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  moderately  curved 
surface 


THEN: 

'.ate  response  as  determined  by  deflection  profiles  was  a  combination  of 
shear  and  flexure  (even  though  early  response  did  not  show  shear  to  be 
such  a  large  factor)  -  Probabi 1 i t y= 1 
and  CLATE  RESPONSE!  IS  GIVEN  THE  VALUE  1 


RULE  NUMBER:  66 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  falls  into  the  vague  area  of 
flexure-shear  =  1 

and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  moderate 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  highly  curved  surface 
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THEN! 

late  time  response  as  determined  by  deflection  profiles  was  initially 
shear  (even  though  early  .response  did  not  show  overwhelming  evidence  of 
this  fact)  followed  by  extensive  bending  into  membrane  mode  - 
Probability*! 

and  CLATE  RESPONSE]  IS  GIVEN  THE  VALUE  1 


RULE  NUMBER!  67 
IF! 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  flexure  *  1 
and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  very  small 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 
centerline  deflection)  to  slab  thickness  implies  a  relatively  flat 
surface 

and  the  difference  between  differential  deflection  ratios  at  the  centerline 
and  at  a  distance  t  *  thickness  from  wall  is  small 


THEN! 

late  time  response  as  determined  by  deflection  profiles  was  predominantly 
diagonal  tension  -  Probability*! 

and  late  time  response  as  determined  by  deflection  profiles  was  predominantly 
flexure  -  Probabi 1 ity=0 
and  CLATE  RESPONSE]  IS  GIVEN  THE  VALUE  1 

NOTE! 

The  second  THEN  statement  is  used  to  override  a  previous  rule  which  would 
have  found  flexure  to  be  dominant 


RULE  NUMBER!  68 
IFi 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  falls  into  the  vague  area  of 
flexure-shear  *  1 

and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  very  small 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 
centerline  deflection)  to  slab  thickness  implies  a  relatively  flat 
surface 

and  the  difference  between  differential  deflection  ratios  at  the  centerline 
and  at  a  distance  t  »  thickness  from  wall  is  small 


THEN: 

late  tine  response  as  determined  by  deflection  profiles  was  predominantly 
diagonal  tension  -  Probabi lity=l 

and  late  tine  response  as  deternined  by  deflection  profiles  was  predominantly 
flexure  -  Probability=0 
and  ELATE  RESPONSE!  IS  GIVEN  THE  VALUE  1 

NOTE: 

The  second  THEN  statement  is  used  to  override  a  previous  rule  which  would 
have  found  flexure  to  be  dominant. 


RULE  NUMBER:  69 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  falls  into  the  vague  area  of 
flexure-shear  =  i 

and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  very  small 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  moderately  curved 
surface 

and  the  difference  between  differential  deflection  ratios  at  the  centerline 
and  at  a  distance  t  =  thickness  from  wall  is  small 

THEN: 

late  time  response  as  determined  by  deflection  profiles  was  a  combination 
of  diagonal  tension  and  flexure  -  Probabi 1 l ty*l 
and  late  time  response  as  determined  by  deflection  profiles  was  predominantly 
flexure  -  Probabi 1 i ty=0 
and  C LATE  RESPONSE!  IS  GIVEN  THE  VALUE  l 

NOTE: 

The  second  THEN  statement  is  used  to  override  a  previous  rule  which  would 
have  found  flexure  to  be  dominant. 


RULE  NUMBER:  70 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  eerly  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  Interface  pressure  data  is  predominantly  flexure  *  1 
and  the  near  wall  deflection  as  determined  by  the  external  program 
DEFLECT. BAS  is  very  small 

and  the  ratio  of  differential  deflection  (near  wall  deflection  minus 

centerline  deflection)  to  slab  thickness  implies  a  moderately  curved 
surface 
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and  the  difference  between  differentia}  deflection  ratios  at  the  centerline 
and  at  a  distance  t  *  thickness  from  wall  is  small 


THEN} 

late  tine  response  as  determined  by  deflection  profiles  was  a  combination 
of  diagonal  tension  and  flexure  -  Probabi 1 i ty=l 
and  late  time  response  as  determined  by  deflection  profiles  was  predominantly 
flexure  -  Probability*!) 
and  ELATE  RESPONSE!  IS  GIVEN  THE  VALUE  1 

NOTE: 

The  second  THEN  statement  is  used  to  override  a  previous  rule  which  would 
have  found  flexure  to  be  dominant. 


RULE  NUMBER:  71 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  (LATE  RESPONSE!  •  0 

and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  shear  =  1 


THEN: 

late  time  response  as  determined  by  deflection  profiles  is  unknown 
because  of  inconsistent  information  and/or  non-intultive  combination  of 
antecedents  -  Probability*! 

NOTE: 

The  most  likely  late  time  response  was  probably  shear  based  on  early  time 
response. 


RULE  NUMBER:  72 
IF: 

complete  failure  of  the  structure  it  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  (LATE  RESPONSE!  *  0 

and  early  time  response  as  determined  by  a  combination  of  Timoshenko  beam 
analysis  and  interface  pressure  data  is  predominantly  flexure  *  l 


THEN: 

late  time  response  as  determined  by  deflection  profiles  is  unknown 
because  of  Inconsistent  information  and/or  non-intuiti ve  combination  of 
antecedents  -  Probabi lity*l 

NOTE: 

The  most  likely  late  time  response  was  probably  flexure  based  on 
information  from  marly  time  response. 


RULE  NUMBER:  73 


IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  [LATE  RESPONSE!  =  0 

and  early  time  response  as  determined  by  a  combination  of  Timpshenko  beam 
analysis  and  interface  pressure  data  falls  into  the  vague  area  of 
f lexure-shear  =  1 


THEN: 

late  time  response  as  determined  by  deflection  profiles  is  unknown 
because  of  inconsistent  information  and/or  non-intuitive  combination  of 
antecedents  -  Probab i 1 i t y= 1 


NOTE: 

A  best  guess  of  late  time  response  Mould  be  a  combination  of  flexure  and 
shear  based  on  early  time  response. 


RULE  NUMBER:  74 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports: 
and  [DAMAGE  LEVEL!  <  0.2 

THEN: 

the  overall  damage  to  the  structure  as  determined  by  deflection 
information  is  none  -  Probability*! 


NOTE: 

The  value  of  the  variable  (DAMAGE  LEVEL!  is  determined  in  the  external 
program  DEFLECT. FOR  and  returned  via  the  program  RETURN. FOR. 


RULE  NUMBER:  73 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  (OAMAGE  LEVEL!  >  0.20 
and  [DAMAGE  LEVEL!  <  0.40 

THEN: 

the  overall  damage  to  the  structure  as  determined  by  deflection 
information  is  slight  -  Probability®! 


NOTE: 


The  value  of  the  variable  [DAMAGE  LEVEL]  is'  determined  in  the  external 
program  DEFLECT. FOR  and  returned  via  the  program  RETURN. FOR. 


RULE  NUMBER:  7 b 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  [DAMAGE  LEVEL]  >  0.40 
and  [DAMAGE  LEVEL]  <  0.70 

THEN: 

the  overall  damage  to  the  structure  as  determined  by  deflection 
information  is  moderate  -  Probability®! 


NOTE: 

The  value  of  the  variable  [DAMAGE  LEVEL]  is  determined  in  the  external 
program  DEFLECT. FOR  and  returned  via  the  program  RETURN. FOR. 


RULE  NUMBER:  77 
IF: 

complete  failure  of  the  structure  is  false  (i.e. ,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
and  (DAMAGE  LEVEL!  >  0.7 
and  [DAMAGE  LEVEL!  <1.0 

THEN: 

the  overall  damage  to  the  structure  as  determined  by  deflection 
information  is  severe  -  Probability®! 


NOTE: 

The  value  of  the  variable  C0ANA0E  LEVEL!  is  determined  »n  the  external 
program  DEFLECT. FOR  and  returned  via  the  program  RETURN. FOR. 


RULE  NUMBER:  78 
IF: 

complete  failure  of  the  structure  is  false  (i.e.,  the  roof  slab  has  not 
been  completely  separated  from  one  or  both  of  its  supports) 
ind  (DAMAGE  LEVEL!  >  1.0 

THEN: 

the  overall  damage  to  the  structure  as  determined  by  deflection 
information  is  very  extensive  -  Probabil »ty=l 


150 


iruifL-ir.-ift: ic,r-»vinuivn.xvi'1iv xmutvg. n x.r. r.y. mvauguauavav »vx.y. v.  «v  w-j  ir-i-ir-: wv.iv  y, y. 


NOTE; 


The  value  of  the  variable  [DAMAGE  LEVEL]  i.s  determined  in  the  external 
program  DEFLECT. FOR  and  returned  via  the  program  RETURN. FOR.  The  value 
of  the  variable  is  computed  as  the  integral  of  the  deflection  (over  the 
length  of  the  slab)  at  15  ms  divided  by  length  times  thickness  (to 
produce  a  dimensionless  parameter). 
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LIST  OF  SYMBOLS 


Ai 

®2 

®3 

B4 

B5 

B6 

B' 

BP 

DOB 

d 

d» 


Eij 

EO 


G 

IF 

L 

P 

R 

Slj 

SE 

• 

t 

U3 

H 

°i 

s 

m(X) 


p 

p' 

fi* 

4 

0 


1th  structural  attribute 

fuzzy  sat  representation  of  no  damage 

fuzzy  set  representation  of  slight  damage 

fuzzy  set  representation  of  moderate  damage 

fuzzy  set  representation  of  severe  damage 

fuzzy  set  representation  of  very  extensive  damage 

fuzzy  set  representation  of  unknown  damage 

fuzzy  set  representation  of  calculated  damage  level 

blast  pressure  gage 

depth  of  burial 

deptn  of  concrete  to  principal  tensile  reinforcement 
depth  of  concrete  to  principal  compressive 
reinforcement 

element  of  expert  matrix  [E] 
strain  gage 

compressive  concrete  strength 

ultimate  steel  strength 

steel  yield  strength 

shear  modulus 

interface  pressure  gage 

roof  slab  clear  span  length 

interface  pressure 

rotational  end  restraint  of  slab 

linguistic  rating  factor 

horizontal  distance  between  shear  reinforcement 
soil  strain  gage 

distance  between  shear  reinforcement 

slab  thickness  or  time 

element  of  user  matrix  [U] 

linguistic  weighting  or  importance  factor 

ith  element  of  a  set 

value  of  1th  membership  level 

deflection  of  roof  slab 

membership  value  of  element  x 

Poisson's  ratio 

percentage  of  principal  tensile  reinforcement 
percentage  of  principal  compressive  reinforcement 
percentage  of  principal  shear  rainforcaaent 
internal  angla  of  friction 
strength  enhancement  factor 
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